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CX. A Study of Charge Exchange by Helium Ions in Metals 


By G. A. DISSANAIKE 
Cavendish Laboratory, Cambridge * 


[Received July 4, 1953] 


ABSTRACT 


A study has been made of the equilibrium fractions of He®, He+ and 
Het in a helium beam traversing metallic foils of beryllium, aluminium 
and silver at energies from 0-13 to 1-1 Mev. The equilibrium fractions 
have been found to be independent of the metals used. The ratio o,/c, 
of the cross sections for electron loss and capture respectively by the 
helium ions is proportional to v”, where v is the ion velocity. In the 
exchange He®= He* the exponent m is equal to 3-08-+-0-10in the energy 
range 0-13 to 0-6 Mev studied. In the exchange He+=Het++ m appears 
to be constant at 5-10-+0-12 for energies between 0:43 and 1-1 mey, 
but decreases steadily with decreasing energy below this region. 

For the He** ion, the speed for which o,=c, is found to be 2-55 times 
that of the orbital electron in the hydrogen atom, and is almost exactly 
twice as great as the corresponding speed for the Het ion. 


§1. IyrropuctTiIon 


A stupy of charge exchange of light ions in passing through matter 
is of interest in the understanding of the processes of electron capture 
and electron loss and of the energy loss of the ions, especially near the 
end of their range where charge exchange is most frequent. 

The subject of charge exchange has been widely studied for light 
ions and the early work has been reviewed by Richardt (1933). The 
first evidence for electron capture by high energy alpha-particles was 
obtained by Henderson (1922), and the phenomenon was studied more 
closely for energies above 0-65 mev by Rutherford (1924), Kapitza (1924), 
Henderson (1925), and Briggs (1927) using natural sources. Considerable 
work has also been carried out below 0-1 Mev and is summarized by 
Massey and Burhop (1952). Recently Snitzer (1953) has investigated 
the charge exchange for alpha-particles in different gases from 0-1 Mev 
to 0-48 mev. Similar work has been reported for hydrogen beams in 
gases by Keene (1949), Montague (1951), Ribe (1951), and Kanner (1951), 
and in metals by Hall (1949). Theoretical interpretations of these 
exchange processes have been attempted using classical mechanics by 
Fowler (1924) and Thomas (1927), and on quantum-mechanical ideas 
by Oppenheimer (1928) and Brinkman and Kramers (1930). 

OR eo ah a oe 
* Communicated by A. P, French. 
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The present work describes experiments on alpha-particles in different 
metals for energies from 0-13 mev to 1-1 Mev. For this energy range it 
has been of interest to consider the results in the light of the simple 
considerations made by Bohr (1948) in a general survey of the exchange 


phenomena. 
§ 2. DEscRIPTION OF APPARATUS 


2.1. General 
The general arrangement of the target and detector systems is shown 
schematically in figs. 1 (a) and 1 (6). A beam of Het ions, stabilized to 
Fig. 1 
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(Vertical section.) 


within +1 kev, from the Cavendish Laboratory Electrostatic Generator 
entered the target box through steel tubes SS which shielded the beam 
from the stray field of the magnetic analyser. Collimating apertures CC 
(No. 32 drill size) positioned the beam centrally on the target foil F 
The beam current was collected and measured in a Faraday cage Fo. 
Particles scattered at 90° to the incident direction (having dndersone a 
single nuclear collision in the foil) passed through aperture C, into the 
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analyser between pole pieces MM, and with the magnetic field zero entered 
the secondary electron multiplier assembly (fig. 1 (b)) which will be 
referred to as the neutral detector. The detector solid angle was limited 
by the aperture C, of diameter 3-5mm, C, merely preventing any 
scattering from the inner walls of the waveguide. The whole detector 
system was carefully aligned with a light source at the target position. 
Ten collodion foils of varying thickness arranged in two sets in the foil 
changer A enabled stopping thicknesses in the range 0-1 to 7-4 mm air 
at N.T.P. to be introduced in the path of the scattered particles. An 
auxiliary mercury diffusion pump system was used to help maintain a 
good vacuum. A Philips gauge, not shown, registered the pressure near 
the secondary electron multiplier. The pressure throughout the 
experiment did not exceed 1-5 10-> mm of mercury. 

The charged particles analysed by the magnetic field were detected 
by an.arrangement similar to the one used for the monitoring, which 
consisted of a zinc sulphide screen Z and photomultiplier tube P. A lucite 
rod L (53 cm long) was used as a light guide so that the multiplier could 
be removed from the influence of the analyser field. Measurements 
were made with two different magnetic analysers. The first, a 90° 
focussing analyser (Burcham and Freeman 1949), was used in the study 
of the He® content and of the He+ and He*++ groups for energies up to 
0-6 Mev, and is the one generally referred to in this description. A 60° 
focussing analyser was used for the comparison of the intensities 
of the He+ and Het+* groups between 0-6 Mev and 1-1 Mev. The analysers 
were calibrated in terms of magnet current for energies up to about 
3:0 Mev in the way described by Burcham and Freeman (1949). The 
calibration of the 60° analyser was extended to 4-4 Mev using the particle 
groups from the reactions *D(dp)?T, C(dp)#®C, °Be(da«’)’Li* and 
*Be(da)’Li. These calibrations were checked during the experiment 
using the alpha-particles from 1°F(p«’)'*O* and the elastically scattered 
protons from thin films of beryllium, silver and aluminium. By deducing 
the alpha-particle energies in the experiment from the positions of the 
Het group (which should give exactly four times the required energy 
value) it was possible to reduce the error in the energy measurement to 
within --5 kev for the energies studied. 


2.2. Target Foils 


For the comparison of the Het and Het contents, thin self-supporting 
foils of beryllium, aluminium and silver of thickness between 0:3 and 
0-7 mm air equivalent at N.T.P. were used. The spread in energy of the 
scattered particles as they emerged from the foils ranged from 50 to 
90 kev. For the study of the He® content, however, where the energy 
distribution of the particles could not be analysed as for the charged 
components, it was necessary to reduce this energy spread considerably 
in order to define the energy within a narrow interval. The metal was 
therefore evaporated on thin collodion film ; when the beam was directed 
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into the backing, the collodion burnt off leaving a film of the metal. 
Silver and aluminium foils of between 0-15 mm and 0-3 mm air equivalent 
at N.T.P. were obtained in this way, producing an energy spread of less 
than 30 kev in the scattered particles. In using such thin foils it was 
important to consider the effect due to the scattering from the carbon left 
by the burnt collodion and also gradually built up during the bombard- 
ment; the particles scattered from carbon have energies smaller than 
those from the metal itself and would therefore give rise to a higher He® 
content. This effect was found to be serious in the case of aluminium 
and the final measurements were made using a self-supporting foil of 
0-2 mm air equivalent at N.T.P. prepared with special care. The foils 
withstood quite satisfactorily the beam currents (up to 50 mya concen- 
trated on an area 0:06 cm?) used in the experiment. 


2.3. Detectors 
(a) Charged Particles ; 

A zine sulphide screen and an E.M.I. photomultiplier (type 5032) were 
used for the analyser detector and for the monitoring system. In the 
latter a thin aluminium foil (~0-2 mm air at N.T.P.) was placed in front 
of the ZnS screen to prevent any light from the target box from entering 
the detector. The efficiency of detection of particles below about 0-6 Mev 
was well below 100° and depended on the energy of the particles. In 
studying the charged components with the magnetic analyser, however, 
the efficiency of detection would be the same for a particular energy in 
both He+ and Het+ groups provided the detecting system were not 
influenced by the stray field of the analyser. The use of the light guide 
ensured this condition and did not appreciably reduce the size of the 
pulses in the detector. To ensure stability a high-tension battery unit 
was used to supply the voltage direct to each dynode of the photomulti- 
plier. 

(b) Neutral Particles 


A secondary-electron-multiplier was used for the study of the He® 
content. Successful use of these counters for the detection of particles 
has been reported by many workers and recently reviewed by Allen 
(1950). The efficiency of detection of alpha-particles and protons above 
about 5 kev has been found to be near 100°/, and independent of the 
energy of the particles. The multiplier used in the present work was of 
the design developed in the Cavendish Laboratory by Baxter (1949). 
The dynodes were of beryllium-copper alloy activated by a thin layer of 
beryllium oxide. The multiplier was encased in steel to ensure adequate 
magnetic shielding. In the experiment about 300-350 volts per stage 
were used on the dynodes. 

The use of the secondary-electron-multiplier to measure the He® content 
required a detailed knowledge of its response to both charged and neutral 
particles and to photons. The full discussion of this problem is rather 
lengthy and is therefore given in the separate section that follows. 
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§3. THe PERFORMANCE OF THE SECONDARY ELECTRON MULTIPLIER 


In the experimental arrangement (fig. 1 (b)) the first dynode surface 
‘looked’ directly at the target and was found to be sensitive to the 
x-rays given off during the bombardment. The pulse-height distributions 
for alpha-particles and for the background radiation are however 
essentially different and it was found possible to bias out the background 
satisfactorily. As an example, fig. 2 (obtained using the Hutchinson— 
Scarrott kicksorter, 1951) shows differential bias curves due to scattered 
alpha-particles of energy 0-6 Mev and the associated background radiation ; 
curve A is obtained without the magnetic field on which deflects the 
charged particles away from the detector, and therefore represents the 
total contribution of particles and background ; curve B is obtained with 


Fig. 2 


Kicksorter counts (arbitrary units) 


is) 15 2s 35 


Bias in volts. 


Differential bias curves for 0-6 Mev alpha-particles and associated background 
radiation. Curve A shows the total contribution due to the particles 
plus background, curve B due to the background alone. 


the magnetic field on and with a stopping foil to cut off the undeflected 
neutral particles, and therefore gives the background ‘contribution alone. 
The relative intensity of this background contribution was found to 
increase with increasing energy of bombardment and was thus more 
difficult to bias out at high energies. In the experiment a bias of 28 volts 
ensured that the pulsés counted were entirely due to the alpha-particles 
for the energies below 0-6 Mev that were studied. 

The use of the detector to compare the intensities of the neutral and 
charged components of the helium beam presupposes the same efficiency 
of detection for all charge states at the same energy. This is a reasonable 
assumption to make since the mean free path for electron capture by the 
alpha-particle is much less than 10-* mm air equivalent for all energies 
used (Rutherford 1924). An experimental proof of the constant efficiency 
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was possible. From the curves of fig. 2 the pulse-height distribution due 
to charged and uncharged alpha-particles together may be obtained. 
At the same energy, the differential bias curve obtained with the magnetic 
field on, but with no stopping foil to cut off the neutrals, gives the pulse- 
height distribution due to the neutral particles alone once the x-ray 
background is subtracted. It was found that the two distributions so 
obtained could be made to coincide quite satisfactorily by multiplying 
the ordinates of one of them by a suitable constant factor, which clearly 
implies a similar response to charged and uncharged particles by the 
detector. 


$4. EXPERIMENTAL PROCEDURE 


The experiment was carried out in two stages. Firstly, the fraction of 
the He® content in the total scattered beam was obtained using the neutral 
detector. The beam was monitored by observation of the scattered 
particles at 90° using the scintillation counter. With the magnet field 
zero, the detector counted all scattered particles, charged and uncharged. 
With the magnet current increased from zero, the counting rate decreased 
suddenly to a constant value which was thereafter independent of the 
magnet current ; this confirmed the fact that the detecting systems were 
not influenced by the stray field of the analyser. The counting rates were 
then observed alternately with and without the magnetic field in several 
repeated runs at a particular energy. In the measurements a stopping 
1oil of appropriate thickness was interposed in front of the detector to 
cut off any particles scattered from the carbon on the targets. 

The intensities of the Het and He** groups analysed by the magnet. 
spectrometer were then compared. A plot of the counting rate against the 
magnet current gave the position and shape of the two peaks. A resolution 
of 2 to 4% in energy was used in the analyser system, and since the energy 
spread of the detected particles was much greater than the energy interval 
accepted by the analyser, the ratio of the heights of the peaks gave 
directly the proportion of Het to He*+ ions in the beam. Several runs 
were made over the top of each peak at a particular energy. The beam 
was monitored in one of two ways. With the 60° focussing analyser, the 
beam current was measured using a current amplifier and integrator 
arrangement; with the 90° focussing analyser, it was found more conven- 
ient to monitor on the scattered neutral particles detected by the second- 
ary electron multiplier. 

In the case of beryllium, two groups of scattered particles, from the 
metal and from oxygen impurity, were observed. The oxygen-peak was 
nearly twice as high as the beryllium-peak. The equilibrium ratio 
obtained here corresponded to a stopping material of beryllium oxide 
rather than to the metal itself. For energies below about 0-5 Mev the 
contribution due to the * tail’ of the He+ group scattered from beryllium 
was appreciable at the position of the He++ group from oxygen, and a 
subtraction procedure was necessary to separate them. 
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§5. ANALYSIS OF EXPERIMENTAL RESULTS 


5.1. Conditions for Charge Equilibrium 

The fractions of He®, He+ and He++ present in the beam can be related 
to the cross sections for capture and loss of electrons. 

Let fo(x), f(w) and f,(x) represent the fractions of He®, Het and He++ 
ions respectively in the beam at a variable depth in the foil, and x the 
number of atoms per unit area in that depth. Leto), and o,,) denote the 
cross sections for loss and for capture of electrons respectively in the 
exchange He® = Het, and o,, and o,, the corresponding cross sections 
for the exchange He+ =~ He*++. Then, if we neglect the cross sections 
for capture and loss of two electrons in a single encounter, we may write 
for the change in composition of the beam in traversing a layer dz, 


dfy/dx=of,—oxho; 
df dao, fo—o fit orfe—Frsr, mera, CL) 
df /dx=02f,—Faf2- 

The mean free path for charge exchange in the He® = He* transition 
may be expressed by 4%),=(¢91-+0 9)! and for the Het = Het+* transition 
by 4%1.=(c,,+0,,)-1. Since the cross sections remain very nearly 
constant for a change of energy in the beam corresponding to one mean 
free path, we may expect charge equilibrium to be established for 7 > 42%, 
or 4x,,. With this condition satisfied the equilibrium composition of 
the beam is given from eqns. (1) by 


Fi(2) [fo(%)=201/6 10; 


(2) 
fo(@)[fi(@)=e 49/001. 

The thickness of the metallic foils used was several times the mean free 
path for charge exchange ; hence the measured intensities mj, v1 and n, 
respectively of He®, He+ and He** ions in the beam emerging from the 
foil were proportional to the quantities fy, f; and f, of eqns. (2). 


5.2. Huperimental Results 


Figures 3 and 4 show the values of the ratios n4/n and n,/n, respectively 
plotted (logarithmically) as functions of the energy H of the alpha-particles. 
The results show no difference in the charge state of the helium beam, 
within the accuracy of the experiment, for the different metals used. 
Figure 3 shows that the experimental points are closely fitted by a straight 
line, implying that oy,/c,) is proportional to H”’, where m’ is constant over 
the whole energy range studied. The value of m’ is 1°54--0-05. 

The plot of log (n4/n,) versus log H, fig. 4, appears to be linear for 
energies 0-43 to 1-1 Mev with a slope of 2-55-+-0-06 over this limited range. 
For energies below 0:4 Mev the results fit closely a relation of the form 

O19/9;=0:019—0-064H-++1-93H2 (H# in Mev). 

The equilibrium charge distribution of the helium beam as a function 

of the energy is shown in fig. 5. The values are taken from the smooth 
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curves through the experimental points for all the metals ner The 
fraction fy of the He® content in the beam obeys closely the relation 
f= —0-122-++0-119/E—0-0033/E?, : 
for energies between 0-2 and 0-6 Mev, and the curve above 0-6 Mev has 
been extrapolated using this relation. 
5.3. Hrrors 

The main error in the experiments arises from the uncertainty in the 
energy measurements, estimated to be within +5 kev. For energies 
above 0-5 Mev the error in each of the ratios ,/% and ny/n, is less than 
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The ratio o9,/0, of the loss to capture cross sections for the exchange 
He® = He* in metals as a function of the alpha-particle energy. 


2%, but increases to about 5% at energies below 0-2 Mev. In the results 
for beryllium a further error may arise for energies below about 0:4 Mev 
as a result of the subtraction procedure used to analyse the composition 
of the beam (see § 4). 

A few systematic errors are also possible (some of them have already 
been considered in §§2 and 3). Effects due to light and background 
radiation have been successfully eliminated, and by the use of a suitable 
collimating arrangement errors due to particles multiply scattered in the 
analyser wave guide have been avoided. Deposits of carbon or thin oil 
films on the target could have had serious effects. The former appeared 
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Fig. 4 
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Fig. 5 
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The equilibrium charge distribution of a helium beam in metals as a function 
of the energy. 
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to be serious only in the collodion-backed foils of aluminium, which were 
therefore not used. No oil diffusion pumps were used near the apparatus 
nor were there any visible traces of oil film on the targets bombarded. 
Checks were however made from time to time at a particular energy to 
test for consistency in the results. No foils were bombarded for more than 
two hours. The equality of response of the neutral detector to both 
charged and uncharged particles has been established, as discussed in § 3. 

In the study of the He® content, the energy spread of the scattered 
particles was between 25 and 30 kev. The spread, however, does not 
give rise to appreciable errors as the variations of the capture and loss 
cross-sections in this energy interval are small and approximately linear. 


§6. THEORY 
6.1. Electron Loss 


Bohr (1948), considering electron loss as a simple ionization process, 
showed that for velocities v of the particle much greater than vy=e?/h 
(the electron velocity in the first Bohr orbit of hydrogen), the cross 
section o, for electron loss should be proportional to (v)/v)" where n 
ranges from 0 to 2 depending on the charge number Z of the material. 
n is equal to 2 for the lightest nuclei, but should decrease as Z increases, 
and be zero for the heaviest nuclei where the atomic field is so strong that 
any collision in which the particle penetrates the atomic region would 
lead to the ionization of the particle. For velocities v approaching vp» 
there is no adequate theory of the interaction as neither classical (orbital) 
pictures nor simplified quantum ideas (Born approximation) can be 
applied. It is however expected that for v~vy, o,~7a,? where a, is the 
Bohr radius (=h?/me?). 

6.2. Hlectron Capture 


For the capture of electrons in lighter substances Brinkman and 
Kramers (1930), using the Born approximation, showed that for v >v9, 
o, (the cross section for electron capture) is proportional to v-!2. For 
heavier substances Bohr (1948) estimated that for v>v», o, should be 
proportional to v®. He used simple statistical considerations together 
with the assumption that capture was most probable for those electrons 
in the atom with orbital velocities comparable with the velocity v of the 
particle. As the velocity v approaches v», ¢, becomes large, but it cannot 
much exceed the geometrical cross section of the atom and therefore 
becomes less dependent on the velocity v. 


§7. Discussion 


It follows from the considerations of §6 that o,/o, should increase 
with increasing Z and, for vv), should be proportional to v™ where 
m~®5 for intermediate Z corresponding to the metals studied. As the 
velocity decreases m would be expected to decrease, and for U~, F,/F¢ 


should be of order unity. In the experiments o,/o, is given for He* ions 
by 7,/n and for He++ ions by n/n}. 
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In fig. 6 the experimental values of m are plotted as functions of the 
quantity (v/v)). For comparison some of the values obtained by Hall 
(1949) for protons in metals (except gold) and by Henderson (1925) for 
high energy alpha-particles in metals are included. The experimental 
results for the exchange Het = Het+ agree in the main with the theoretical 
expectations, though Henderson’s values of m for alpha-particles of higher 
energy (v >3-5v)) show a steady decrease with increasing velocity. On 
the other hand, the exchanges He® = Het (this work) and H® = H+ (Hall 


Fig. 6 


xe) 2:0 3:0 v/% 4:0 


The exponent m (in o/o,«v™) plotted as a function of v/v 9. The curves 
represent the following exchanges in metals: (a) He®=He* (present 
work), (6) Het=Het+ (present work), (c) Het=Het+ (Henderson), 
(qd) H°=H* (Hall). 


1949) exhibit no such agreement. It is significant that for both protons 
and alpha-particles the charge exchange behaviour is independent of the 
atomic number Z of the metal used. It may be that, in a metal foil, 
the charge transfer properties are determined by the lattice and electron 
band structure rather than by the electronic configurations of individual 
atoms. If this were so the usual theoretical treatments which deal only 
with single collision effects would be inapplicable. 
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The results obtained in the present experiments may also be compared 
with the results of Snitzer (1953) for the charge exchange of alpha- 
particles of energies between 0-1 and 0-48 Mev in the gases hydrogen, 
helium, air and argon. These do, by contrast, exhibit a definite depen- 
dence on Z. Snitzer’s results also show that in gases m increases with 
increasing velocity v for the exchange He? = He*, but is constant at a 
value between 5-1 and 6-3 in the exchange Het = He*++. The equilibrium 
charge distribution of a helium beam in hydrogen has an energy depen- 
dence remarkably close to that in metals, although the absolute values of 
fo: fx, and f, are not quite the same. 

The salient features of the charge distributions for both metals and 
gases are summarized in table 1. The velocity of the particle for which 


Table 1 


Stopping Eu Eos Ex» fi Ey 


material (Mev) (Mev) (Mev) | (max) Yo1/%o Y12/%o 


2-54 
2-58 
2-84 
2°56 
2-42 


Metals 0-157 0-370 0-645 0-66 


Hydrogen | 0-148 0-370 0-660 0-76 
Helium 0-145 0-430 0-800 0-72 
Air 0-098 | 0-340 0-650 0-76 
Argon 0-155 0-330 0-580 O20 


HORS 
S&S © bo bo bo 
CO = bo oO 


Eo,, Eo2, and Hy, are the energies corresponding to n»=n,, Nyp—=N, and ny=Ny 
respectively. 91, V;, are the velocities corresponding to H,, E45. 
Vo=—e?/h=2:19 108 cm sec-!. f,(max)=maximum value of Het 
fraction in beam, H, the energy at this value. 


o,=o, is significant since, according to any simple theory of electron 
transfer, the processes of capture and loss become equally probable when 
the linear velocity of the ion becomes about equal to the orbital velocity 
of the transferred electron. It is interesting to note that the results on 
the passage of ions through metals seem to bear out this hypothesis. 
Table 2 shows the extent of the agreement, the experimental values being 


. 


Table 2 
Orbital velocity of 
Exchange Sed aout transferred el 
(in e?/h units) : oS tae ca eres 
(in e?/h units) 
H? —<H+ 0-95 1-00 
He® Het 1-25 1-34 
He+t—Het+ 2°55 2°35 


taken from the present work and from that of Hall (1949), and the electron 
orbital velocities being defined by the known values of the corresponding 
lonization potentials. From the work of Snitzer (1953) it appears that a 
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similar correspondence can be made when the stopping material is a gas 
although the agreement is certainly not as good. 
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SUMMARY 


In the dilatation field of an edge dislocation in a metal the Fermi 
energy and ground state of the conduction electrons are both altered. 
To preserve a constant level of the Fermi surface everywhere, a charge 
shift occurs by a redistribution of the electrons and an electrical dipole 
is formed on the dislocation. Mott’s estimate of the effective charge 
on a solute atom in copper is used in a calculation of the electrical 
interaction of such an atom with an edge dislocation. For the copper- 
based alloys considered, the calculations suggest that this electrical 
interaction is between 1/3 and 1/6 of the corresponding elastic interaction. 


§ 1. INTRODUCTION 


From time to time evidence has been presented to show that the 
strengthening of metals by alloying elements in solid solution depends 
on the difference in valency of the solute and solvent atoms as well as 
on the difference in size. For example, Dorn, Pietrokowsky and Tietz _ 
(1950) showed that silver hardens aluminium twice as intensely as does 
zine, even though it produces no lattice parameter charge whereas zinc 
produces an expansion. Other alloying elements showed similar effects, 
although to express the effects of atomic size and valency in a single 
empirical formula it proved necessary to assume that the valency of 
aluminium was 2. In a similar investigation of copper alloys containing 
zinc, gallium, germanium, axid arsenic, it was shown by Allen, Schofield 
and Tate (1951), that all the alloys studied gave the same stress-strain 
curve at the same electron concentration, irrespective of atomic concen- 
tration. In interpreting their results the authors suggested that high 
positive charges on the solute atoms in solution might be the cause of 
hardening. This interpretation is not quite certain since Meyer (1951) 
has observed that the product ec is fairly constant for the alloys which 
gave the same stress-strain curve; here c¢ is atomic concentration and 
«=(1/a)(da/dc),_9, where a is the lattice constant. He suggests that 
lattice distortion could be the cause of hardening in these alloys. This 
does not of course rule out a valency effect since the extra electrons 
which remain localized to the higher valency atoms in solution may 
produce a local expansion of the lattice round these atoms. 


* Communicated by the Authors. 
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The purpose of the work described below was to examine the position 
from a theoretical standpoint ; to see what electrical interaction might 
be expected between a dislocation and a charged solute atom, and how 
this would compare with the corresponding elastic interaction. That 
an edge dislocation in a metal should possess an electric field is to be 
expected on general grounds, and the view has already been expressed 
by Mott (private communication) that there should be an electrical 
interaction with a solute atom. Landauer (1951) has already used the 
idea of an electrical dipole round an edge dislocation in another 
connection. 


§ 2. CHARGE POTENTIAL RouND AN EpGE DIsLocATION 
Following the usual Wigner-Seitz approximation we write the energy 
of an electron in a state with a wave number & in an unstrained lattice as 
h?k 
where Hy, is the energy of the ground state of the electron. In a 
homogeneously strained lattice the corresponding expression is 
| hk? 
— — 2 
E(k)=HE,+#,4+ oo oe a (2) 
where (H,)+,4) is now the energy of the ground state, 4 is the 
dilatation, and 
Ho = (CH OA) So: . . . . . . . (3) 
This distortion changes the wave number at the surface of the Fermi 
distribution from k, to k)(1-+-4)~-'%. The total energy of the n conduction 
electrons in the lattice is thus 


n | Bot 814+ 


A formula for #, is 


3 H2k_2(1+ A)23 
EON) i a 


meen reyes 1 Mmmenttie sae het.’ e- (8) 
where V(r) is the value at the boundary of the Wigner—Seitz cell of the 
periodic lattice potential (Bardeen 1938, Hunter and Nabarro 1953). 
Numerical values given by Mott and Jones (1936) for [V(7))—,] are 
2-3 ev for silver and 1-3 ev for sodium. For metals where the ion-ion 
interaction can be neglected, #, can be found quite easily from the Fermi 
energy ; in such cases the equilibrium lattice parameter is determined 
by the condition that the derivative of the total energy of the conduction 
electrons with respect to 4 should vanish at J=0. Applied to (4) this 
gives tae 
= é five ‘ (6) 

5 2m 
Since the effect of the ion-ion interaction is only to increase the equili- 
brium spacing a certain amount (about 10% for copper), from the value 
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at which the energy in (4) is a minimum to a value where the slope of the 
E,(r) curve is not greatly altered, the formula in (6) remains a fair 
approximation to #, in cases where the ion-ion interaction is not 
negligible ; we shall therefore use it in the following estimates. 

To deal with the inhomogeneously strained lattice in the field of a 
dislocation we first regard the crystal as divided up into a number of 
small pieces, each homogeneously strained to the average value in that 
part of the field. So far, each piece is electrically neutral. However, 
the level at the top of the Fermi distribution now differs from one piece 
to another. To correct this, and make the Fermi surface level everywhere, 
the electrons redistribute themselves so that there is a shift of charge 
from piece to piece (Hunter and Nabarro 1953). 

Before this charge shift, the energy at the Fermi surface, in a piece 
where the dilatation is J, is 


irk 
Ey +H, A+ —™ (144), 
which to the first order in J is 
aN ae Cad (ae 
Bot SY + (B- 35,5) Ae on Re. ee 


The third term here measures the difference between the local level of 
the Fermi surface and its average value. This difference is removed by 
the charge shift. Hence the charge potential at this place, due to the 
rearrangement of the electrons is, substituting for Z, from (6), — 


2], 2 
4 h'kg 
15 2m 


A; . & 


apart from a constant which measures the average dilatation through 
the lattice and vanishes for a dislocation in linear elasticity theory. 
This charge potential is the same as the electrical interaction energy of a 
unit charged ion placed in this part of the field. 


§3. INTERACTION WITH A SOLUTE ATOM 


For an edge dislocation the shift of electrons is towards the expanded 
side of the dislocation field, and away from the half-plane. The dilatation 
is greatest in the region within one atomic spacing of the edge of the 
half-plane, and its value here is about 1/8. Hence the electrical 
contribution to the binding energy of a solute atom and a dislocation 
is of order 

1 h*k,? 
30 Din otek ae ae 


per electronic unit of charge on the atom. For a unit charged ion in 
copper this is about + ev. 

According to the calculations of Mott (1936, 1937 ), the effective charge 
on a 2-valent ion in copper is about 0-075 (2—1) unit electronic charges, 
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We shall use this result, although it may be seriously in error when the 
ion is in a heavily strained region of the lattice, such as exists at the 
centre of a dislocation. The electrical binding energy of the ion with an 
edge dislocation should thus be about 0:02 (z—1)ev; that is, about 
0-02, 0:04, 0-06, and 0-08 ev respectively, for zinc, gallium, germanium, 
and arsenic, in copper. The corresponding elastic binding energy, from 
the formula of Cottrell (1948), is as follows : 


Zn Ga Ge As 
€ 0-062 0-078 0-090 0-150 
Elastic binding } 


: 0-125 0-163 0-188 0-272 
energy, in ev 


Here « has been calculated from the lattice parameter measurements of 
Hume-Rothery, Lewin and Reynolds (1936), and Owen (1947). 

The elastic binding energy is thus some 3 to 6 times greater than the 
electrical one in these copper-based alloys. It appears that in these 
alloys the role of the electrical interaction is secondary to that of the 
elastic interaction. Its importance in other alloys cannot be ascertained 
until more is known of the effective charges on solute atoms in them. 
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ABSTRACT 


A convenient method of representation is proposed (§ 2) for data on 
z-meson decay configurations, applicable when the unlike outgoing 
m-meson is not distinguished. The relation between the spin and parity 
of the 7-meson and the distribution of decay configurations is obtained 
for some simple cases. The hypothesis that the 7- and y-mesons are 
identical requires a non-zero spin for this particle and the available 
data on t-meson decay does not exclude this possibility. However, 
observations in which the unlike outgoing 7-meson is not distinguished 
are relatively ineffective in discriminating between the various possibilities. 
The distortions which strong meson—meson attraction may produce in 
r-decay configurations are discussed ; the present data offers no evidence 
on this effect. 


§1. IyrropucTION 


THE existence of a t-meson whose mass is 978+-6 m, and which decays 
into three z-mesons is now well established (Brown ef al. 1949). Although 
mass measurements on the outgoing particles in high energy cosmic 
ray stars suggest (Daniel and Perkins 1953) that only a small proportion 
(<20%,) of the heavy mesons produced could possibly be 7-mesons, the 
number of r-decay events observed in photographic emulsions has 
become quite considerable, presumably because of the very characteristic 
7-decay pattern. About fourteen clear examples of the process 
T*->n*-+nt+-+n~ are now available. In each event the angles between 
the outgoing meson tracks, and also the individual meson energies, have 
been measured, but it is not known which track is produced by the 
meson whose charge is opposite to that of the 7-meson (we shall generally 
refer to this 7-meson as the ‘ unlike’ meson). In § 4, it will be seen that 
the frequency distribution predicted for the decay configuration depends 
rather less sensitively on the nature of the. 7-meson when this additional 
information is lacking ; however, since such a complete identification 
of the decay process may well not be generally possible until r-mesons 
may be artificially produced, it seems worthwhile to present an analysis 
appropriate to the data available at the present stage and to discuss 
what one may hope to learn concerning the 7- and 7-mesons when more 


ee ee 
* Communicated by Professor R. E. Peierls, 
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events have been observed. The r-meson properties of first interest are 
those related to the conservation laws of physics, the mass, charge, 
Intrinsic spin, parity and (possibly) isotopic spin. The charge of the 
7-meson is known to have magnitude |¢| and the experimental data on 
its mass has recently been discussed by Rochester and Butler (1953). 
Concerning the isotopic spin properties, one type of observation which 
could lead to some information has been discussed recently (Dalitz 1953). 
A phenomenological description of the decay for assumed spin and 
parity may also allow the determination of other parameters relevant 
to the 7-meson, which may provide guidance in a more detailed theory 
of the 7-meson. The existence of any strong interactions between the 
outgoing 7-mesons would be expected to influence the decay configuration 
considerably (Brueckner and Watson 1952), so that a study of the decay 


events may also provide some evidence concerning the meson-meson 
interaction. 


§2. THE SPECIFICATION AND ANALYSIS or 7-MEson Decay EVENTS 

In the decay of a 7-meson at rest, the three outgoing 7-mesons have 
zero total momentum and a total kinetic energy of H=(m,—3m,)c?, 
so that, apart from the spatial orientation, the specification of the decay 
configuration requires two parameters. A convenient choice would be 
the energy of the unlike meson and the magnitude of the difference 
between the energies of the like mesons. However, it is not possible to 
distinguish between the outgoing mesons in the present data, so that a 
choice of parameters which is symmetrical between the three 7-mesons 
is more appropriate. Denoting the total kinetic energy of the 7-mesons by 
E and the individual kinetic energies by ¢€,, €9, €3, define 


a =e,—H/3 (esl; 2, 3). 5 © 5 5 . 5 (2.1) 
These «, then have zero sum, thus X,a,=0. If H and the «, are known, 


the decay configuration is specified. If the unlike meson is unidentified, 
two symmetrical functions, 


Y= — (a4 5 +o 9%g 109% ,)—=— LHe, DHAyeotg, . . (2.2) 
may conveniently be used to specify the configuration uniquely, since, 
for given (Y, Z) the a, %», a are the three roots of the cubic equation 

ere Vico Lia O0) 8 sn GQ yah ad « (208) 
Y and Z are limited firstly by the condition that (2.3) has three real 
solutions, which requires that Z?<4Y%/27, and secondly by the require- 
ment that the outgoing mesons have total momentum zero. A necessary 
and sufficient condition that the momenta p,, Po, p; satisfy the triangular 
inequalities p,+p.>P3, Pot Ps2>Pv Pst PiZzP2 which are necessary 
and sufficient to allow a choice of p,, Ps, ps satisfying py+P2+P3,—29, 
is that 


FE? 
(py + P2t+Ps)(P1-+P2—Ps)(P1—P 2+ Ps)(P1 + P2— Ps) = 16m," (5 —¥ ) Za 
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so that Y, Z are limited by the inequalities 

Y<H412 and Z<v/(4Y3/27)< #9108. . . . . (2.4) 
As the uncertainty in the measurement of individual meson energies is 
often quite large, it will be more appropriate to use the ratios of the 
meson energies, which may be deduced from the more accurately known 
angles between the meson tracks. Accordingly, co-ordinates (A, 6) which 
are functions of these ratios, may be defined by 

Y =)H#?/12, Z=3/2 sin H3/108=yH3/108,, . . . (2.5) 
where it is convenient to define the quantity 7=A*”? sin@. The inequalities 
(2.4) then correspond to 0<A<1 and —2/2<0<7/2. 

The kinetic energy of an outgoing 7-meson cannot exceed 24/8, 
which is close to 50 Mev since the experimental data gives H=76-5+3 Mev. 
Consequently it is reasonable to adopt a non-relativistic description, 
since the resulting errors are at most a few per cent, small compared with 
other uncertainties in the data. The expression for the energy 
distribution of the three 7-mesons may be written in the general form 


F (Eq, €g, €3)8(€y + €2+€3—£)3(P1+ Pot Ps)d*p, Eps d*ps, 
where F(e1, €5, €3)=F'(€g, €1; €3) since the like mesons (1, 2) are indistinguish- 
able. Non-relativistically this may be reduced to 
827m," F'(€1, €9, €3)0(€y +€g+€3—LE) de, deg deg. 
If the unlike particle cannot be identified the same configuration may 
occur in three ways and the total frequency is then the sum of the three 
separate frequencies, thus 
[87?m,?(F(€1, €2, €3) + F(€g, €g; €1) 

+ (és, €1, €2))] 8(€; +e. +e,—H) de, degdeg. . . «. (2.6) 
The function within the square brackets is a symmetrical function of 
€1, €9, €g and thus of «,, %, «3, so that it may be expressed as a function 
of # and of Y, Z and hence of A, 0, say 4(A, =A, n). If we write 
C=a,+%.+a5, then (2.6) may be written as 
O(%1, X, Hs) 

0(f, A, 7) 
The ratio of the volume elements da, da,dz,; and d¢ ddd», namely 
O(%1, X», %3)/0(C, A, n), is equal to the magnitude of 
(%1—%)(%2—%)(%3— 1) 1296/29, 

which equals 4/(4Y%—27Z?) 1296/H®=364/3)3/2 cos 6/H?. Integrating 
over ¢, (2.7) becomes, apart from some numerical factors, 


b(X, 0) dd dO. ee MER We ah ih 2): 


The following considerations lead to a very convenient realization of 
the co-ordinates defined above. For an equilateral triangle YUV (see 
fig. 1), the sum PL+PM-+-PN of the perpendiculars to the sides from 
any point P within the triangle equals the altitude of the triangle. Hence 


B(A, 0)3(L) dary day dag=(A, 8) 8(C) dL OA. dna acct, Coke 
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a decay event may be uniquely specified by a point P such that the 
perpendiculars (PL, PM, PN ) are proportional to the meson energies 
(€,, Eg, €3). Interchange of two of these energies corresponds to the 
reflection of P in a corresponding altitude of the triangle so that, 
corresponding to the incomplete identification of the mesons, the event 
is represented a point in each of the six sub-triangles of YUV. The 
distribution in any sub-triangle is obtainable from that in AOV by 
successive reflections in the altitudes, so that one need consider only 
that ordering of the meson energies for which P lies within AOV. The 
Cartesian co-ordinates (a, y) with respect to the axes OX, OY are given by 


t= /3(ey—€5)/H, y=(2eg—€—€g)/E, 
where the altitude of the triangle has been taken as 3 units. It may 


now be readily verified that the polar co-ordinates (r, ¢) of P are 
(A'?, —@/3). The points representing decay events are therefore confined 


Vv 


to the interior of the circle inscribed to the triangle YUV, since momentum 
conservation requires A<1. In plotting decay events on this diagram, 
the momenta have been calculated from the observed angles between 
the meson tracks, and the energy ratios from the squares of these 
momenta, so that the representation expresses angular relations in the 
decay and is not modified by relativistic effects. When the co-ordinates 
are given, the corresponding configuration of the decay process may be 
quickly obtained by reference to this diagram. 

If the function (A, 6) is independent of 4, @ which is the case, for 
example, when F is independent of ¢€,, <2, <3, the experimental points 
should be distributed randomly over the area AOB, since the element 
of area is rdrdé—da d6/6. The limited data available is displayed in 
fig. 2 and although there appears to be some tendency to favour the 
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upper and outer part of AOB the distribution obtained is not yet 
significantly different from a random distribution. The corner A 
corresponds to decay events in which one 7-meson is rather slow and a 
tendency to favour this corner may perhaps result from some _ initial 
experimental bias favouring detection of this particular decay configura- 
tion. The distribution to be expected according to various hypotheses 
will be obtained in § 4, where it will be shown that the present statistics 
are insufficient to allow discrimination between the alternatives 
considered. 
Fig. 2 


8 
The data from 13 7-meson decay events. 


$3. CONSIDERATIONS ON THE NATURE OF THE 7t-MESON 


The questions of most direct interest concerning the t-meson relate 
to its intrinsic spin, its parity and to the applicability of the notion of 
isotopic spin. If the 7-meson is some kind of meson—nucleon complex 
occurring in consequence of the interaction between the meson and 
nucleon fields, it would follow that the 7-meson would have a definite 
isotopic spin, if charge independence be established for the meson—nucleon 
interaction. If the r-meson needs specification by some new field coupled 
directly to the meson and nucleon fields, these interactions may still 
satisfy charge independence and the 7-meson may still have a definite 
isotopic spin. Some evidence concerning the applicability of isotopic 
spin to the 7-meson may be obtainable from the study of the branching 
ratio for the decay processes t*>7*-+-a++-n7- and +r+37++79+79, as 
discussed earlier (Dalitz 1953). At the present, there is perenne no 
‘evidence requiring that the 7-meson should satisfy charge independence : 
nevertheless, in this and the following section, we shall indicate what 
additional consequences would follow if this notion were applicable. 

The configuration of the three mesons resulting from z7-meson decay 
will depend on the spin j and parity w of the z-meson. The configuration 
is completely specified by the momenta p,, p, of the like particles and 
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the matrix element for decay to this configuration must be an irreducible 
tensor of rank j and parity (—1)%w, formed from p,, p, and symmetrical for 
interchange of the suffices 1, 2, since the 7-meson satisfies Bose statistics. 
If w=—(—1) such a tensor is easily found (e.g. Y,,,’(p,)-+ Y,,/(Ps)) 
so that the decay is then not forbidden by angular momentum and parity 
conservation. If w=-+(—1)’, then for every j except j=0, such a 
tensor may be formed, for example, by multiplying the tensor just formed 
for CS 1, w'=—(—1)’"1) by the pseudovector p, p,(p,2—p.2) and 
extracting the irreducible tensor of rank j from the product. For j=0, 
however, we are required to form a pseudoscalar from the vectors pj, Po, 
which is not possible. The 3 7-decay, therefore, is forbidden by angular 
momentum and parity conservation only for a scalar heavy meson. 

The selection rules for decay of a heavy meson into (a) two 7-mesons 
or (6) one z-meson and a photon, which are allowed as far as energy and 
momentum conservation are concerned, are well known. For case (qa), 
they may be simply obtained by noting that, since the final state is 
specified by the relative momentum p of the mesons, the matrix element 
for decay of a particle of spin j must be given by Y,,,/(p). Since this 
tensor has parity (—1)’, the 27-decay is therefore forbidden by angular 
momentum and parity conservation for w=—(—1)’. If isotopic spin is 
a good quantum number, then since interchange of the mesons 
corresponds to p-+—p and the mesons satisfy Bose statistics, the 2 7-decay 
is allowed only for even or odd j according as the isotopic spin is even or 
odd. For states with even j and odd isotopic spin, the 2 7-decay may 
still proceed unfavourably through electromagnetic interactions which 
violate charge independence, though one might then expect that the 
decay +++7*++7°+y would be more probable. For case (b), the final 
state is specified by the relative momentum p and the photon polarization 
vector e which satisfy e . p=0 since the photon is a transverse wave, and 
the matrix element for the decay to this state must be a linear function 
of e. For every j (except j=0) such a matrix element may be formed 
for parity w=(—1)’ or —(—1)’ by multiplying Y,,’-\(p) by exp or e 
respectively and extracting the irreducible tensor of rank 7. However, 
no non-zero scalar or pseudoscalar can be formed, so that angular 
momentum conservation forbids (7, y) decay for a heavy meson of spin 
j=0. In view of the number of 3 7-decay events observed, it has been 
natural to suppose that the 7-meson is pseudoscalar, since the 3 7-decay 
would then be the only allowed decay process involving 7-mesons alone 
and other decay processes involving photons (e.g. r+>7=-++7°+-y) might 
well be less probable. 

Recently, however, a number of events described as y-meson decay 
have been found from which a charged 7-meson of unique energy 11.0 Mev 
results (Menon and O’Ceallaigh 1953). If it is supposed that the neutral 
particle produced in this decay is a 7°-meson, these events may be 
interpreted as the decay of a particle of mass about 980 m,. Also the 
analysis of V°-particle decay (Armenteros ef al. 1951, Barker 1953, 
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Thompson ef al. 1953) strongly suggests that there exists a V.°-meson 
which decays into two 7-mesons with a Q-value corresponding to a V,° 
mass of 962 m,. It is very plausible that this V.° meson is the neutral 
counterpart of the y-meson and an attractively simple representation of 
the data is offered by the further hypothesis (Menon and O’Ceallaigh 1953) 
that the 7- and y-mesons are identical, the events observed representing 
two alternative decay schemes of the (7, x) particle. In this case the 
(r, x) particle could not have spin j=0 but must be a particle of parity 
(—1) if its spin is j(>1)—if isotopic spin is rigorously applicable, j is 
restricted to odd values. On the other hand, the hypothesis that the 
neutral particle in y-decay is a photon is still compatible with the 
hypothesis of a (r, x) particle and also requires that this particle have 
440. Ifthe parity of this particle is (—1)’, then both 27- and (7, x)-decay 
may compete with 37-decay, although, for j even, the 2 7-decay will be 
relatively unfavoured if charge independence is a good approximation. 
If the parity is —(—1)’, the y-decay would have to be (z, y)-decay. 
The hypothesis that the 7-meson is a particle of spin 1 is therefore of 
particular interest, and if the hypothesis y=7 is to be maintained, we 
must examine whether the consequences for the 37-decay pattern are 
then in agreement with the data. 


§4. Discussion oF SoME ParTICULAR HyPOTHESES 


The decay of the z-meson into the z-meson states of momenta 
Pi; Po, Ps (in this and the next section, the suffix 3 will refer to the unlike 
meson) will be described by a matrix element M(p,, py, py). Corresponding 
to the analysis of the final state into partial waves of definite angular 
momentum, this matrix element may be written as the sum of a number 
of terms of definite form with coefficients depending on the energies of 
the mesons. It will be assumed that only low angular momenta are 
important and that these coefficients may be expanded in powers of the 
meson energies. ‘This supposes the region of production to be small 
compared with the meson wavelengths and neglects the effects of any 
strong interactions between the mesons after the decay. These effects 
will be discussed briefly in §5, following the viewpoint of Brueckner and 
Watson (1951) that the main features of these interaction effects will not 
depend on the details of the production process. 


(1) The Hypothesis of a Pseudoscalar r-Meson 
Since there is no change in intrinsic parity, M (Py, Po, Ps) is a scalar 
whose general form is 
M(P,, Po Ps)=Mo(p1, Po, Ps) +P, - PoMs'(P1, Po, Ps) 
+A(P1 + PsP2 + P3—$P1 + PoPs”).M5"(p, Po Ps) +++. » (41) 
of which the successive terms describe the emission of s-wave mesons, 


of one s-wave and two p-wave mesons, and of one d-wave and two p-wave 
mesons respectively. The summations are over all permutations of (123) 
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and M(p,, p2; ps) is required to be symmetric for interchange of (12), 
because of the Bose statistics. By expanding the scalar coefficients in 
powers of p,?, p,”, p,2, making use of the energy conservation relation 
(1? +P2?+p3")=2m,H, momentum balance p,-+ Potp3;=0 and the 
symmetry of M, (4.1) may be reduced to the form 


M=A-+dmg+ba,%5-+co42-4 aS Loe 200) 


A, a, b, c being constants and «,, «5, «3 being defined by (2.1): In (4.2), 
A contains contributions from all partial waves, a arises from the energy 
dependence of the s-waves and from s—p interference, b and c from 
the p-waves alone or from higher effects, and so on. The function of 
interest is 4(A, 0) which is proportional to 


2(AA*+ 0,22 (aA*)+a,2(aa*+22(cA*)) + 20,0,2%(bA*)) 
=3A A*+ A(2aa*+ 2A(2cA*—bA*))H2/12+. . . . . (4.8) 


Although the term a in (4.2) allows a considerable range of variation in 
the energy distribution of the unlike meson, the cross terms between 
a and A vanish in the symmetrized function (A, 0) whose (A, 0) dependence 
arises from higher order terms. The fact that the unlike meson cannot 
be identified leads here to an insensitive dependence of the observable 
phenomena on the details of the theory. Consequently, apart from 
meson—meson interaction effects, the decay may be expected to be 
isotropic. Also, if a is negligible, the orbital state of the outgoing mesons 
is totally symmetric (of the representation [3]) which implies a branching 
ratio of + for the 7-decay involving neutral mesons if isotopic spin is a 
good quantum number (Dalitz 1953). The term a, and possibly the 
effects of meson—meson interactions, will lead to an admixture of [21] 
orbital symmetry which would increase this branching ratio. 


(2) The Hypothesis of a Vector t-Meson 


The matrix element must be a pseudovector and, including the 
possibilities up to a final (p?d) or (sd?) configuration, its form is uniquely 


M(P4, Ps, Ps)=BpiX Po(Pi?—P2”)) - - - + (4-4) 
Averaging over the orientations of the initial 7~-meson then leads to 
P(e, €9; <3) =4m,'B?(€y—€ 9)" [€3(2H — 3e3) — (ey —€ "1/3, - (4.5) 
which corresponds to, a distribution function 
PURO Ahem e cata ets Ge. 3a» , (4.6) 


Since the matrix element (4.4) is zero if two of the 7-mesons have parallel 
paths or if the like mesons have equal energies, it may be easily under- 
stood why the centre and edge of the circle in fig. 1 are unfavoured. 
The matrix element (4.4) leads to a final state of orbital symmetry [21] 
which corresponds to equal probability for the two modes of 37-decay 
for the 7-meson if the notion of isotopic spin is applicable. In this case, 
since the only decay observed with certainty to yield a single charged 
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m-meson is the y-decay, the hypothesis of (7, x) identity would be 
compatible only if this decay is the predominant mode of decay for the 


particle. 
(3) The Hypothesis of a Pseudovector r-Meson 


The matrix element must be a vector. The simplest decay is emission 
of two s-wave and one p-wave particle for which M(p,, Pz; Ps) is 
proportional to ps, giving 

Fley, €4€) =Cegr ae ee 


so that the distribution function is the constant 
b(A; 0) == 80771 CE ee ee ee (4.8) 


This matrix element leads to a final state of orbital symmetry [21], its 
symmetric part 4(p,+p.+p;) being zero owing to momentum conserva- 
tion. Remarks similar to those at the end of the last paragraph apply 
also to this case except that here the y-decay must be of the (7, y) type. 


(4) Higher Spin Values 


For spin j, even the simplest forms for M rapidly increase in number 
with j, so that the predicted distributions offer a considerable range of 
possibilities. For the cases 2+, 3—, however, the simplest / is unique 
and leads to the expressions (l—A) and (1—A)(1—+4A) for the function ¢, 
respectively. 

For cases (1), (2) and (3) the energy distribution of the unlike meson 
is plotted in fig. 3. To emphasize the loss of sensitivity resulting from 
the incomplete identification of the decay, the energy distributions 
predicted* for the slowest and fastest meson and for all mesons, 
irrespective of charge, are plotted in fig. 4 and compared with the 
available data. This effect is particularly marked in the pseudovector 
case, where the distributions of fig. 4 agree with those for the pseudo- 
scalar case, despite the marked difference in the distribution for the 
unlike meson. The data available at present are clearly inadequate to 
discriminate between the various possibilities discussed here. 


* Given the function (A, 7), the energy distribution of all particles is given by 


v(1—y*) 
Ple= [PO (a2 +42, y(y®—B2)) de, 
where y=(3e—L)/E. The distribution Q,(¢) of the slowest particle is given by 


Qs(€)=P(e) for «<E/6, since ¢ is then necessaril Ei 
for H/6<e<E)3, 1ecessarily the slowest particle, whereas 


—yv3 9 9 9 € 
Qsle= [> "° Oa +y2, yly?—B22) dar, 
corresponding to integration over the sector OCB of fig. 1. T istributi 
tion over g. 1. The distribution 
Qs(<) of the fastest particle is given by Q;(e)=P(e) for ¢>H/2, whereas for 


E/3<e<E/2, Qs(€) is found by integrating ® from 0 to --y4/! 
4 x ? : . y to ey F 
differ in form only in virtue of any a Sas of ®, ae Oe ee 
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§5. Some Posstpt—E Mrson—Meson [ntTEeRaction EFrrects 
In the emission of 7-mesons in cosmic ray stars, a marked correlation 
in the directions and energies of pairs of 7-mesons has been found by 
Danysz, Lock and Yekutieli (1952), who tentatively suggested that these 
meson pairs may result from the decay of a short lived £°-meson (see 


Fig. 3 


o fa € 26g 


The energy distribution of the unlike z-meson in r-meson decay for several 
T-meson spin values. 


Fig. 4 
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Energy distribution for all 7-mesons from 7-meson decay : comparison with 
histogram for 13 events. 
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Energy distribution of the slowest and fastest 7-meson in r-meson decay. 
Experimental histogram based on 13 events. 
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also Danysz et al. 1953). However, Brueckner and Watson (1952) have 
remarked that the existence of a strong attractive interaction between 
z-mesons could provide an alternative explanation for these correlations 
and have pointed out that, in any case, the study of correlation effects 
near the threshold for multiple meson production would provide 
quantitative evidence on the nature of the meson—meson interaction. 
A similar point of view has already lead to some information on the 
neutron-neutron interaction (Watson and Stuart 1951) from the 
distortion of the y-ray spectrum emitted in the capture of 7~-mesons by 
deuterium due to the neutron—neutron correlations resulting from their 
interaction. 

In r-meson decay, three 7-mesons are produced with low kinetic energy 
and the existence of strong attractive meson—meson interactions might 
be expected to distort the distribution of decay patterns in a characteristic 
way. The kinetic energy of the relative motion of two mesons (1 and 2, 
say) is given by 

A=(p,—P2)?/4M =H —3e,/2, 

so that, since the interaction is most effective between mesons of low 
relative momentum, this effect would favour a configuration in which 
the third meson has energy near to the maximum 2#/3. The corresponding 
angle between the paths of the interacting mesons is about 4/(34)/5, 
or less, i.e. about 30° for 1~2 Mev (a value for 4 which is typical of the 
cosmic ray data at very much higher total energies). The parameters 
(A, 0) of these configurations have values close to A=1, 6=7/2, so that 
the distribution of the decay events would then crowd towards the 
corner B of fig. 1. 

Following the discussion of Watson (1951), the matrix element leading 
to a state in which two mesons have small relative momentum k has a 
k-dependence essentially given by a factor (sin 8)/k, when it is supposed 
that meson—meson scattering is strong for small & in only one partial 
wave (assumed here, for simplicity, to be the s-wave) of phase-shift 8. 
If, to illustrate the effects quantitatively, it is supposed that the basic 
process is otherwise isotropic and that only unlike mesons attract 
strongly, the matrix element is essentially C(sin 843/K13)(Sin 5o3/ko3), and 


P(E4, €g 5 €3)=C*/(a+ (py—Ps)*)(@2+ (Po—ps)?), . . (5.1) 
if it is supposed further that k cot 8=«x—=h/a, a being the zero-energy 


scattering length. With »=«?/2M#, the corresponding symmetrized 
function is then 


H(A, 0) =C?(1 +2) M2B*(u(3-+2p.)?-+3(1+4)(1—A) + (1—A?? sin 6)) 
=C?/M?H?(9u+3(1—A)+1—A%2 sin f), . . . (5.2) 
since the values of « of interest are small. If the interaction effect is 
strong for relative energies of up to ~J, then pu~2A/E, which is 0-05 


for J=2mev. The expression (5.2) will not be a good approximation 
for points not close to B since both k,, and ky will then be large and the 
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approximations leading to (5.1) are not valid. Even for a point close to B, 
one of the three corresponding configurations will have both ky, and Ko 
large, but the error thus made will be small compared with the large 
mecrease in the probability of the other two configurations. Hence if 
the z-meson is pseudoscalar and there is an attractive interaction 
between mesons comparable with that necessary to explain the 
correlations in cosmic ray stars, a considerable fraction of the 7-decay 
events should lie quite close to B. The present data (fig. 2) shows no 
evidence for any very strong effect of this kind. 

The vector meson deserves separate consideration since its matrix 
element favours configurations in which the relative momentum of each 
meson pair is high. For this case, 


[e3(2H — 3e3) —(€y—€g)” (eq —€ 9)” 


o4 a” 
(7 +28 —3e,] (5 +2H—3e,] . (5.3) 


P(A, 8) =4(1—A)(A+-$m)/ (92+ 3(1—A)+(1—n)), 

assuming yu is small. This expression remains bounded as p>0; though 
the meson—meson interaction increases the probability of an event near 
to B, the effect does not crowd a high proportion of events into this 
corner. Away from this corner, (5.3) shows a considerable 6-dependence, 
which depends, however, on the incorrect assumptions underlying (5.1) 
when k,,, ky, are large ; a genuine 6-dependence will, of course, occur near 
the outer edge of AOB. The present limited data of 7-meson decay does 
not conflict with the hypothesis of a strong meson—meson attraction in 
the case of a vector 7-meson (or of a (2+) or (3—) 7-meson). It should 
finally be noted that in all cases, the upper half of the region AOB should 
suffer very little distortion as a result of this interaction effect. 


Prey, €g 5 €3)= 


§6. CONCLUSIONS 


The data available at present on the decay of the 7-meson offers no 
significant evidence for any lack of isotropy in the decay process nor for 
any distortion of the decay patterns as would be characteristic of a 
strong attraction between 7-mesons. On the other hand, the data does 
not at present exclude the possibility that the 7-meson is a vector particle, 
in which case the absence of evidence for any correlations characteristic 
of strong 7-7 attraction would not be significant, their effect not being 
particularly marked in the vector theory. However, theories giving 
quite different matrix elements have been found to lead finally to very 
similar results, so that a comparison with the evidence would not be 
sensitive to the r-meson spin and parity. This insensitivity results from 
the need for averaging over three mesons whose energies are limited by 
energy and momentum conservation—a knowledge of the energy 
distributions for the unlike meson alone would provide a more sensitive 
indication of the meson spin. In this respect, the r-meson decay contrasts 
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sharply with the decay event ,.>e+2v in which less quantitative data is 
available from each event but in which the electron spectrum has a fair 
dependence on quite detailed -meson properties. If it should prove 
possible to analyse 7-meson decay events in more detail (e.g. if the 
z-mesons which decay were all of positive charge, as is conceivable, 
identification of the 7~-meson in each event would suffice) a more direct 
comparison with the functions F(€,, €2 ; 3) will become possible and a 
more sensitive analysis of the data be feasible. At that stage a more 
detailed study of these functions will be justified. 


In conclusion, I am pleased to thank Dr. E. P. George, Dr. P. Hodgson, 
Dr. M. G. K. Menon and particularly Dr. W. O. Lock for information 
concerning the experimental data on t-mesons and to thank Professor 
R. E. Peierls for his encouragement and advice during this study. 
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ABSTRACT 


Evidence for a y-ray line at about 12-5 Mev, in the 7Li(py) spectrum has 
been put forward recently by several workers. Although the measure- 
ments differ amongst themselves it appears that a line is present to about 
10% of the intensity of the well known 17-6 Mev component. Experiments 
in this laboratory have shown that, in addition to the broad 14-8 Mev line. 
there are other transitions from the 17-6 Mev state of Be to lower even 
states. The results were obtained by observing the emitted y-ray 
coincidence with one of the «-particles resulting from the breakup of the 
even state. The measured «-particle ranges show that transitions occur 
to §Be levels at 4-09, 5-31--0-06 and 7-51--0-06 mev. The corresponding 
y-rays have energies of 13-54 Mev, 12-32-L0-06 and 10-12-L0-06 Mev and 
intensities of 2%, 2% and 0-5% of the total intensity. Backgrounds of 
soft x-radiation from the target backing, and energetic «-particles from 
the reaction “*Li(pa)*He complicated the experiment necessitating the 
performance of check experiments to ensure correct interpretation of 
the results. In particular the time of flight of the «-particles was used 
to give a rough indication of their energy and demonstrate the reliability 
of the coincidence counting. 


§ 1. PREvioUS OBSERVATIONS BY OTHER WORKERS 


THE gamma-ray spectrum from the reaction ‘Li(py)*Be was studied by 
Walker and McDaniel (1948) using an electron pair gamma-ray spectro- 
meter. They found two lines, a narrow one at 17-6 Mev and a broad one, 
2mev wide, at 14:8 Mev, with an intensity ratio of 2 to 1 when the 
bombarding proton energy was 460 kev. No other component was 
detected. Three particle stars from the reaction “C(y, 3x) produced by 
the radiation from ’Li(py) have been observed by Nabholz, Stoll and 
Waffler (1951) and by Goward and Wilkins (1952). The former reported 
evidence for a line at 12-6 Mev with an intensity of about one tenth of the 
17-6 Mev line, while Goward and Wilkins reported the energy to be 12-3 Mev 
and the intensity one half that of the 17-6 Mev component. Recently, 
Titterton (1953) used the radiation from this reaction to induce the 
disintegration 7Li(yt)*He and the results of the experiment gave evidence 
for a line of 12-5 Mev energy with intensity about one twentieth of the 


17-6 Mev line. 


* Communicated by Professor EK. W. Titterton. 
SER. 7, VOL. 44, NO. 357.—OCT. 1953 4B 
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McDaniel has stated in a private communication that a line at 12-5 Mev 
having a width less than 1 Mev would have been observed if the intensity 
were one tenth of the one at 17-6 Mev. Ifthe width were 2 Mev about twice 
that intensity would have been required for detection. 

Burcham and Freeman (1950) have shown that the 14-8 Mev component 
is produced by a transition to the 3 Mev level in*Be which then disintegrates 
into two alpha-particles. The following gamma-alpha coincidence experi- 
ments were carried out to determine whether weaker gamma-ray lines are 
present, and if so, to which levels of the Be nucleus the transitions occur. 


§ 2. APPARATUS 


The target chamber and counters were arranged as shown in fig. 1. 


Fig. 1 
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The arrangement of the target chamber, absorber supports and counters. 


The alpha-partieles were detected by a scintillation counter in which 
the scintillator was a layer of stilbene, 6 mg/cm? thick evaporated on to a 
glass slide. This was attached to the end window of the photo-multiplier 
tube which itself formed the vacuum wall across the tube B. Alamein 
foil of 0-66 cm air equivalent prevented light from the target and low 
energy particles from entering the crystal. The solid angle of acceptance 


of the counter was 47/138 sterad, the cryst : 
> al being 4-4 
of the target. J £ cm from the centre 
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The scintillator in the y-ray counter was a Nal crystal 1 cm thick by 
2-cm square placed with its centre 4 cm from the target in the plane of the 
beam and the «-counter. Aluminium absorber foils of graded thickness 
were mounted on frames D, supported from the two plungers A and C 
which could move into the end of the target chamber through Wilson 
seals. These foils could be interposed between the target and the alpha 
counter, so that the air equivalent of the total thickness of absorber 
between the target and counter could be varied from 0-66 cm to 7 cm in 
steps of about 0-1 cm. 

Big.*2 
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Block diagram of recording units. 


A block diagram of the recording equipment is shown in fig. 2. The 
fast coincidence circuit was similar to that described by Bay (1951). The 
pulses into the coincidence circuit were clipped to 6 x 10~° sec using a 
shorted transmission line, space charge limiting in the multiplier ensuring 
that the pulses reached their limiting value within this time. 


§ 3. EXPERIMENT 


Separated 7Li targets of 350 g/cm? superficial weight deposited on a 
silver backing were bombarded with a beam of 450 kev protons. as 
The bias of the discriminator in the «-counter were set so as to reject 
the large number of small pulses due to x-rays from the target and the 


4B2 
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timing unit was arranged to stop the scaling circuits after 10° y-ray counts 
had been recorded. ; 

The coincidence counts recorded in this interval are plotted in fig. 3 
against a value for the equivalent range which includes the range in the 
crystal required to produce a pulse of sufficient amplitude to be detected. 
«particles from polonium were used to calibrate this scale. ; 

The minimum absorber which could be used was 0-9 cm air equivalent 
because with less absorber, pulses due to the soft x-rays and protons 
scattered from contaminants on the target caused excessive currents in the 
photo-multiplier tube. The chance coincidence counts were about one- 
tenth of the true coincidences. It was necessary to check the balance of 
the fact coincidence circuit frequently, since the pulses producing the 


Fig. 3 
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Showing the range of the alpha-particles detected in coincidence with a 
gamma-ray. 


coincidences were small compared with the background pulses and, there- 
fore, small pulses from the coincidence circuit had to be counted, although 
a slight misbalance could give such an output from a single large back- 
ground pulse. Each point on the curve is the average of four observations 
of 100 seconds duration, during which time 10° gamma-rays and some 
1-5 x 10° alpha-particles were counted. About 95° of the alpha-particles 
were produced by the reaction *Li(pa)*He and had a range of about 
9cm. Hence the alpha counting rate and therefore the chance coincidence 
rate changed but little over the range of fig. 3. In spite of the checks on 
the fast coincidence circuit there was a spread in the coincidence counts 
for any particular absorber. 

Counts were therefore taken for a sequence of absorbers and this was 
reversed and repeated several times. Snedecor’s F test applied to the 
many readings taken at each point showed that the standard deviation 


of the readings was not significantly greater than that expected from a 
Poisson distribution. 
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§ 4. Eneray Levers or 8Be. 


Three groups of alpha-particles, having ranges of 1:19 cm, 1-60 em, and 
2-50 cm were detected in coincidence with a gamma-ray. These could 
only have been produced by a cascade reaction of the type 

"Li+ *H>*Be*—8Bet+y+Q,, et. ee) 

Ue Sd si 2a se I te 8 ae eer) 

Devons and Hine (1949) have confirmed the expected lifetime of about 
10°'° seconds for the 17-63 Mev level of *Be, formed by the resonant 
capture of the proton. Therefore, in these experiments, before the ’Be 
disintegrated as shown in (2) it was in motion due to the combined effect 
of the momentum of the incident proton, plus the recoil from the emitted 
gamma-ray. Since in the experiment the gamma-ray was detected in the 
opposite direction to the alpha-particle, the resultant motion of the nuclei 
which produced coincidence counts was uniquely determined. In order 
to calculate the levels of §Be involved, the measured energy of the alpha- 
particles was converted to the equivalent energy in the centre-of-mass 


Fig. 4 
Direction of incident 
Beam 


a- particle 


The momentum diagram for an alpha-particle detected at 90° to the direction 
of the beam and at 180° to the direction of emission of the associated 
gamma-ray. 


system. The gamma-ray energy was then equal to the difference between the 
energy of the *Be level from which the «-particles arose and the 17-63 Mev 
of the compound nucleus. Thus the energy conversions were simplified 
by assuming the energy of the level to be twice the measured energy of the 
alphas and so obtaining a value of the gamma-ray energy which was 
sufficiently accurate to calculate the recoil. The momentum diagram, 
fig. 4, represents the momenta for one of the alpha-particles before and 
after disintegration. The recoil energy received from the proton was 
0-055 Mev. The momentum gained from the proton by one of the alpha- 
particles in the combined nucleus was represented by AB and was 
proportional to $4/55. The recoil momentum from the emission of the 
gamma-ray was represented by AC (about 44/10), and the resultant 
momentum by AD. Then if AE was the momentum of the alpha-particle 
detected, DE was the momentum of the particle in the centre-of-mass 
system. 

In table 1 the measured range of the alpha-particles is shown. The 
range-energy relations found by Bethe (1950) were used to obtain the energy 
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values shown in column 2. These include an extra 40 kev to compensate 
for loss in the target. The energy in the centre-of-mass system 1s also 


listed. So 
The energy released, Q», in the dissociation of the *Be nucleus with an 


excitation energy H above the ground state is 
Q,—H-+0-09 Mev, 


the 0-09 mev being the energy difference between the ground state of *Be 
and two free alpha-particles (Ajzenberg and Lauritsen 1952). Allowance 
was made for this when calculating the excitation energies shown in the 
table. 


Table 1 
Energy of alpha- —_ Energy of alpha- Corresponding 
Measured particles. particles. excitation 
range Lab. system C. of M. system energy of §Be 
(cm) (Mev) (Mev) (Mev) 
1:19+0-1 2-36 +0-16 2-20 4-32 +0-32 
1-60 +0-1 2-96 40-14 2-79 5:50 +0-28 
2-50 +0-1 4-40 +0-10 3-91 7:73 +0-2 


The accuracy of measurement of the excitation energy was sufficient 
to identify the source of the 1:19cm and the 2-50 cm groups as the 
4-09 Mev and 7-5 Mev levels of *Be. The 1-60 cm group confirms the 
existence of a level which was suggested by Titterton et al. (Calcraft and 
Titterton 1951, Brinkworth and Titterton 1951, Titterton 1953). The 
greatest uncertainty in the value of the excitation energy, and consequently 
the energy of the associated gamma-ray, is due to the difficulty in 
determining the zero for the range measurements. The relative ranges of 
the groups are about four times more accurate. For this reason the energy 
values quoted in table 2 were based on the value 4-09 +0-018 Mev deter- 
mined by Salmon (1951) for the energy of the level of SBe which apparently 
disintegrated into the alpha-particles observed with a range of 1-19 cm. 


Table 2. Corrected Energy Values 


Corrected Energy of Corresponding Gamma-ray 
range particle excitation of ®Be energy 
(cm) (Mev) (Mev) (Mev) 

1-12 2-24 4-09 (assumed) 13-54 
1-53 +0-02 2-87 +0-03 5°31 +0-06 12-32 +0-06 
2-39 +0-02 3°95 + 0-02 7-51 +0-06 10-12 +0-06 


Using this calibration the alpha-particles from the well-known level at 
2-9 Mev would have been centred at 0-86 cm range in fig. 3, a fact which 


suggests the cause of the rise observed in the coincidence counts when the 
thinnest absorbers were used. 
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§ 5. INTENSITY OF Gamma-RADIATION 


The number of coincidence counts © detected from n transitions each 
be of which produces a gamma-ray and two alpha-particles would be 


C=e, X 26, Xn 


where «,—the absolute efficiency of the gamma-ray counter, assumed to 
constant for gamma-rays in the energy range 10 to 18 Mev and 


€,—absolute efficiency of the alpha-counter 
—=1/138. 


If the n transitions were accompanied by N gamma-radiating events, the 
number of counts JV, in the gamma-ray counter would be 


N,=(N-+n) ¢,. 


The fraction n/(N-+n)=C/(2e, N,) then expresses the intensity of the 
-weak gamma-ray component relative to the total. Table 3 lists the 
intensities, expressed as a percentage, calculated in this way from the 
data of fig. 3. 


Table 3 
Energy of 
radiation Intensity 
(Mev) (%) 
13-54 2:0+0°3 
12-2 1-9+0-3 
10-1 0-5 +0-1 


The values quoted for the intensity and the accuracy depend upon the 
assumption that all the transitions to the levels in *Be which emitted 
products into the cone of acceptance of the counters, were detected as 
coincidence events. Although this is likely to be true it is not possible 
to be sure of it. It is very unlikely that the intensities could be as much 
as twice those observed, unless the angular correlations were markedly 
anisotropic. 


§ 6. OPERATION OF THE Fast CoINCIDENCE CIRCUIT 


Because such a small fraction of the counts in each of the detectors was 
registered as coincident events, it was necessary to ensure that these 
could not have been produced by faulty operation of the coincidence 
circuits. 

Figure 5 shows the results of three tests which were applied. The 
coincidence rate is plotted against the delay produced by extra cables 
inserted in the y-counter channel. 

Curve (a) was obtained with a ®°Co source placed between the counters 
and gives a value for the delay caused by the cable used in the v-channel 
being longer than that used in the y-channel, namely 11 mus. 
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Curve (b) was obtained for the y-rays and «-particles from the reaction 
under study. The 5 mys displacement of the peak of (b) with respect to 
(a) compares well with the value 4-4 mys calculated as the time required by 
a 2-5 Mev a-particle to travel to the crystal. The asymmetry of the curve, 
it being wider to the left of the peak, would be due to the shorter time of 
flight of the more energetic «-particles. 

Leaving the alpha-counter in place, the gamma counter was removed 
from its position near the target chamber and a Co source was placed 
near to it to make the Co gamma-counting rate equal to that produced 
by the gamma-rays from the 7Li(py)*Be reaction. The E.H.T. supply 
to the photo-tube was then increased until the output pulses were of the 
same amplitude as those produced by the gamma-rays from the *Be. 


Fig. 5 


Coincidence Counts per 10° Y¥-ray Counts. 


Delay in Pray channel - myuwsec 


The change in coincidence counts produced by varying the delay in the gamma- 
ray channel, with a fixed delay of 11 ms in the alpha-particle 
channel. 

Curve (a) for gamma—gamma coincidences of gamma-rays from Co, 
(6) for alpha-gamma coincidences of radiation from *Be*, 
(c) for alpha~gamma coincidences of gamma-rays from ®°Co and 
alpha-particles from *Be*. 


> 


The amplitude of the pulses and the number entering the coincidence 
circuit were then substantially the same as for the main experiment. In 
this case no more than the expected chance coincidences were observed 
(curve (c) fig. 5) although there was a small reduction in the rate for delays 
greater than 25 mus, probably due to unexpected losses introduced by the 
25 mus delay cable. 

These experiments taken together confirm that the coincidences observed 
from *Be were due to alpha-particles and gamma-rays from the target and 
not to some spurious effect in the circuits. 
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§7. THE AupHA-PaRTICLE SPECTRUM AND THE PULSE AMPLITUDE 
SELECTION 


Having identified the y-« transitions by coincidence measurements the 
pulse amplitude selector was set to record only those pulses within a limited 
amplitude range so that the «-counter pulses produced by the 9 Mev 
particles from the reaction 7Li(px)He were rejected. The absorption 
curve obtained is shown in fig. 6. It was not possible to resolve the 
a-particle groups from the *Be levels because of the background from the 
pile-up of x-ray pulses and pulses produced by protons scattered from 
target contamination. This background varied from point to point across 
the surface of the target so that reproduceable counts could not be 
recorded. According to the coincidence counts, about 5000 of the 23,000 
counts shown at 1-0 cm in fig. 6 would have come from the disintegration 
of the excited *Be nuclei produced by the y-ray transitions. 


Fig. 6 


WO” Ge ray counts. 


Alpha -particle counts per 


10 15 20 | 2°5 cm. 
Range in air equivalent to fofol absorber. 


Pulses from the alpha-particle counter within the amplitude limits of 25 to 
165 volts. On this scale the amplitude of the pulses due to the alpha- 
particles from *Li(pa)*He was 360 volts. 


§ 8. Discussion oF RESULTS 


The work of Walker and McDaniel with the pair spectrometer would not 
have revealed the presence of lines occurring with an intensity as low as 
those reported here. On the other hand in the photo-disintegration 
experiments discussed above, the 13-54 and 12-32 Mev lines would be 
difficult to resolve and could be expected to appear as a broad line at 
about 12:8 mev. The total intensity would have been estimated as 
7°%, of the 17-6 Mev line. These figures are in reasonable agreement with 
the values 12-6 mev and 10% found by Nabholz et al., and 12-5 Mev and 
5% quoted by Titterton. They are, however, in marked disagreement 
with the high intensity (50°) reported by Goward and Wilkins. However, 
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the relative intensities would be expected to depend upon the beam energy 
and the angle between the direction of observation. Measurements have 
been made at a beam energy of 700 kev and others with an angle of 120° 
between the counters. Each confirm that the relative intensities do change 
but not enough to explain the result of Goward and Wilkins. Measure- 
ment of the y-« angular correlations should make it possible to determine 
the angular momentum quantum numbers for the levels concerned and 
enable a test to be made of the alpha-particle model for Be. 
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ABSTRACT 


Large-angle scatterings which have been found on 126-6 m of the tracks 
of relativistic particles in Ilford G5 plates exposed at a depth of 60 m.w.e. 
underground, are analysed, and are shown to be scatterings of j.-mesons. 
A cross section of (2-3-+1) x 10-?8 cm?/nucleon is derived for scatterings 
greater than 7° by u-mesons with pB>100 mev/c in the cosmic radiation 
at 60 m.w.e., and a cross section of (3-4-+1-5) x 10-2? cm2/nucleon for a 
similar scattering by -mesons with pf in the region (100-600) Mev/c. 
These cross sections are shown to be rather too high to be explained in 
terms of a pure Coulomb scattering at a finite nucleus. Agreement 
between theory and experiment can be achieved, however, on the 
assumption of a point charge nucleus. Various mechanisms are examined 
for a possible explanation of the effect without success, and it is concluded 
that the scatterings probably represent a weak nuclear interaction of 
po-mesons. 


§1. INTRODUCTION 


It is well known that a slow p-meson has only a very weak interaction 
with matter. Wheeler (1949) has shown that the «-meson is capable of 
travelling distances of the order of 200 cm in nuclear matter before being 
captured. However, during the past thirteen years, evidence has been 
presented by various authors for the scattering of .-mesons as they passed 
through metal plates, in a manner which could not be explained in terms 
of Coulomb forces and a finite nucleus (Wilson 1940, Code 1941, Sinha 
1945, Shutt 1946, Sahiar 1951, George and Trent 1951, Amaldi and 
Fidecaro 1950, Whitemore and Shutt 1952, Leontic and Wolfendale 1953). 
The experiment described in I of the present series of papers provided a 
good opportunity to investigate this problem further using the photo- 
graphic plate technique. As shown in I, the observations refer to 126-6 m 
Seen Oeigee tae ae en une tee ee ee ye 
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of straight minimum ionization track occurring in [ford G5 plates, 
coated, exposed, and processed at a depth of 60 m.w.e. underground. 
In I and II, the knock-on electrons, tridents, electron pairs, and stars 
found during the scanning of this length of track were analysed. In the 
present paper, the large-angle scatterings are described and discussed. 


§2. EXPERIMENTAL RESULTS 


A large angle scattering is observed as the sudden deflection of a straight 
track at a point (see I, fig. 2 (d)). A list of all such events recorded is 
given in table 1, where «° is the spatial deflection. No significant energy 
change is noticed at the deflection. 


Table 1 
Track before deflection Track after deflection 
kven¢)| eee ee Detlccion 
No. | Lengthmm | pf Mev/c | Length mm | pf Mev/c a? 
1 4 110+20 15 190 +35 5-7+0:-1 
2, 36 275 +20 16 280 +30 7-4+0-2 
3 13 550 +60 1-2 >300. 8:0 40-5 
4 16 130 +28 9 100 +20 8-1+0°5 
5 7-2 120+30 1-5 ~100 9-0+0-1 
6 5 210+40 16 140 +20 9-1+0-2 
7 20 340 +50 Td 210 +25 10-8 40:5 
8 17 160 +25 1-2 ~100 34:4-+40-1 


All the tracks were grain counted on both sides of the deflection, and 
none of the counts were significantly different from counts on high-energy 
electron pairs on the same plate. The results of the grain-counting are 
shown in fig. 1. The plates were not developed sufficiently uniformly to 
distinguish between * plateau ’ and minimum grain density. 


§3. IDENTITY OF THE PARTICLES 


From fig. 1, it is clear that of the known particles, the tracks could have 
been caused by 7-mesons, -mesons or electrons only. It will be shown 
that 7-mesons and electrons are very unlikely to have caused the events. 

As a result of a scan of 27 m of 7-meson track, with pf in the range 
(150-200) Mev/c, Bernardini e¢ al. (1951) gave a mean free path of about 
25cm in G5 emulsion for the production of stars, inelastic scatterings 
(where the grain density changes by a factor of more than three) and 
stoppings in the emulsion (cf. Lock 1952). None of these events which are 
characteristic of 7-mesons have been observed in the present experiment. 
Again, from the underground star rate given by George and Evans (1951), 
it can be shown that the percentage of 7-meson tracks in the relativistic 
radiation at 60 m.w.e. is only 0-05% (George, private communication). 


Jmin 
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Thus, of the total track length followed, viz. 126-6 m, only 6 cm would be 
expected to be due to z-mesons. Hence it is extremely unlikely that the 
large angle scatterings are due to 7-mesons. 

It was estimated in II that about 70 cm of the track length examined is 
due to electrons with energy greater than 100 Mev. It might appear at 
first sight that the large angle deflections could be caused by ‘ Bremsstrah- 
lung’ from these electrons. However, a 100 Mev electron is already a 
highly relativistic particle, and a deflection of the order of 7° is unlikely 
from such a process unless the electron were to lose almost its entire energy. 
If the events are to be explained in terms of some other form of electron 
interaction, then the cross section for such a process would have to be 
almost twice the geometric cross section and no such interaction is known. 


Fig. 1 


PB Mev/c 
Grain density as a function of pf. 


The tracks will therefore be assumed to be due to j.-mesons. The possi- 
bility that the events are caused by the decay of u-mesons can be dismissed 
at once, since the time taken by a relativistic particle to travel 126-6 m is 
about 4:5 10-7 sec, which is shorter than the lifetime of u-mesons even 
without time dilatation. 


§4. Momentum TRANSFER 
The momentum transfer involved in the events described in table | is 
given in table 2, together with the corresponding energy transfer to a 
proton, and the range of such a proton in Hford G5 emulsion. 
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Since no recoil track was observed in any of the events, it is clear from 
table 2 that at least four of the deflections could not have been at hydrogen 


nuclei. 


Table 2 


_ ee eal 


Range of 
Energy transfer | proton in 


Event | Momentum transfer to a proton Tlford G5 


No. BE Mev emulsion 

in pu 

il 10 0:07 <l 

2 34 0-6 al 

3 70 2°6 60 

4 1 0-08 =<“, 

5 12 0-08 al: 

6 24. 0-2 1:5 

7 46 1-1 17 

8 66 2:3 50 


§5. SCANNING EFFICIENCY 


It is important to examine the geometrical factors which influence the 
finding of these events. It is clear that relatively small deflections will be 
largely missed, whereas large deflections will almost certainly be detected. 
If the primary track remains flat in the field of view before and after the 
deflection, it is estimated that a deflection of less than 3° would not be 
detected at all, and that the scattering would have to be greater than 4° 
before the event is observed with high probability. This is the most 
favourable case for detection. The most unfavourable case arises when 
the track deflects in a plane at right angles to the surface of the emulsion. 
In this case the only way of detecting the deflection is to notice the sudden 
change in the flatness of the track being followed. Now a track is followed 
in the first place because it appears flat in the field of view. From the 
track length distributions (I, fig. 1), it is seen that the number of tracks 
followed drops sharply when the track length drops below 5mm. Thus it 
is reasonable to assume that if a track which is flat at first suffers a sudden 
deflection in a plane at right angles to the surface of the emulsion, so as to 
cause it to move subsequent to the deflection with an inclination corre- 
sponding to a track of length less than 5 mm, then the deflection would have 
a large probability of being observed. Clearly, a deflection which causes 
a track moving towards the surface (say) to turn back towards the glass 
would usually be detected quite easily. Most of the deflections that 
actually occur will lie between the extreme cases described above. 

A track of 5 mm length has a steepness of about 6° with respect to the 
surface of an unprocessed 400 1 or 6004 emulsion. It follows then that 
deflections less than 6° would have an appreciable chance of being missed, 
whereas deflections larger than about 7° should be detected with high 
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efficiency. Hence in evaluating cross sections, only deflections greater 
than 7° will be considered. 


§6. Errective Track LENGTH 


When a track suffers a large deflection, it ceases in general to be flat in 
the emulsion. The consequent shortening of the track results in a smaller 
probability ofits being followed. This effect is largely compensated for by 
tracks being scattered into the emulsion. The effective track length can 
therefore be taken as the length actually followed—i.e. 126-6 m. It is 
unlikely that this estimate is in error by more than 25%. 


§7. THe EXPERIMENTAL Cross SEcTION FoR Larcu-ANGLE SCATTERINGS 


A mean free path of 18 m of [ford G5 emulsion is deduced for a deflec- 
tion greater than 7° by a u-meson with pf greater than 100 Mev/c in the 
cosmic radiation at 60m.w.e. This corresponds to a cross section 
o(a>7°)=(2-3-+1) x 10-28 cm?/nucleon. However, since all the events 
recorded are associated with low energy tracks having pf less than 
600 Mev/c, a cross section referring to -mesons of these energies gives a 
better picture of the situation. According to the u-meson energy spectrum 
at 60 m.w.e. given by I, eqn. (4), the fraction of tracks with pf in the 
region (100-600) Mev/c is about one fifteenth. A cross section 
o(a>7°)=(3-4+1-5) x 10-2? em?/nucleon is then obtained for a deflection 
greater than 7° in G5 emulsion by u-mesons with pf between 100 and 
600 Mev/c. 


§8. THE ExPECTED NUMBER OF LARGE-ANGLE SCATTERINGS FROM A PURE 
CovuLoMB INTERACTION 


If x(0) d@ cm?/nucleus is the differential cross section for the scattering 
of a fast charged particle (charge ze and momentum p) through an angle 6 
in the range d@ by a nucleus of charge Ze, then if the charge on the nucleus 
is concentrated at a point, the well known Rutherford formula applies : 


2nZ*z7e4 d(sin 0/2) 
See eee te 
x(9)d pB?c? sin®-Of2 — (1) 
The integral cross section for a deflection greater than «° is given by 
te mater 1 
Ha)=|_ x(0)d0= ppc Free ris |. ate Lins (2 a) 


If p is in Mev/c and «=7°, eqn. (2 a) becomes 

I(7°)=1-74 <10-73(Z2/p28?)-em*/nucleus. . . . (20) 
It may be noted that for Ilford G5 emulsion Z?=405 and the average 
number of nucleons per nucleus= 27. 

The finite size of the nucleus was taken into account by Williams (1939) 
by assuming the nuclear charge to be distributed uniformly over a sphere 
of radius 7,=1-5x 10-4‘ cm (the classical nuclear radius). Williams 
showed that the Coulomb potential inside the nucleus could then be 


represented by 
V(R)=Ze/R[1—exp (—2R/r,)]. 
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Using this potential distribution, it was shown that for a finite nucleus, 
the Rutherford formula becomes 

27Z2z?e4 d(sin 6/2) 
ppc? _—sin® 6/2 


, 1 pay ir pean tee 
Kaji FOF] Beatie 02 
The variation of K with 0/@,,.. is shown in table 3. 

Using eqns. (2b) and (3), the differential energy spectrum of the 
u-meson flux at 60 m.w.e. given in I eqn. (3), and also the composition of 
Ilford G5 emulsion, the expected number of large angle scatterings greater 
than 7° can be calculated. 


x (0)d0= Bo S ae 


where 


Table 3 
6/ Omax K 6/0 max K 
0:3 0-84 5 0-09 
0-5 0-64 2-0 0-04 
1-0 0:25 0-01 


Since only tracks which appeared to be at minimum grain density were 
followed, only mesons with pf greater than 100 Mev/c (or kinetic energy 
greater than 50 Mev) are of interest. The u-meson flux has therefore been 
divided into energy bands each extending over a range of 100 Mev from 
50 Mev upwards, and the expected number of large angle scatterings 
calculated from each of these energy bands. _The results are given in 
table 4. It is found that scattering at silver and bromine nuclei of the 
emulsion provide practically the entire cross section for the scattering. 
In fig. 2, O,,,. is plotted against pf for silver and bromine nuclei. 

It has been shown by Williams (1939) that the effect of a finite nucleus 
becomes apparent at scattering angles greater than h/pr, (i.e. 0,,./2). 

In the present experiment, five events were observed with «>/h/pr,. 
These are the points above the dashed line in fig. 2. The expected 
number of large-angle scatterings falling into this category is 5-5 for a 
point charge nucleus and 1-5 for a finite nucleus. A similar calculation 
can be made for scattering angles greater than 2h/pr,. In this case, the 
expected number for a point charge nucleus is 2-1, and 0-12 for a finite 
nucleus. The observed number of events was 2. These results are 
summarized in table 5. 

It can be shown, Regener (1951), that the probability that the events 
are solely due to Rutherford scattering at a finite nucleus of the Williams’ 
type (with r,=1-5x 10-8 4-+), is less than 5%. If the assumed value of 
r, 18 too large by 25%, the probability is increased from 5%, to 16%. 

Defining an anomalous scattering as one which cannot be explained in 
terms of Rutherford scattering at a finite nucleus, a cross section 
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o(anomalous)=(1-5-++1) x 10-2? cm2/nucleon is deduced for j.-mesons with 
pp between (100-600) mev/c. 


Table 4 
Observed Expected number of 
Range of Corresponding number of scatterings >7° 
K.E. range of pB scatterings §_—————_——_—_—___ 
Mev Mev/c aah Point charge} ‘Finite 
nucleus nucleus 
50-150 90-220 4 6 3 
150-250 220-340 2 1-6 0-5 
250-350 340-450 0 0-8 0-1 
350-450 450-560 1 0:5 0-05 
450-550 560-650 0 0-35 0-02 
550-650 650-750 0 0-25 — 
650-750 750-850 0 0-2 — 
750-850 850-950 0 0-15 — 
850-950 950-1050 0 0-1 ~- 
Total 7 10 3°7 
Fig. 2 


iE 


IO 


pB. MeV/c, Es 


6... and Omax/2 (full and dashed lines respectively), as a function of pB for 
a0 bromine and silver nuclei, 
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$9. ResuLtts oF OTHER WORKERS 


Recent cloud-chamber experiments by Whitemore and Shutt (1952) 
and by Leontic and Wolfendale (1953) indicate clearly the tendency for 
the large angle scattering of »-mesons in lead to fit the calculations better 
when no cut-off is applied for finite nuclear size. The earlier work on this 
subject is summarized in the paper of Leontic and Wolfendale (1953). 
If this picture of a point charge nucleus has any significance, the cross 
section per nucleon for a large angle scattering should be proportional to 
Z and hence should be lower in a light material. Between lead and [ford 
G5 emulsion the expected difference in the cross section per nucleon is only 
a factor two, and the statistical accuracy of the experimental results is not 


Table 5 


Expected number 
Observed No. 


Range of scattering angle of events Point charge Finite 
nucleus nucleus 
vis 7 10 a 
=T-Land L=nlpr, 5 5:5 1-5 
>7° and >2h/pr, ps 2-2 0-12 


sufficient to make this difference significant. Recently Walker (1953), 
working with -mesons of energy greater than 1-5 Bev at sea-level, has 
failed to observe a single deflection greater than 5° at a one inch carbon 
plate after 21.500 traversals. He concluded that the cross section for 
scattering greater than 5° in carbon is less than 2x 10~?° em?/nucleon. 
The expected value for carbon for a point nucleus is 10-®° em?/nucleon, 
For lead one would expect a cross section of the order of 


5x 10-79 em?/nucleon 
for a deflection greater than 5° by u-mesons in the stated energy range. 


Thus Walker’s results appear to give a slight indication that the scattering 
is Z dependent per nucleon (or Z? per nucleus). 


$10. PosstsLE EXPLANATIONS OF THE ANOMALOUS SCATTERING 


In the analysis presented in §8, the effect of a finite nucleus was cal- 
culated by considering the nuclear charge Ze to be spread uniformly over 
the volume of the nucleus. Such a model can strictly be used only when 
the de Broglie wavelength of the incident particle is longer than the mean 
distance of protons inside the nucleus. When this wavelength is of the 
same order of magnitude or shorter than such a distance, as it is in the 
case with the particles under consideration, the nuclear protons have to be 
considered as individual centres of Coulomb forces distributed inside the 
nucleus. One has therefore to consider the effects of coherent and 
incoherent scattering within a nucleus: in the former case the nuclear 
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protons are not excited, whereas in the case of incoherent scattering the 
nuclear protons undergo a transition from the ground state to an excited 
state. 

A simple approach to the problem of coherent scattering within a 
nucleus could be made by first giving each of the charged nucleons in a 
complex nucleus the same status as a point charge proton. Then ‘using 
eqn. (26) and the fact that the number of charged nucleons per cm®* of 
G5 emulsion is approximately 1-8N, (where Ny is Avogadro’s number), 
it can be shown that the expected number of large angle scatterings greater 
_ than 7° in the present experiment is only 0-3, a value at least 10 times 
smaller than the observed number of anomalous events. It is difficult to 
visualize a plural scattering mechanism within a nucleus which would 
increase the cross section by such a large factor. 

The effects of incoherent scattering have been considered in some detail 
by Amaldi et al. (1950 a, b) for light nuclei only. The problem was shown 
to be similar to that of the scattering of an electron of a few kev energy or 
x-rays by anatom. In particular, Amaldi e¢ al. considered the scattering 
of 200 Mev and 600 Mev yz-mesons at lithium and carbon nuclei, using two 
possible extreme values of the electromagnetic radius of the proton. They 
deduced cross sections for large-angle scatterings which were of the same 
order of magnitude as those calculated on a point charge picture. How- 
ever, the incoherent scattering cross section is only Z dependent per 
nucleus. Thus if the cross sections of Amaldi et al. can be extrapolated 
to heavy nuclei such as silver and lead, the cross sections for incoherent 
scattering would be expected to be considerably smaller than the corre- 
sponding cross sections for Rutherford scattering worked out on a point 
charge model. Such an extrapolation has been made by Whitemore and 
Shutt (1952), who estimate that for lead the cross section drops to one- 
tenth the value calculated for a point nucleus when 7, the electromagnetic 
radius of the proton, is taken as 2-1 x 10-14 cm, and to one-hundredth the 
value for a point nucleus when 7, is taken to be 1-4 x 10~¥ cm. 

It would therefore appear that the explanation of the phenomenon is 
not to be found in terms of coherent or incoherent Coulomb scattering 
within a nucleus. 

Another possible explanation of the effect is in terms of a photo-nuclear 
interaction between the ‘ virtual’ photons of the electromagnetic field 
of the incident -meson and the short-range nuclear forces. George and 

‘Evans (1950) explained the stars produced by fast .-mesons underground 
in terms of such an interaction. They deduced a cross section, o, for an 
energy transfer greater than « of the form o=Ko,,{log,(#/<) ? cm/nucleon, 
where H is the energy of the incident .-meson, K is a constant, and o,, 
is the cross section for nuclear disintegration by a high-energy photon. 
Substituting H=14 Bev for the average energy of y-mesons at 60 m.w.e., 
o—10-28 em?/nucleon from the experimental work of Miller (1951) and 
¢=150 Mev as the energy required to produce a 3-prong star, they obtained 
a cross section o=4 x 10-39 cm2/nucleon. In the present experiment, 
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the average energy of the u-mesons under consideration is about 200 Mev, 
and the energy transfer to a nucleon is about 2 Mev. Thus (E/e) remains 
the same as for George’s calculations. However, the experimental 
cross section for anomalous scattering is o~1-5 x 10~°? cm?/nucleon. 
Hence to bring experiment and theory into agreement o;,, would have to 
be increased to about 3X 10-26 em2/nucleon for small energy transfers of 
the order of 1 or 2 Mev, as against 10-8 em?/nucleon for larger energy 
transfers of the order of 150 Mev. There is no evidence whether such an 
assumption is reasonable or not. This mechanism is Z-dependent per 
nucleus, so that if experiments with light and heavy elements were to 
support a Z?-dependence, this process would be excluded. 


§11. CoNcLUSION 


The large-angle scatterings of .-mesons are shown to be consistent with 
theory if the scattering nuclei are imagined as point charges. Since such 
a picture is contrary to present theories of the structure of the nucleus, 
the agreement seems to be accidental. It is possible that the scatterings 
represent a new type of nuclear interaction. 


ACKNOWLEDGMENTS 


This investigation was carried out under the direction of Dr. G. D. 
Rochester, whose valuable advice is deeply appreciated. The authors 
thank Dr. A. W. Wolfendale and Mr. B. Leontic for communicating 
their results prior to publication. The authors’ thanks are also due to 
Dr. G. N. Fowler and the members of the photographic emulsion 
group at Manchester for helpful cooperation and discussion. Mr. P. H. 
Fowler and Dr. M. G. K. Menon of Bristol University have made valuable 
comments on the paper. 


REFERENCES 


AMALDI, E., and Frpgcaro, G., 1950, Nuovo Cim., 7, 535. 

AMALDI, E., Frpecaro, G., and Martant, F., 1950.a, Nuovo Cim., 7, 553; 1950 b, 
Lode Tore 

BERNADINI, G., Boora, E. T., Leprerman, L., and Tintort, J. H., 1951, Phys. 
Rev., 82, 105. 

Cove, L., 1941, Phys. Rev., 59, 229. 

GEorGE, E. P., and Evans, J., 1950, Proc. Phys. Soc. A, 63, 1248. 

GrorGE, KE. P., and Trent, P. T., 1949, Nature, Lond., 164, 838. 

Lxontic, B., and WoLrenbaLE, A. W., 1953, Phil. Mag., 44, 1101. 

Lock, W. O., 1952, Ph.D. Thesis, Bristol University. 

Miter, R. D., 1951, Phys. Rev., 82, 260. 

REGENER, H. V., 1951, Phys. Rev., 84, 164. 

SanzaR, A. B., 1951, Proc. Indian Acad. Sci., 34, No. 3A, 201. 

Suutt, R. P., 1946, Phys. Rev., 69, 261. 

Srnwa, K., 1945, Phys. Rev., 68, 153. ‘ 

Watker, W. D., 1953, private communication. 

WHEELER, J. A., 1949, Rev. Mod. Phys., 21, 133. 

Wurrteworg, W. L., and Suurv, R. P., 1952, Phys. Rev., 88, 1312, 

WiitaMs, KE. J., 1939, Proc. Roy. Soc. A, 169, 531. 

Witson, J. G., 1940, Proc. Roy. Soc, A, 174, 73. 


Le | 


CXV. The Anomalous Scattering of u-Mesons in Lead 


By B. Lrontic and A. W. WoLFENDALE 
The Physical Laboratories, The University, Manchester* 


[Received April 25, 1953] 


ABSTRACT 


A multi-plate cloud chamber has been operated at sea-level under a 
thick lead absorber and the scattering of penetrating particles— 
presumably j.-mesons—has been investigated. In common with the 
results of other workers, a component of scattering has been found which 
is not consistent with the distribution expected from simple Coulomb 
scattering from a nucleus in which the charge is distributed uniformly 
throughout the volume (‘ solid’ nucleus). We find that the magnitude 
of this anomalous component is similar to that expected from the Coulomb 
interaction with a nucleus of effective radius much smaller than that of 
the ‘solid’ nucleus. Since it is very unlikely that the electrostatic 
charge of the nucleus can be considered to be concentrated at the centre 
of the nucleus, we conclude that the anomaly is probably due to a short 
range non-electric interaction between p-mesons and nucleons. 


§ 1. INTRODUCTION 


THE interaction of very slow p-mesons with matter has been studied 
by a number of workers and has been found to be very weak (Conversi 
et al. 1945, 1947, etc.). In the present work the interaction of fast 
j-mesons with lead nuclei is found to be relatively large. 

Data on the scattering of y-mesons in lead is analysed statistically 
and the magnitude of a scattering component, additional to that expected 
from simple Coulomb scattering from a solid nucleus, is determined. 
The results of other workers are analysed using the scattering distribution 
found by us and are found to be mainly in agreement. 

The scattering distribution should give evidence for the particle- 
structure of the nucleus or the short range forces between ;.-mesons and 
nucleons. 


§ 2. EXPERIMENTAL ARRANGEMENT 


The experiment was carried out in the laboratory at sea-level using a 
cloud chamber of illuminated volume 40x 25x25 cm? containing six 
2 em lead plates (fig. 1). The chamber was illuminated through a rear 
window by two flash discharge tubes placed at the foci of two parabolic 
mirrors, and stereoscopic photographs were taken by twin cameras, the 


axis of photography making an angle of about 20° with the direction of 
IRs eg pales st i eked ok a eves re ae 
* Communicated by Dr. G. D. Rochester. 
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illumination. The expansion for the chamber was produced by the 
symmetrical lateral movement of two rubber diaphragms on opposite 
sides of the chamber. Electromagnetic release valves, connected in 
series, controlled the movement of the diaphragms ; the delay between 
the passage of a particle through the chamber and the completion of 
the expansion was only 8 m sec, ensuring reasonably fine tracks. 

A narrow lead filter 100 cm in height was situated above the chamber 
and particles which had traversed the filter were selected by the counter 
telescope A BC. The telescope was situated wholly above the chamber 
in order to give no geometrical bias for events occurring inside the 
chamber. An upper limit to the number of primary and secondary 
protons and 7-mesons reaching the chamber has been calculated for an 
absorber of this thickness (Wolfendale 1953), and is quite negligible. 


Fig. | 


100 CM. Pb. 
FILTER 
~*~ 
. B 
CHAMBER c 
CONTAINING | = ee 
LEAD PLATES SSS Sa) 
| Sa ee | Ss OR 10 


r 
COUNTER Sia =D 
10cm 


SGACEQe at 


, 
The experimental arrangement. 


For part of the experiment an absorber 3 was present underneath 
the chamber, and below it there was a double tray of counters. This 
tray was connected to an indicator lamp by way of the appropriate 
circuits so that the further penetration of mesons appearing in the 
chamber could be determined, 

For the analysis single tracks were selected whose positions indicated 


that the particles had penetrated the filter and were in fact the particles 
responsible for the expansions. 
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§ 3. EXPERIMENTAL RESULTS 


2969 photographs were taken, each showing at least one particle 
penetrating all the plates. The projected angles of scattering of selected 
particles in some of the lead plates were measured using a microscope 
fitted with a goniometer eyepiece. On account of track distortion above 
the top and below the bottom plates only the deflections in the centre 
four plates were measured. Measurements were limited to cases where 
the root mean square deflection for the four traversals was greater than 


Table 1 
hy No. Oy No. 
2=3° 69 8-9° 3 
3-4° 44 9-10° 2 
45° 18 10-11° 0 
DOR 12 25 ] 
6=7° 16 TP2=13° 


3 

Total 174. 

2°, since scattering angles appreciably less than this magnitude were 
uncertain due to distortion. The selection of tracks to be measured 
was made by eye and a number of tracks were probably missed having 
root mean square projected deflections 0) >2°. However, it will be shown 


later that this inefficiency is not important. The frequency distribution 
for 0) is given in table 1. 


§ 4, ANALYSIS OF THE EXPERIMENTAL RESULTS 
4.1. Data with the Chamber and Indicator Tray 


In order to detect any anomalous scattering, the Coulomb scattering 
spectrum must be known and this can be calculated if the momentum 
spectrum is known from independent observations. This, however, is 
not the case and it is not possible to derive a Coulomb scattering spectrum 
by simple methods from the observed scatterings of table 1, since they 
are few in number. Consequently no attempt has been made to derive 
the magnitude of the anomalous scattering in this way from the complete 
data. 

625 photographs were taken in which the single meson seen in the 
chamber penetrated a further 5 or 10 cm Pb, and discharged a counter 
in the tray D below the chamber. Six photographs showed large-angle 
scattering of more than 10° in one of the centre four plates, and in each 
case the momentum of the meson after the interaction was > 250 Mev/c 
(from the known minimum range). In this momentum region the 
multiple Coulomb scattering has a fairly low average value, compared 
with the occasional large scatterings which may or may not be anomalous. 
It is therefore possible to infer an approximate multiple Coulomb 
scattering spectrum from our measurements, and on the assumption that 
this is Gaussian (Williams 1939, and others) we can compute the number 
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of large scatterings to be expected and compare it with the observed 
number. 

The effect of occasional large single scattering will, in general, be 
superposed on the multiple scattering distribution. This single scattering 
follows a power law in 0, and if the nucleus has a finite size it will be cut 
off at a maximum angle given by 


6, a (Williams 1939, Olbert 1952) . . . . (1) 


where 27A is the De Broglie wavelength of the meson and d=0-5 ro Als 
is the ‘radius’ of the nucleus. A is the atomic weight of the element 
and 7, is the classical electron radius. The true root mean square 
projected deflection 6,5 due to multiple scattering in an absorber of 
thickness ¢ is given by 


14:8 : : 
Oamus = IS \/t radians (Rossi and Greisen 1941, and others) . (2) 


where p is in Mev/c and ¢ is in radiation lengths (1 radiation length in 
lead=0-52 cm). 
From eqns. (1) and (2) we find 
Go ae ee 
Orus vt ®) 


For a 2 cm lead plate the maximum angle for single scattering is nearly 
the same as the RMS angle of multiple scattering and the contribution 
of single scattering to the distribution is very small. (In general, when 
the cut off is not considered, the single scattering ‘ tail’ is only important 
for deflections greater than about 3 Oxys.) The resultant scattering 
distribution in this experiment should, therefore, be approximately 
Gaussian. 

On the multiple Coulomb scattering distribution we expect only about 
one scattering having a deflection greater than 10° and five events remain 
as ‘anomalous’. This yields a cross section for anomalous scattering 
for angles greater than 10° with our particular spectrum 


O 49°~(15-0-7) 10-8 em?/nucleon, 


Table 2 gives the deflections, minimum ranges and corresponding 
minimum momenta of the particles scattered through more than 10°, 
together with the calculated probabilities that they form part of the 
multiple Coulomb scattering distribution (MCS). 

It will be noticed that the number of events decreases rapidly as the 
deflection @ increases so that the cross section for 06> 20°, say, is possibly 
an order of magnitude less than the corresponding cross section for > 10°. 


4.2. Statistical Analysis 


An alternative and more satisfactory treatment of the data has been 
carried out using all the measured tracks and considering the relative 
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magnitudes of the four deflection measurements made on each particle 
traversing the cloud chamber. The analysis is statistical in character 
and the method was suggested by Professor M. S. Bartlett. 


Table 2 
eee ee eee were ea ee SIME gee ee 
Probability of 


Event | Deflection @ | Min. range | Min. momentum : : 
deflection being 


No. degrees in cm Pb in Mev/c due to MCS °% 
] 18 13 285 0-4 
2 Al 18 345 10 
3 12°5 18 345 3 
4 Hal 11 250 40 
5 10 14 305 20 
6 in 13 285 20 


The main problem in the interpretation of the results consists in 
determining to what extent the anomalous component (non-Gaussian 
or Gaussian of different magnitude) contributes in producing occasional 
large-angle deflections. 

We first consider the case where no anomalous component exists, the 
scattering being simple Gaussian. 

Let 6)?=mean square angle from the four observations on a single 
particle ; ' 

6,,2—=Maximum square angle. 

We define then g=6,,7/4 6,7. 

Cochran (1941) has discussed the distribution of g and has given as 
the probability distribution of g 


8 
Polg)=-g'"%(1—g)'" forge}... . - - (4) 


The number distribution is then 
NAg)=NoP 9), 


where J, is the total number of events (for g>0). 

The advantage of this method of analysis is that g is independent of 
the true mean square and also of the observed mean square (for Gaussian 
observations). 

It is found that the distribution is independent of selection of 6)? so 
that rejecting small values of 6,2 will not affect appreciably the anomalous 
samples since these show up as large angles. 

The values of g have been calculated for the tracks having 6)>2° and 
the experimental distribution is given in fig. 2. The theoretical distribu- 
tion calculated from eqn. (4), normalized to the same number of particles 
having g>0, is drawn for comparison. To determine the effect of track 
distortion the distribution of g for 6,>4° has also been determined and 
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normalized to the same number of particles. This is also given in fig. 2. 
The form of the distribution is very similar to that for ,.>2°. It can be 
assumed therefore that the distortion has little effect on the shape of 
the distribution for either category, since we expect the distortion effects 
for 0)>4° to be quite negligible. es 

The discrepancy between the expected and experimental distributions 
is large particularly in the region of high g and the statistical x* test 
shows that the two distributions are significantly different, so that an 
anomalous effect clearly exists. 


Fig. 2 


25 


THEORETICAL 
aes NO ANOMALY 


EXPERIMENTAL | 


INTERVAL 


NO. PER 0-05 


ie) 


x 


Comparison of the experimental distribution of ‘g’ with the theoretical 
distribution for multiple Coulomb scattering only. 


The probability of anomalous scattering can be estimated on the 
following assumptions : 

(a) the chance « of an anomalous deflection is small, so that the 
probability of more than one such deflection in four observations is 
negligible ; 

(6) the distribution of the anomalous scattering is Gaussian with a 
larger standard deviation oj, which bears a constant ratio to the normal 
Gaussian distribution o, for particles of different momenta. 


on" 1 
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The latter assumption means that the probability for anomalous 
scattering varies with momentum in an identical manner to that of 
multiple and single scattering. Since the excess of events having large 
values of g occurs for all measured values of 5, i.e. for all momenta less 
than about 1 Bev/c, assumption (b) is at least approximately true. 

Assume further that an anomalous deflection provides the maximum 
reading (this is only strictly true in the limit as p+0). The distribution 
of g with one anomalous observation in each set of four deflections can 
be determined. 

This has been evaluated by Bartlett as 


zal Geet 


20S aan eet Betas aa (8) 


The chance of one anomaly in a set of four is 4e and so the resultant 
probability distribution of g is 


P(g)=(1—4e) Pag) + 4eP ug) 
and the number distribution 
N(g)=N P(g). 
By varying the parameters « and ., the curve of best fit (minimum ;?) to 
the experimental distribution (from fig. 2, 0)>2°) has been found and is 
_ given in fig. 3. 
The optimum values of the parameters are 

«=0-083, 

je==0-07. 
Since the optimum value of ¢« is small, assumption (a) is valid, A large 
number of pairs of « and u were tried and y* determined for each pair. 
A set of pairs that gave a constant y? corresponding to a 90°, chance of 
the true value lying within the set, were taken and the limiting 
distributions determined. 

We identify the multiple Coulomb scattering with the normal Gaussian 
distribution. The distributions are given in fig. 4. The abscissa 6/@x us 
is independent of the momentum of the particles for single scattering 
so that the probability distributions for a point nucleus (i.e. omitting the 
cut off mentioned on p. 1103) can also be plotted. Of the single scattering 
distributions given in fig. 4, one is that calculated by Williams (1939) 
and the other is the result of Snyder and Scott’s (1949) more detailed 
theoretical investigation. 

Our experimental distribution shows that the large-angle scattering 
observed for particles having 6)>2°, i.e. momenta less than about 1 Bev/c, 
is definitely in excess of the Coulomb scattering for a solid nucleus and 
of magnitude similar to that for a point nucleus. Our curve of best fit 
lies rather above that calculated for single scattering from a point nucleus, 
but this is not necessarily important because our experimental points are 
raised by the inclusion of small distortions. 
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The effective radius of the lead nucleus for this scattering may ie ue 
order of 1/10 of the value given by the usual relation ir 0-50 fe he ee 
still give results consistent with our experiment within the statisti 


errors. ae 
Examples of anomalous scattering are shown in Plates 44 an 


Fig. 3 


INTERVAL 
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Comparison of the experimental distributions of ‘g’ with the ‘ best-fit’ 


distribution, considering both multiple Coulomb scattering and anoma- 
lous scattering. 


§ 5. DETERMINATION OF THE Cross SECTION FOR ANOMALOUS 
SCATTERING 


If the so-called anomalous scattering is treated as a type of single 
Coulomb scattering, then the cross section will in general depend on the 
particular particle spectrum and the minimum angle of scattering accepted 
as anomalous. The cross section would in that case fall off rapidly with 
increasing angle of scatter, whereas if the scattering is due to short range 
nuclear interactions where the angle of emergence of a particle after 
interaction is probably symmetric in the centre of mass system, the fall 
off should be less rapid (Amaldi and Fidecaro 1950). 

From the best fitting curve (fig. 3) we can deduce the cross section for 
anomalous scattering as a function of the ratio of the minimum angle 
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accepted as anomalous to the RMS angle of multiple Coulomb scattering. 
The cross section can also be determined as a function of the minimum 


- momentum transfer in the anomalous scattering event, from the following 
considerations. 


Fig. 4 


INTERVAL 


UNIT 
io) 


NO PER 


SINGLE 
(POINT | 


Coulomb scattering distributions for ‘ solid’ and ‘ point’ nuclei compared with 
the distribution giving the ‘ best-fit’ to the experimental results. 


Consider a particle of momentum p (Mev/c) being scattered anomalously 
through a small angle 6. The momentum transfer for small deflections 
is approximately 
6 


Oaus 


0 
p=pl= perms Onis = 29-0 


from the equation 


14-8 

Oaus= =p Vt 

RMS DB V/ 
assuming that b~1. 


Therefore the momentum transfer is a function of 0/@gys and hence of 
the cross section, The approximate values are given in table 3, 
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Table 3 


Cross section (em2/nucleon) in | Approximate momentum 


4/Orats Pb for deflections>6 transfer Mev/c 
3 D3 6 One 85 
6 Bote aes 170 
9 6-9 x 10-29 255 


4-4 « 10~%° 


§ 6. COMPARISON WITH OTHER OBSERVATIONS 
The results of other experiments on anomalous scattering (AS) are 
given in table 4 where the cross sections (c) found are in most cases 
compared with the values expected from our results. It will be noticed 
that most of the results are consistent with our interpretation. 


§ 7. Discussion AND CONCLUSION 

It has been shown that most of the experimental results on «.-mesons 
scattering are consistent with the point nucleus hypothesis. Considerations 
of binding energy, nuclear stability, magnetic moment, etc., show, 
however, that the individual nucleons in a nucleus are not in fact concen- 
trated at a point but are distributed fairly uniformly over the nuclear 
volume. The effective radius of the nucleus for short range forces has 
been determined from neutron cross sections by De Juren and Knable 
(1950) and others. It is reasonable to assume that the effective radius 
for Coulomb interactions is of similar magnitude, and in fact the scattering 
of moderately fast electrons by nuclei (Lyman ef al. 1951) shows that 
the nucleus can be represented by a sphere of uniform electrostatic 
charge density having a radius about 20°, smaller than the value 
determined from the neutron interaction cross sections. 

It appears then that the similarity between the experimental results 
and those expected for the single Coulomb scattering from a point nucleus 
is quite fortuitous and that the anomaly is due to incoherent Coulomb 
interactions or to short range non-electric interactions. 

It is interesting to note here that the calculation of the coherent and 
incoherent Coulomb scattering of u-mesons by light nuclei, carried out 
by Amaldi et al. (1950), gives scattering distributions similar to those 
expected from coherent scattering by point nuclei. Whittemore and 
Shutt have estimated very approximately the maximum contribution 
from this incoherent scattering in the case of lead to be <10°% of the 
observed effect so that if the calculations are correct it seems likely that 
a large fraction of the anomalous scattering is due to non-Coulomb short 
range forces. Further information on the forces responsible for the 
anomaly might be given by the sign dependence of the scattering. 
Whittemore and Shutt find no difference within rather large statistical 


errors and it is important that this aspect should be examined in greater 
detail, 
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SUMMARY 


The results of measurements on nuclear disintegrations caused by 
50-125 Mev protons in photographic emulsion are given. The relative 
numbers and the angular and energy distributions of the secondary 
alpha-particles and protons have been found and the collision mean free 
path for protons in emulsion deduced. 


§1. INTRODUCTION 

In this paper the results of an investigation of the nuclear disintegrations 
or “ stars ’ produced in photographic emulsion by 50-125 Mev protons are 
presented. They refer to the mixture of elements constituting photo- 
graphic emulsion, and are consequently difficult to interpret theoretically. 
In a subsequent paper an attempt to separate the disintegrations into 
those in the light elements in the emulsion (carbon, nitrogen and oxygen) 
and those in the heavy (silver and bromine) will be made. Apart from a 
few general remarks, all theoretical discussion of the results presented 
here is postponed until that paper. 


§2. EXPERIMENTAL PROCEDURE 


Ilford C2 Nuclear Research emulsion, 200 microns thick, was exposed 
to the external proton beam of the Harwell cyclotron. The beam passed 
through blocks of aluminium of various thicknesses, which reduced its 
energy from an initial value of 147 Mev to 50, 75, 100 and 125 Mev 
respectively for the different exposures. The thicknesses of these blocks 
were calculated using the range-energy curves of Smith (1947) for protons 
in aluminium. The protons entered the emulsion at an angle of about 
three degrees and consequently, on the average, they passed through a 
considerable thickness of emulsion before collision. The average energies 
of the protons at collision were found from the inclinations of their paths 
and the rates of energy loss in emulsion to be 45, 68, 94 and 122 Mev 
respectively. The total spread of the proton energies was about 6 Mev. 

The emulsion was sufficiently well developed for the tracks of the 
incident protons to be visible but the tracks of secondary protons of over 
30 Mev could not be found with certainty as they are generally inclined 


* Communicated by the Author. 
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FS the plane of the emulsion. The relation between grain density and 
rate of energy loss was obtained from the tracks of incident protons and 
from those of stray protons ending in the emulsion. It is shown in fig. 1. 


Fig. 1 


50 
Proton Energy in MeV. 


—— N=N, {\-exp Z/o)} 
(Barkas, 1951) 
No = 2:25 grains per micron 
Jy = 20-4 KeV. per micron 


N=grain density in grains per micron. 


2 Breet ASS 7, 10 
I=rate of energy loss in Kev per micron. 
Grain density as a function of rate of energy loss for C2 emulsion. 


The emulsion was scanned at high magnification with an oil-immersion 
objective, and two hundred events found for each energy of incident 
proton. Events due to collisions between incident protons and hydrogen 
nuclei in the emulsion were excluded. 


§3. THE Sizz DistRIBUTION 


The distributions of star sizes are given in table 1 for the four energies 
of incident protons. 


Table 1 
Proton Number of tracks 
energy 1 2 3 4 5 6 te 
45 38 53 53 49 rl ~_ ~ 
68 46 63 56 23 10 2 = 
94. 15 68 70 30 13 4 - 
122 18 61 51 44 >All 4. il 


The mean size M and the standard deviation W of these distributions 
was calculated for the events with two or more tracks. One-track events 
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were excluded as the difficulty of finding them makes their frequencies 
uncertain. These values of M and W are compared with the results of 
other observers in figs. 2 (a) and (6), and show good agreement. 


Hig. 2 


(a) 


GERMAIN, 1951 


PROTON 
eee Vg LEES er 27, 1955 
PRESENT WORK 


STARS 


sure A TITTERTON, 95! 


+ FISHMAN & PERRY, 
1952 


Jt-MESON MENON, MUIRHEAD 
STARS AND ROCHAT, 1950 


50 100 150 200 
Energy of incident particle in Mev. 


Mean star size as a function of energy of incident particle. 


(2) - 


50 100 150 200 
Energy of incident particle in Mev. 
Width of size distribution as a function of energy of incident particle. 
4D2 
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§4, EneRGy DISTRIBUTION OF THE SECONDARY PARTICLES 


The ranges of all the secondary particles emitted from the nuclear 
disintegrations and ending in the emulsion were measured, and their 


Fig. 3 


20 ONE- AND TWO-TRACK STARS 


(a) 


Corrected number of particles per Mev. 


Oo @ 4 6 8 10 12 14 16 18 20 


Energy in Mev. 
Energy distribution of alpha-particles from one- and two-track stars. 


THREE-TRACK STARS 


(3) 


Corrected number of particles per Mev. 


ce) 2 A 6 8 Te) \2 \4 I6 18 20 


Energy in Mev. 
Energy distribution of alpha-particles from three-track stars. 


energies found from the range-energy relations. The tracks of alpha- 
particles were separated from those of protons by inspection. It was not 
possible to discriminate between protons, deuterons and tritons and so 
all singly-charged particles were grouped together. Similarly, particles 
of charge greater than two were grouped with the alpha-particles. 
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The energy distributions of the alpha-particles are shown in fig. 3 
for stars of different sizes for all energies of incident protons. They are 
grouped in this way rather than as a function of incident proton energy 


Fig. 3 


FOUR-TRACK STARS 


(C) 


fe) 2 ot IPT ES 12 Att ee ee OO 


Corrected number of particles per Mev. 


Energy in Mev. 
Energy distribution of alpha-particles from four-track stars. 


No 
On 


FIVE- SIX- AND SEVEN-TRACK STARS 


Do 
Ze; 


(d) 


On 


Corrected number of particles per Mev. 


O 


a 4 6 8 lO l2 14 I6 13 20 
Energy in Mev. 
Energy distribution of alpha-particles from five- six- and seven-track stars. 


because it was found that the distributions of alpha-particles from stars 
of various sizes were significantly different from each other while those 
from stars due to protons of different energies were not. A peak at about 
10-15 mey, which may be interpreted as due to alpha-particles emitted 


Corrected number of tracks per Mev. 
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from the heavy nuclei silver and bromine, is visible in the distribution 
corresponding to one and two track stars, but not in the others. It is 
likely that nearly all of the alpha-particles of less than about 9 Mev are 
emitted from the light elements of the emulsion. If this is so, then fig. 3 
shows their energy distributions, apart from superimposed distributions 
due to alpha-particles from the heavy elements, which are an appreciable 
fraction of the whole only for the one and two track stars. 

The energy distribution of the secondary protons of less than 10 Mev 
were found to be rather insensitive both to star size and to energy of 
incident protons. The results for all stars were therefore added together 
to give the distribution shown in fig. 4. 


Fig. 4 


fo) 
ie) 


8 


S) 


fo 
ie) 


*O ) 2 3) 4 ) 6 Ys 8 S fe) UI 


Energy in Mev. 
Energy distribution of protons from stars of all sizes. 


The usual corrections for particles passing out of the emulsion were 
made in all cases. 


§5. Tur ANGULAR DiIsTRIBUTION OF THE EMITTED PARTICLES AND 
THE RELATIVE NuMBERS OF ALPHA-PARTICLES AND PROTONS 


A convenient measure of the angular distribution of the emitted 
particles is the forward fraction F, which is defined as the fraction of 
all the emitted particles which go in the forward direction in the 
laboratory system. This was found as a function of star size for all 
particles together, and then for the alpha-particles and protons separately 
and is given in table 2. The fraction P, of all the emitted particles having 
charge two or more is also given in the table. 
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These results are in agreement with those of Seed (1952) who found a 
forward fraction of 0-79 for all emitted particles from stars produced by 
100 Mev protons. 

P,, increases sharply from a value 0-28 for one and two-track stars to 
about 0-55 for three- to seven-track stars. This may be interpreted in 
the light of previous investigations of the disintegrations caused by 
cosmic radiation, in which it was shown that P, for stars in the light 
elements was about 0-50 and in the heavy about 0-25, if the one and two- 
track stars are mainly due to disintegration of the heavy elements while 
all larger stars are mainly due to the light elements. This conclusion is 
similar to that reached by Lees, Morrison, Muirhead and Rosser (1953) 
from an examination of disintegrations caused by 130 Mev protons. 


Table 2 


Star F F fF IE» 
size Alpha-particles Protons All particles 
: 0-77 +0-22 0:79 +0-13 0-79 +0-11 0-28 +0-04 
2 0-72 +0-10 0-62 +0-06 0-65 +0-05 0-28 +0-03 
3 0-69 +0-07 0-59 +0-06 0-63 +0-05 0-52 +0-04 
+ 0-71 +0-06 0-62 +0-07 0-67 +0-05 0-61 +0-05 
6 0-75 +0-09 0-76 +0-10 0-75 +0-06 


The values of F do not, however, show a corresponding discontinuity, 
which indicates that the forward collimation of particles from the heavy 
elements is not markedly different from those from the light elements. 
This conflicts with the concept of evaporation from heavy nuclei which 
are moving slowly forward under the impact of the primary particle for, 
if this were so, the particles from the heavy nuclei would show only a 
slight forward collimation. It is in agreement with the results of 
Bernadini, Booth and Lindenbaum (1952) who concluded, from an 
examination of disintegrations produced by 350 Mev protons, that even 
the low energy particles from excited heavy nuclei are to a great extent 
the result of a nuclear cascade and not of an evaporation process. 

All tracks shorter than four microns, which are mainly due to the 
recoil of the residues of heavy nuclei, and also tracks of protons of over 
30 Mev, were excluded when F and P, were computed as the difficulty 
of finding them makes their frequencies uncertain. These tracks are 
strongly collimated in the forward direction. 


§6. Tur Cottision Mean Free Pato OF PROTONS IN EMULSION 


The relative values of the collision mean free paths for production of 
events of two or more tracks by protons in emulsion were found from the 
numbers of events of two or more tracks per unit volume of emulsion 
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and the relative beam currents in the four exposures. These were 
converted to absolute values by counting the number of tracks of incident 
protons in the 45 Mev exposure. This could not be done for the other 
exposures owing to the difficulty of finding all the tracks of high energy 
protons. ee 

The results are shown in fig. 5 as a function of incident proton energy. 
The values found in the present experiment are rather lower than those 
of Germain. This is possibly because more short tracks due to recoil 


Fig. 5 


150 


ete) 


Mean free path in CM. of emulsion. 
n 
fe) 


b GERMAIN, 1951 


J PRESENT WORK 


O 50 100 ISO 200 250 
Energy of incident protons in Mev. 


Mean free path for protons in emulsion producing events of two or more tracks. 


fragments were observed in the fine grain C2 emulsion used here than in 
the G5 emulsion used by Germain. As a result, more events are counted 
as two-track stars, resulting in a lower mean free path. 
Lees et al. (1953) have pointed out that mean free paths found by ‘ area ’ 
scanning tend to be considerably longer than those found by the more 
accurate ° track’ scanning and they attribute this to the difficulty of 
making properly overlapping scans when using high power objectives with 
very small fields of view. This effect also could account for some of the 
discrepancy between the present results and those of Germain. 
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SUMMARY 


The probability distribution function m(n, t) of the number of registered 
events n in a time interval ¢ for a counter in the case when the registered 
and unregistered events are followed by different dead times, is derived 
by using the method of product densities formulated recently by one of 
us (R.). 


§1. STATEMENT OF THE PROBLEM 


A COUNTER is supposed to register random events such that the probability 
of an event in any small interval of time of length 4t (At is an infinitesimal 
quantity) independently of previous events is «4¢+O(4t)?, or simply 
At (provided we choose the unit of time suitably). The probability that 
exactly n events occur during a time interval of length ¢ is given by the 
Poisson law e~‘(é"/n!). Due to resolving time, the counter is unable to 
register all events. It is customary to treat two cases :— 


Type I.—After each registration the counter is locked for a constant 
time a. An event is registered if and only if no registration has taken 
place during a time a preceding it. 


Type II.—The counter registers an event if and only if no event has 
occurred during the preceding time interval of length, a. Here an 
event occurring at a moment when the counter is locked prolongs the 
inoperative period. In theory the counter can remain locked indefinitely. 

In a more general type of counter we assume the registered and 
unregistered events are followed by two different ‘ dead times ’, a and 6 
respectively. 

The problem is to determine x(n, t), the probability that » events are 
registered in time ¢. While this has been solved in the case of counters 
type I and type II (Feller 1948), the authors have not seen any 
derivation of 7(n, t) in the case of a general type counter, though this 
problem was discussed by Klemens (1948). In this paper we shall 
derive an exact expression for 7(n,t) by using a technique developed 
recently by one of us (R.). We shall first formulate the method in a 
general manner and apply it to the problem of counters. 
a a ee ee re 


* Communicated by the Authors. 
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§2. THE ForMULATION oF THE METHOD. 


Recently one of us (R.) has formulated* a ‘ theory of product densities ’ 
in connection with stochastic problems where we have to deal with a 
stochastic variable representing the number of ‘ particles’ distributed 
in a continuous infinity of states characterized by a parameter. The 
word ‘particle ’ was used merely to facilitate the understanding of the 
problem, but the method developed is quite general and can be used in 
a wide variety of problems involving the statistical distribution of a 
discrete number of ‘ particles’ or ‘ entities’ in a continuous infinity of 
states. 

Consider events occurring along the time axis. We can define product- 
densities of events along the time axis. Let the product density of 
degree n of events be f,,(t, t2,...,¢,). Then according to the theory of 
product densities f,,(t,, t2,..., t,) dt, dt,, ... dt,, represents the probabil- 
ity that an event occurs in the interval between t and ¢,+dt,, one 
event between f, and f,+dt.,..., and one between ¢, and ¢t,,+dt, 
irrespective of the occurrence of events elsewhere. Let us assume we 
deal with a class of product densities defined by 


Frltrs tas +++ t= fila—t) -.--filr—tra) - + - » (I) 


provided 7,,7,,.-. 4%, are. ordered (f,<t,<i,<:...<t,){ where f, is 


the product density of degree one. This represents a class of stochastic 
processes examples of which are easily found. For instance 


1. Poisson events (a trivial example), 
2. Electron counter problems, 
3. ‘ Renewal ’ problems, 


can be cited as illustrative of such processes. 

If z(n, t) represents the probability that events have occurred in 
time t, what happens between ¢ and t-+dt is not only dependent on the 
fact that » events have occurred in time ¢ but upon the ‘ prehistory “— 
i.e., when the last event happened. Hence the usual method of expressing 
a(n, t+-dt) in terms of z(n, t) is unsuitable for such problems. We shall 
show that the product density technique or the method of regeneration 
points yields the required results. 


A. Product Density Approach 


Let G(u, t) be the generating function corresponding to 7(n, t) and by a 
result proved in a previous paper f 


* For a complete formulation of the theory, see Alladi Ramakrishnan, 1950, 
Proc. Camb. Phil. Soc., 46, 595 ; 1952, Ibid., 48, 451. 

+ If t,, t,,..., 4, are not ordered, f, still exists. /, is a product density and 
hence is symmetrical in the variables ¢,, t,,..., t,. We can therefore order 
them and then apply (1). 

¢ Alladi Ramakrishnan, 1952, Proc. Camb. Phil. Soc., 48, 451. 
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n t t t t 
ane at, | di wore ee i T(t taste at agldhege 
ot” 0 0 0 
Since f,, is symmetrical in ¢;, fy,...,¢, we can write the above as 
t t Pee to 
nt | dt, {"atya | di var: .| THRE EA 
0. 0 0 0 
t ty, to d a 
an! | di, {"dta----[ flbdfile—h)---filtetead d+) 
0 0 0 


Let the Laplace Transform* (L.T.) with respect to the variable ¢ of the 
function f,(t) be #(s)= | e-s'f(t) dt. 
0 


t 
Then since [“ 2 = [fi (7) dr, 


the L.T. of [| is (8) 5 a ee re ae 
teed Paes 8 
and of oe | is n! — (4) 
If g(u, s) be the L.T. of G(w, t) 
ikea) ee eee (5) 
3° T—@—1 de) 
If p(n, s) is the L.T. of x(n, t) 
1 [(s)]" 
p(n, = 5 lee top ae On 


In any stochastic problem of this type we must determine either #(s) 


or f,(t). 


B. Approach by the ‘ Method of Regeneration Points ’ 


Let us consider the above class of stochastic processes from the point of 
view of the method of regeneration points. The relation (1) between the 
product densities of degree 1 and n is equivalent to stating that if we are 
given that an event occurred at t=0, the probability that an event occurs 
between ¢ and t+-di (whether or not an event occurred between 0 and ?) is 
determined by the function f,(t) dt and is independent of what happened 
before t=0. It also follows that, given an event happened at t=0, the 
probability that the neat event happens between t and t-+-dt is also inde pen- 
dent of what happened before t=0. The probability can be expressed by 
the function (—d¢(t)/dt)dt where ¢(t) represents the probability that no 
event occurs between 0 and ¢ given that an event occurred at t=0. 


t 
Usually ¢(t) can be expressed as exp | - J a(7) ir | where «(7) dz has the 
0 


obvious meaning of the probability that an event occurs between 7 
and +-+dr given that no event occurred between 0 and 7 and that an event 


* In this paper, L.T. means Laplace Transform. 
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occurred at t=0. According to the method of regeneration points we 
at once write the integral equation 


n(n, in t=— | mt iA eM it PE ees 
where ¢'(7)=0¢(7)/57 and 8(n)=1 for n=0 and zero otherwise. . . (8) 
t 
G(u, t)\=— \, u Gu, t—) $"(7) dr+-¢(0), 
_ g(%, s)=—g(u, s)(so(s)—1) wto(s), . . .'. . . . . . (9) 
where o(s) is the L.T. of 4(#) 
g(u, 8)= 2 (10) 
1+-u(so(s)—1) 
Ag(w, 8) (s)__ 1 1—so(s) 
ee | Ou | -=+ sa(s) (11) 
1 1 
or o(s)=—. 1+4(s) (12) 
; a cle 
ie He = TTY)’ mh 
a result derived by the product density method. 
Thus p(n, s) the L.T. of z(n, t) is 
ay CO) ae n 
a fled (1—so(s)) A opty ears Cia (14) 


In any stochastic problem we must determine o(s) or ¢(s) i.e., d(t) or f,(t) 
in order to obtain p(n, s), i.e., 7(n, t). We shall now oats the above 
results to the problem relating to counters. 


§3. THE SOLUTION OF THE COUNTER PROBLEM 
In the following discussion we assume (i) that the initial condition is 
specified as follows: At t=0 an event is registered. We require the 
probability z(n, t) that n events are registered in time ¢. Omitting the 
initial event n(n, t)=0 at t=0 for nO, 
(0, 0O)=1, 
(ii) that ¢(¢) represents the probability that no event is registered in 
time ft, 
(iii) that f,(¢) represents the product density of degree one of registered 
events. 
For type I counter it is more easy to compute ¢(¢) rather than f(t) 
d(t)=e"*™ for t>a 
d(t)=1 for t<a 
L.T. of f(t) is 
e724 


1 
WO 9 aH 
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L.T. of z(n, t), the probability that n registrations occur in time ¢ is 
a e—nas 1 e—% | 
a(s)(1—o(s) ° S) = (s+1)" E = s(s-+1) 


1 e-nas e(nt+las ] 
= Lene a} 


n—-1 . 
Thus cae t) =H(t—na) et" —_ —H(t—(n+1)a) e-(-nt g 
t 1h 
t= (m+ 1a)" 
n! 
t—na e—tn—-1 t—(n+1)a e—tyn 
— nt oe See 
each. or 08 


where H is the Heaviside unit function, i.e., 
Ale)=it 20.40 (2) — Oe 
For the counter of type II it is easy to calculate f,(¢) the product density 
of the first degree relating to Paap ik events. 
fii) di=H(t—a)e* dt. = ee 


where H is the Heaviside unit aah The above expression results 
from arguing that an event is recorded when it is not preceded by an event 
in time interval a and is not also preceded by the event occurring at t=0 
within a time intervala. ‘Thus 


1 
ab(s) the L.T. of fi ()= 3s eustD * (18} 
And this yields 
—1)e~* (t—a u—1)¥ eX (t—Na)¥ 
gypsies ee eae ae ret) 
where WN is anumber such that (V+1) a>t>Na 
m=N —ma 
a(n, t= x (—1)"-" (") “— ] (t—ma)" <> a Fae as 


a result derived in a previous paper.* 

We can now derive the L.T. of x(n, t) when it represents the probability 
of n registrations occurring in time ¢ for a more general type of counter in 
which the registered and unregistered events are followed by two different 
‘ dead times ’, a and b respectively. We proceed as follows : 

At t=0an Seen is registered ; we ask : what is the probability ¢(¢) that 
no event is registered up to time ¢t? Let us assume a>b. We can write 
P(t) as 

(t)= {H(t)—H(t—(a—b))}+-¢"(t—(a—b)) H(t—(a—b)). . (21) 
The first ae expresses the condition that no event is registered for 
t<a—b. ¢(t)=1 for t<a—b. 
i 
* Alladi Ramakrishnan, 1951, J. Roy. Stat. Soc. B, 18, 131. 
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¢'\(t—(a—b)) represents the probability that no event is registered in 
time ({—(a—6)) on the hypothetical assumption that we are dealing with a 
counter of type II with dead time b. As has been previously mentioned, 
for type II counter we do not usually compute ¢"(t) but only f(t), the 
product density of registered events of a type II counter with dead time b. 
But the L.T.’s of $"(¢) and f,"(¢) are connected by the relation (12). Thus 
L.T. of d()=L.T. of {H(t)—A(t—(a—b))} 

Fle eoted (¢——(@--0)) teem . wake 4) (22) 

We know that the L.T. of f,"(é) is 


1 
5 ev(SF1) * 


Thus we have L.T. of f,(t) the product density of registered events of the 
generalized type of counter given by 


¥(8)= — Saath pg abe] * OH copie. oe) 


Put b=0 we get counter type I. 

Put a=b we get counter type II. 

Put a=b=0 we get a Poisson distribution for 7(n, t) characteristic of an 
ideal counter. 

The same result is obtained for b>a following exactly the same 
arguments. 

We have, for the L.T. of z(n,¢) for general a and 6, the expression 


er = = (s ecto e(a—b)s)—(nt1) (s e%stb 1 e(4—b)s__ 1), : (24) 


Taking the inverse 


o+io ] 
Aas t) 1 | E (s etstb_t e(4—b)s)—n __ = (s osrt Pp glen tnta | eds. 


~~ Qari o—io0 
(25) 
e O(s+1) —n 
1 (s certty-t 1 : } 


foo} ad —rb(s+1) 
aie (See eee li (x ') ae , 
s r 


1 
Since —(s e%t!+-el4—6)s)-n— 
s 


% | 


- 
r=0 s 


es (oe mtr—l\ estne-rd) 
a(n, t)= | a ds jetn(—1y ( ae ) ee (ert 1) 


Qt a—ino r=0 r Cut ay 
n-+-r es(t—(n+1)a—rb) 
x(1y( re Pee MRS Webelos isin ote. oa, (20) 
Remembering that 
G+10 @sz n—1 
al ; ee en ee eee eS z>0, and zero forz<0, . . (27) 
271 o—jo Be (n—1)! 
we get 
Na n+r—l)\ (t—na—rb)"*" 
a(n, t)= ¥ (1 ( y )\ ae b(wr) 
f= 
x —l)" ee (Ge Naa b)T ET j-bintren es (28) 
sae r (ntr+1)! 
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where N,, N, are integers such that 


a ne pe 1 rte 
i 7 -1<N,< and ae ee 


Particular cases of the above solution for b=0, and for a=b agree with 
those already given. 

When a0, 40, we note that z(n, t) is given by a finite series ; also 
a(n, t)=0 if na>t. Hence the moments of n are calculated from the 


definition 


N 
w= LX nnn, t)= 2X n'n(n,t) where N 


r=0 r=0 
is an integer such that (V+1)a>t>WNa. oN eet gt NEG te md 
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CXVIII. The Onset of Convection by Thermal I nstability in Spherical 
Shells—(A Correction) 


By 8S. CHANDRASEKHAR, F.R.S. 
Yerkes Observatory, University of Chicago* 


[Received July 6, 1953] 


In a paper (Chandrasekhar 1953) recently published in this Magazine, the 
author has considered the criterion for the onset of convection by thermal 
instability in spherical shells and has described a variational procedure 
for solving the underlying characteristic value problem. Unfortunately 
in the numerical evaluation of the final results an error has occurred due 
to the oversight that @’,,4/o(«;) and ©’, 1/2(«;y) which occur in equations 
(21), (22) and the subsequent ones are not the same as those defined in 
equations (14) and (15): actually in equations (21) and (22) (and in all 
later ones) @';,1/9(«;) and @’,,4/0(%;n) are 1/a, times the values given in 
equations (14) and (15). 


Table 1. The Characteristic Numbers C, for Various Values of / and 7 
l 4 n=0-4 9=0°5 7=—0-6 7=0°8 
J 4-999 x 10? | 1-682 x 104 | 4-188 x 104 | 1-403 x 10° | 7-789 x 106 
2 EO007. 10> sel-09 1s 10S 022181 >< 104) 6-133>¢10* | 2-753 x 10° 
3 8-881 x 10? | 1-196 x 104 | 1-924 104 | 4-424 x 104 | 1-500 x 10° 
+ 1-400 x 104 | 1-585 x 104 | 2-146 x 104 | 4-076 x 104 | 1-005 x 10° 
5 2-121 x 104 | 2-227 x 104 | 2-673 x 104 | 4-313 x 104 | 7-656 x 10° 
6 3°143 x 104 | 3-492 x 104 | 4-945 x 104 | 6-368 x 10° 
Dy 4-365 x 104 | 4-629 x 104 | 5-933 x 104 | 5-651 x 10° 
8 6-125 x 104 | 7-292 x 104 | 5-270 x 10° 
) 8-027 x 104 | 9-057 x 104 | 5-109 x 10° 

10 1-039 x 10° | 1-128 x 10° | 5-104. 10° 
ll 1-325 x 10° | 1-401 x 10° | 5-223 x.10° 
12 1-669 x 10° | 1-732 x 10° | 5-448 x 10° 
13 2074 x 10° | 2°126 x 10° | 5-767 x 10° 
14 2-549 x 10° | 2-590 x 105 | 6-178 x 105 
15 3:099 x 10° | 3-131 x 10° | 6-678 x 10° 


The corrected values of C, are given in table 1; this table replaces 
table 1 of the paper to which reference has been made. The dependence of 
C, on 1 and 7 is illustrated in fig. 1 ; it will be noticed that this figure does 
not have the unsatisfactory feature of the figure based on the erroneous 
calculations : the curves for the different ’s do not cross but lie systema- 
tically one above the other. 

The principal conclusions of the paper are not affected by this correction. 


* Communicated by the Author. 
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Fig. 1 
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The criterion (C7) for the onset of convection in the mantle as a spherical 
harmonic disturbance of order / for various thicknesses of the mantle. 
The different curves are distinguished by the values of the fraction (7) 
of the radius of the sphere occupied by the inviscid core. 


REFERENCE 
CHANDRASEKHAR, 8., 1953; Phil. Mag. [7], 44, 233. 


eelioten 


CXIX. Spallation of Vanadium, M anganese, and Cobalt with 
187 mev Protons 


By 8. G. Rupstam 


The Gustaf Werner Institute for Nuclear Chemistry, University of Uppsala, 
Uppsala, Sweden* 


[Received June 17, 1953] 


SUMMARY 


This paper describes measurements of the yield of various spallation 
products from vanadium, manganese, and cobalt. From the experiments 
it is concluded that in the spallation of medium weight elements with 
high energy protons, the total cross section per mass number decreases 
almost exponentially with increase in the difference between the actual 
mass number and the mass number of the target nucleus. It is also 
concluded that for a given mass number the cross sections of the isobars 
are approximately represented by a Gaussian distribution. 


§1. INTRODUCTION 


A REACTION in which a high energy particle hits a nucleus and knocks 
out a number of low mass particles from the nucleus is called a spallation 
reaction. If the residual nucleus is radioactive, it is often possible to 
measure its yield and so determine the cross section for the spallation 
reaction, which gives this residual nucleus. The purpose of the present 
investigation was to determine the cross section distribution, i.e. how the 
total reaction cross section was distributed among the product nuclides. 

Previous spallation investigations of medium weight elements have 
been carried out for arsenic with 180 Mev deuterons (Hopkins and 
Cunningham 1948), copper with 340 Mev protons, 180 Mev deuterons, and 
380 Mev alpha particles (Batzel, Miller and Seaborg 1951, Miller, Thompson 
and Cunningham 1948), copper with 400 Mev neutrons (Marquez 1952), 
cobalt with 240 mev protons (Wagner and Wiig 1952), iron with 340 Mev 
protons (Rudstam, Stevenson and Folger 1952). 

Ideally a spallation investigation should be carried out with a single 
isotope target. For this reason vanadium, manganese and cobalt were 
chosen as target materials for the present experiments, although the 
vanadium target contained 0-25% of °°V. 


* Communicated by T. Svedberg. 
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§2. EXPERIMENTAL 


(a) Procedure 


The target materials were analysed chemically (V, Mn, Co) and 


spectroscopically (V). 


to interfere with the measurements. 


Table 1 


The analyses did not show any impurities likely 


i 


Element Main purification method Weighed as 
Cr Extraction of perchromic acid with diethyl ether. | BaCrO, 
V Precipitation with cupferron. Extraction with | V.0; 
cupferron and chloroform. 

sie Precipitations with potassium iodate and with | TiO, 
cupferron. 

Se Extraction from 50% potassium thiocyanate | Sc,O3 
solution with diethyl ether. 

Ca Precipitation with oxalic acid.. Destruction of | CaC,0,.H,O 
the oxalate and hydroxide-scavenge. 

K Hydroxide-, sulphide-, and oxalate-scavenges. | KCIO, 
Precipitation with perchloric acid. 

Cl Distillation of hydrogen chloride. Precipitation | AgCl 
of silver chloride. 

S Precipitation of barium sulphate, which was dis- | BaSO, 
solved in a saturated potassium carbonate 
solution. 

RP Precipitations with molybdic acid and with | Mg,P,0, 
magnesia mixture. 

Si Carrier: ammonium fuosilicate. Precipitation | SiO, 
of silicic acid with boric acid. Dissolution in 
alkali. 

Mg Precipitation with diammonium ortho-phos- | Mg,P,0, 
phate. 

Na Hydroxide-, sulphide-, and oxalate-scavenges. | NaCl 


Precipitation with zink uranyl acetate. Pre- 
cipitation of sodium chloride with alcohol 
saturated with hydrogen chloride. 


—>- | .<— | \, = = = 7 kik, =. 5 laclanjg ila iE DninGn cr a 


The targets were bombarded with 187 Mev protons in the ‘circulating 
beam of the Uppsala synchrocyclotron. After irradiation the target was 
dissolved in nitric acid containing carriers for the elements to be studied 
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in the experiment. The carriers were then separated and purified very 
carefully, and finally precipitated in a form suitable for gravimetric 
analysis. The precipitate was spread as uniformly as possible on a 
weighed aluminium or platinum dish (diameter 13 mm) and then dried or 
ignited. The dish was weighed again and mounted on a small aluminium 
holder for particle counting. 

An end-window Geiger—Miiller counter in an ordinary lead house 
set-up was used for counting the activity present. 

The separation and purification methods used are indicated in table 1. 


(b) Bombardments 
A summary of the bombardments is given in table 2. 


Table 2 
Bombard- Duration of 
ment bombardment Elements isolated 
No. (min) 
Vi il Cries 
We 80 Cl 
Ve3 60 Cri, Ca, K 
V 4 61 We, Gi S.1Ps ts 
V5 52 K,S 
V6 20 V, Ti, Se, Si 
AW 3 Cre Vemscw barNa: 
V8 9 V, Ti, Mg 
V9 5 V 
V 10 60 We Cig GIL IS, 1s IN) 
AW ahi 10 V 
Val2 10 VEC 
Mn 1 143 CR OL 1 
Mn 2 61 AR, UC 
Mn 3 65 Lae SC @ae hers 
Mn 4 76 Grekysi 
Co 1 60 (Cie WV I S05 18S) 12, Si 


In most of the vanadium bombardments the nuclide 4°V was used as 
a monitor. The absolute cross section of this monitor was determined 
by comparing the yield of it with the yield of *4Na from an aluminium 
foil covering the same area as the target. The cross section for the 
reaction 27Al (p, 3pn)24Na was taken to be 10-8 millibarn (Stevenson and 
Folger). The nuclide “Ti was chosen as a monitor in the manganese 
bombardments. 

§3. RESULTS 

The cross sections were computed from the measured data in the 
usual way. Corrections had to be applied to the counting rates for 
coincidence losses, background activity, geometrical differences, absorp- 
tion in the window and in the air between the sample and the window, 
backscattering, self-absorption and self-scattering, counting efficiency, 
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chemical yield, and decay (Burtt 1949). The correction factors for the 
backscattering were obtained from the papers by Burtt (1949) and by 
Yaffe and Justus (1949). A smaller correction factor was used for 
positrons than for electrons (Seliger 1952). The correction factors for 
self-absorption and self-scattering were estimated from the work of 
Nervik and Stevenson (1952). 

Owing to the great number of corrections which had to be applied, 
very accurate results were not to be expected. The measured cross 
sections were estimated to be accurate to 50%. 

The cross sections for the spallation of vanadium, manganese, and 
cobalt are listed in the tables 3, 4 and 5 respectively. The cross sections 
in table 3 and table 4 are, with few exceptions weighted means of 2-4 
determinations. 


Table 3. Spallation of Vanadium 


Cross Cross 
Nuclide section Nuclide section 
(mb) (mb) 
49Cr 1:8 431K 2-5 
48Cr C:23* 42 1c 3°9 
18 28 29@) 0-30 
47VV 5:7 38C] 0-74 
dBi 6-1 24C] 0-03] 
48Se 4:8 358 0-61 
47Se 10 iB 0-32 
416S¢ 20 32P. 0:34 
45cm 8-4 318i 0-18 
4388c+-“48c 13 -™Mg 0-028 
47Ca, 0-11 =4Na 0-030 
45Ca 3°6 22Na 0-077 


* Computed from the amount of 48V formed. 
+ Uncertain since the activity of ?2Na was low. 


It is difficult to separate the chlorine isotopes by particle counting. 
since their half-lives are nearly the same. The cross sections of the 
chlorine isotopes given in table 3 have been kindly submitted to me by 
Mr. EK. Arbman of this institute, who used a betaspectrometer for the 
determinations. . 


§4. Discussion 


Devoting the first part of the discussion to the determination of the 
cross section distribution for vanadium, it is clear that this can be 
described if two terms are known, namely (a) the total cross section per 
mass number (=0,(M)) and (b) the cross section as a function of the 
atomic number (—o(M, Z)) for a given mass number. 


o(M)=Zo(M, Z). 
Z 
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Tn order to determine o,(M/) it is necessary to make certain approxi- 
mations. Every mass number in the actual region of the nuclide chart 
contains at least one stable isobar, for which the cross section has to be 
estimated. To be able to carry out these approximations it is assumed 
that the cross sections are independent of the individual properties of the 
product nuclides. This is reasonable when the energy of the bombarding 
particles is high. Then it is to be expected that the cross sections of 


Table 4. Spallation of Manganese 


Cross section Cross section 
Nuclide relative to Nuclide relative to 

4577] 47S. 
Or 2-2 48S¢ 0:6 
eo ik 4°Sc¢ I 
=i) 0-007 46Sc 2-0 
25Ca 0:25 44Scm (kay 
431 0-21 43S c¢+448¢ 1:3 
eK 0-39 
32@ 0-008 

340] 2380] 0-013 
32 0-034 
318i 0-027 
Table 5. Spallation of Cobalt 
Cross section Cross section 
Nuclide relative to Nuclide . relative to 

49Cy 47Se6 
49Cr 1 48Sc 0:35 
487 The 47S¢ 1 
45T 7] 0:45 46Sc 4-1 
42K 0-16 44Scm 1-7 
32P 0-025 43S¢ -+44S¢ ley 
318i 0-010 


nuclides belonging to the same isodiaphere should vary monotonically, 
i.e. decrease monotonically with increasing distance from the target 
nucleus. That this is really the case will be evident from fig. 1, which is 
a plot of the cross sections of all the investigated nuclides as a function 
of the mass number. Curves are drawn through nuclides belonging to 
the same isodiaphere and the isotopic number J is noted. It is clear 
from the figure that the cross sections vary in a quite regular way. It 
should be possible to obtain a fairly good approximation of the cross 
section for any nuclide by interpolating or extrapolating in the curve 
with the same value of J as the nuclide in question. 
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The total cross section o,(M) of a certain mass number can now be 
estimated in the following way. Suppose the mass number to be studied 
is 34. The cross section of 24Cl (and of all the neutron deficient isobars) 
has been measured. The cross section of 34S and of all the neutron 
efficient isobars can be extrapolated from the curve with /=2 in fig. 1. 
The sum of these cross sections should give o,(34). In this way the 


Fig. 1 
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curve in fig. 2 has been constructed. This figure gives o,(M) in a 
logarithmic scale as a function of M. The decrease of the cross section 
with decreasing mass number is seen to be nearly exponential between 
M=48 and M=32. 

It now remains to determine the variation of the cross section with 
the atomic number for a given mass number. Consider the mass number 
47. The cross sections of 47V, 47Sc, and 47Ca have been measured. 
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Figure 2 gives o,(47) so that the cross section of 47Ti can be obtained by 
subtraction. The same procedure can be undertaken for the mass 
numbers 43, 45, and 48. The variation of the logarithm of the cross 
section as a function of the atomic number for these four mass numbers 
is shown in fig. 3. 

It is true that some of the measured cross sections include the cross 
sections of parent nuclides. However, it is evident from fig. 3 that the 
daughter nucleus generally has a much larger cross section than its parent, 
the cross section curves of fig. 3 having their maxima close to the stability 
line. For that reason the measured cross section of a nuclide does not 
differ much from the primary yield of that nuclide. 


Fig. 3 
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The experimental points in fig. 3 can be fitted rather well to parabolas 
as is shown in the figure. These parabolas can be written : 


log o(M, Z)—log o(M, Z,)=—p(Z—Z,)?> ~. . . - (1) 
(p, o(M, Z,), and Z, are constants). 
This means that the distribution of the cross sections of isobars with 
the same mass number is approximately a Gaussian one : 
o(M, Z)=o(M, Z,) exp [—(p/log e)(Z—Z))?]. wae (2) 
All the parabolas in fig. 3 have been constructed with the same value 
of p, namely p=0-78, which gave a fairly good fit to the experimental 
points. The Gaussian distribution (2) has a variance of 0-53 atomic 
number units, if the value 0-78 is chosen for p. 
In the following estimations of cross sections the distribution (2) 
with p equal to 0-78 was arbitrarily extended down to the mass number 22. 
Now, it would be possible to compute the cross section of a given 
nuclide from the eqn. (2) if o(M, Z) and Z)(M) were known as functions 
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of M. From the parabolas in fig. 3 the following values of Z,(M/) and 
a(M, Z,) can be obtained (table 6). 


Table 6 
M rap o(M, Z,) (mb) 
43 20-1 13-5 
45 21-1 26 
47 21:9 47 
48 22:3 60 


Z,(M) is plotted as a function of M in fig. 4. The points lie close to 
a straight line which passes through the origin. From this figure Z,(M) 
can be estimated for any mass number. 


Fig. 4 
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The distribution (2), integrated from the smallest permitted vaiue of 
Z to the largest, would give a rough estimate of o,(M). The limits of 
the integration can be extended from — co to + oo without introducing 
serious errors since small and large values of Z give negligible contributions 
to the integral 


a(M)=Zo(M, Z)~ | o(M, Z) exp [—(pllog €)(Z—Zy)?] az 


=o(M, Z,) I, 
(I is equal to 1-32 for p=0-78). 

Since I is independent of the mass number, o(M, Z,) is approximately 
proportional to o(M ). A plot of o(M, Z)) as a function of o,(M) should 
then give a straight line through the origin with a slope equal to 1/1-32. 
Figure 5 shows that this is so. 

From the figs. 2 and 5 o(M, Z,) can now be computed as a function 
of M (fig. 6). 


The eqn. (2) and figs. 4 and 6 describe completely the cross section 
distribution for spallation of vanadium. 
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The measured cross sections are compared in table 7 with the cross 
sections calculated from eqn. (2) and the figs. 4 and 6. In the fourth 
column of the table the difference between computed and measured cross 
sections is given as a percentage of the latter. 

The agreement is fair. This is of course to be expected since many of 
the measured cross sections of table 7 were used for the construction of 
the figs. 4 and 6 and for the determination of p in eqn. (2).* 


Fig. 5 
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The methods of calculation indicated can now be used to estimate the 
cross sections of the spallation products of manganese and cobalt. 
Assuming that the cross section distribution is the same as for vanadium, 
eqn. (2) is again used. The straight line in fig. 4 passes through the 
points (0; 0) and (23-8; 51). (Target nucleus : 33V.) The corresponding 
lines for manganese and cobalt are now assumed to pass through the 


EE a SL haa Sa 
* Note added in proof.—The cross sections of the spallation products of 
vanadium can be described by the formula : 
o(M, Z)=exp [0:358M —12-82—1-80(Z—0-466.M)?] mb. 
This formula can be used down to 27Mg and gives cross sections not very 
different from the calculated cross sections in table 7. 
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Nuclide 


Table 7 


Measured 
cross section 


Calculated 
cross section 


48(1y 
487 
47VV 
45Tj 
48S. ¢ 
47S 
16S. ¢ 
448¢ | 449m 
43S.¢ 
47g 
45a 
43K 
ak 
390] 
38] 
34Q)] 
358 
33p 
s2p 
319} 
27Mg 
24N a, 
2Na, 


Difference 


(%) 


+13 
—2] 

0 
—30 


Table 8 
; Measured Calculated Diff 
Nuclide cross section | cross section : ee 
o/o Ti o/o ©Ti ( %o) 
ar 2-2 3-0 +36 
“Ca 0-007 0-006 —14 
“Ca 0:25 0-51 +108 
SK 0-21 0-18 14 
me 0-39 0:63 +62 
"P 0-034. 0-074 +117 
3181 0:027 0-037 +37 
a/o Sc Ae 47Se 

uSe 0-6 0-2 —67 
16S¢ 2.() 2-6 +30 

44S e +4485 em 1:8 

“Sc 0-5 
43S e +448c +44Se¢m 9.5 2.3 8 
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points (0; 0) and (25-8 ; 55) and (0; 0) and (27:8: 59) respectively. The 
curve in fig. 6 is used except that the abscissa is increased by four mass 
units for manganese and eight mass units for cobalt. 

The measured and computed cross sections for the ‘spallation of 
manganese and cobalt are shown in tables 8 and 9 respectively. 

The agreement between measured and calculated cross sections is 
quite good. From this it can be concluded that the cross section distri- 
bution for spallation of vanadium, manganese, and cobalt with high 
energy protons is rather well characterized by an exponential decrease of 
the total yield per mass number with increasing distance from the mass 
number of the target nucleus and by a Gaussian distribution for a given 
mass number. 


Table 9 
Measured Calculated Dithorencé 
Nuclide cross section | cross section (4) 
alo Cr o/o Cr o 
SAN 1-2 2-2 +83 
baat 0-45 0-32 —29 
2K 0-16 0-10 —38 
32P 0-025 0-043 ete 
318i 0-010 0-022 +120 
alo **Se ojo *"*Se 
aoc 0-35 0-15 —57 
1686 4-] 3-5 —17 
44S ¢ + 448em 3°5 
430 Ihe 
4336 1449 ¢ | 445¢m 3-4 4-6 +35 
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CXX. Electron Microscopic Observations of 111 Faces of ZnS Crystals 
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G. VANDERMEERSSCHE 


Lab. voor Kristalkunde, Rozier, 6, Centrum voor Electronenmicroscopie, 
Plateaustraat, 22, Gent, Belgium* 


[Received July 6, 1953] 


THE object of this communication is to draw attention on two features, 
one concerning growth fronts, the other slip lines, which were observed 
by electronmicroscopic examination of replicas of two (111) faces of a 
large and macroscopically nearly perfect ZnS crystal. 

A previous optical and interferometric study (Votava, Amelinckx 
and Dekeyser 1953), the results of which will be published elsewhere, 
had revealed that the topography of these faces was very complex. 
They present, next to growth fronts and slip lines, areas exhibiting a 
complex substructure, delimiting closed or nearly closed areas of the 
order of 0-001 mm2, which themselves are further subdivided into similar, 
but much smaller regions. 

The observation of these details was made possible, or at least greatly 
increased, by the slight natural corrosion of the crystal, which proved 
also to bring out clearly details on the electronmicrographs. 

It is assumed that the crystal, originating from Jugoslavia, has 
grown at rather high temperature (hydrothermally or from the vapour 
phase) and that rather great internal stresses developed owing to tempera- 
ture fluctuations and impurity content during the final stages of growth 
and after growth had ceased. The replicas were methylmethacrylate— 
SiO obtained in the usual way, and then shadowed with palladium. 
The photographs are negatives, black is shadow. 

Figure 1 (Pl. 46, figs. 1 and 2) represents lines which are interpreted 
as growth fronts. This assumption is based on two arguments: (i) 
their main aspect is identical with lines on other electron-micrographs of 
the same face (not reproduced here), and where we are dealing unambigu- 
ously with growth fronts, as could be deduced from interference with 
obstacles and dislocations ; (ii) they present features, which will be 
commented upon, and which have also been observed, on a different 
scale, on spiral growth fronts on (0001) faces of some carborundum 
crystals (Amelinckx 1952). 

These features consist of lines such as AB which are confined between 
two (rarely more) successive growth fronts, making a nearly constant 
angle with them. The lines between consecutive pairs of fronts form a 


* Communicated by the Authors. 


Electron Microscopic Observations of 111 Faces of ZnS Crystals 1143 


geometrical pattern. In idealized form, this pattern looks like an array 
of quadrangles on the sides of which limited transverse lines are drawn. 

As the segments AB are situated between two growth fronts, it would 
seem likely to consider them as misfit lines separating, in a same front, 
parts which have a different composition (polytypes). We also can accept 
that such a discontinuity can be ‘washed’ away by a following front. 
But we cannot explain by this hypothesis the pattern following which these 
lines are aranged as a whole. 

On the other side, if we consider the lines AB connected with features 
in the interior of the crystal, one can difficultly see why they should be 
confined between growth fronts. So, we find it better to delay any 
interpretation, and to stress that the similar observations on carborundum 
indicate that we are dealing here with something which is not accidental, 
at least for substances with this type of lattice. 

Figure 2 represents sets of lines which are different from the growth 
fronts described so far. Figure 3 (Pl. 47, figs. 3 and 4) is an enlargement of 
part of it. They have been interpreted as slip lines because: (i) the 
corrosion along them occurs in a special way ; (ii) their spacing is irregular 
and their height variable (as deduced from the shadow). 

Spots like A are hills on the replica, they correspond to etchpits. in the 
crystal. As can readily be deduced, the corrosion was confined along the 
slp lines; etch channels, most of which are more or less resolved in 
isolated pits, are seen to alternate with areas where nothing happened. 
This discontinuity can be very well seen along lines as BC. 

If we accept that etchpits develop most easily at points where screw 
dislocations emerge (Gevers, Amelinckx and Dekeyser 1952, Horn 1952, 
Gevers 1952, 1953), we deduce that the lines are defined by the emergence 
of groups of such dislocations, resulting from the action, in the same 
(or nearly the same) plane of successive Frank and Read generators. On 
some occasion they act in a different set of (111) planes (DE) (cross slip). 
This picture corresponds rather well with the ideas put forward by 
Mott (1952). 

The spacing of the lines is rather small, 4. as a mean. Some are closer 
together, and for these something like a fine structure may be present 
(GF). This points towards deformation at temperatures well above 
‘room temperature ’. 

Figure 4 originates from a different (111) face, things are more compli- 
cated here: slip has taken place on two different sets of planes. The 
discontinuity of the shadow shows here very strikingly the rows and 
grouping of etchpits. Some segments are more or less bent. The 
nature of intersection is that one line seems to stop as the transverse 
one is approached, and to start again on the other side (A, A’) sometimes 
with a slight deviation. 

The concentration of dislocations is here extremely great, the opposite 
situation being the case of a well annealed single Al-crystal considered 


by Amelinckx (1953). 
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CXXI. Hxperiments on the Deposition of Iodine 131 Vapour 
onto Surfaces from an Airstream 


By A. C. CHAMBERLAIN 
Atomic Energy Research Establishment, Harwell* 


[Received May 7, 1953] 


ABSTRACT 


Experiments are described on the deposition of ‘ carrier-free ’ Iodine 131 
vapour from a moving airstream onto surfaces. With copper and filter 
paper surfaces and at low or moderate Reynolds numbers, the deposition 
is controlled by aerodynamical factors, and the analogy with heat transfer 
by forced convection holds. The detailed deposition pattern produced by 
an obstacle on a flat plate is shown by a radioautograph. 


$1. InTRODUCTION 


THE purpose of these experiments was to show that the deposition of a 
carrier-free radioactive vapour of iodine from air streams onto certain 
surfaces is controlled by diffusion across the laminar layer of air flow 
over the surfaces, and is governed by the same laws as the diffusion of 
heat in forced convection. 

Pasquill (1943), with plane surfaces of saturated paper, and Ranz and 
Marshall (1952) with droplets of water and other liquids, have shown that 
evaporation from saturated surfaces can be predicted by diffusion theory, 
and the same should be true of deposition of a vapour on a surface, 
provided that the physical or chemical conditions at the surface are such 
that the vapour pressure just above it is zero, so that the surface acts 
as a,‘ perfect sink ’. 

The use of ‘ carrier-free’ radioactive isotopes makes this condition 
easily realizeable, since it is possible to work with less than a mono- 
molecular layer of activity on the surface. Clean copper, and filter 
paper surfaces have been found to act as perfect sinks for the deposition 
of ‘ carrier-free ’ #1I vapour, at the activity levels used in these experi- 
ments. 

The use of radioactive tracers in this way has a possible practical 
application, since by making a radioautograph of the deposited activity 
very fine detail can be seen of the variation of deposition with the 
characteristics of the boundary layer. 


* Communicated by the Author. 
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§2. MeTHOD 
2.1. Source of #41 Vapour 


Iodine 131, prepared in the Harwell pile by neutron bombardment of 
tellurium, was extracted as a solution of elemental iodine in carbon 
tetrachloride. No carrier iodine was added, but iodine impurities in 
the tellurium resulted in there being about 5 micrograms stable iodine 
per millicurie of }"I. 

The iodine solution was dried by standing over concentrated sulphuric 
acid, and sprayed up as a mist of very fine droplets by means of a 
Collison spray operated by dried and filtered compressed air. 

This method of generating !I vapour was not very efficient, as there 
was a considerable loss of 1#1I by adsorption in the spray, but none better 
could be found. Since all the measurements consisted of comparisons 
of the activity deposited on some surface with the activity present in the 
airstream, the exact amount of !I vapour generated was not required 
to be known. 


2.2. Estimation of 11 in Airstream and on Surfaces 


The quantity to be measured was the time integral of the activity per 
cc. of air, or dosage as it will be called, and it was estimated by drawing 
the vapour at a known rate through bubblers containing 20 ml. of 5% 
NaOH, to which a few milligrams of KI had been added as carrier. 
Two bubblers in series were used at each sampling point, and the ratio 
of activity found in the bubbler first in series to that second in series 
generally lay between 20: 1 and 50: 1. 

The deposited activity was measured either by direct counting of a 
cut out sample under an end-window counter, with allowance for back- 
scatter, or by soaking the sample in NaOH solution with KI carrier and 
counting in a liquid counter. 


2.3. Miniature Wind Tunnel 


The wind tunnel used was of square section, 4in. by 4in. The very 
small size had the disadvantage that the area available for working was 
small if the region near the walls of retarded air flow was to be avoided. 

The vapour from a Collison spray, and a much larger volume of air, 
was drawn into the mouth of the tunnel and passed firstly through an 
orifice of 2in. diameter to promote mixing, and secondly through a 
honeycomb of } in. mesh to remove any swirl. A relatively long settling 
length enabled a uniform distribution of wind velocity and activity con- 
centration to be attained before the working section, in which was placed 
the trial surface, was reached. The inlet tubes of the sampling bubblers 
projected through the floor of the working section, four sampling positions 
usually being used. 

Pairs of sampling positions 5 em apart in the working section usually 
gave results agreeing to within 10%. There was a slight fall off in 
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dosage as the sampling position was moved down the tunnel due to 
adsorption of the 1I vapour on the walls of the tunnel, but this amounted 
to no more than about 0-3%, per cm in the working section. 

The wind velocity was measured by a vane anemometer at the front 
end of the tunnel, calibration being effected by means of a pitot static 
tube connected to a Chattock—Fry gauge. 


§3. Resutts or WiInD TUNNEL EXPERIMENTS 
3.1. Definition of Velocity of Deposition 
The velocity of deposition, v,, of a vapour to a surface is defined by 


Rate of deposition per cm? per sec : 
7 Volumetric concentration per cm? (1) 
Amount of activity deposited per cm2 
- Time integral of volumetric concentration, or dosage per cm?. 


3.2. Deposition on Flat Filter Paper Surfaces 

In table 1 are shown the values of v, obtained with a flat plate covered 
with filter paper. The experiments are marked H or V according as the 
plate was placed horizontally or vertically in the tunnel. 

The deposition on the upper and lower surfaces of the plate when it 
was placed horizontally was found to be equal to the accuracy obtained 
in the experiment. Activity deposited on the filter paper was quite 
firmly fixed, the half-life of samples exposed to the air being only slightly 
less than the normal 1*"J half-life of 8-0 days. 

The reproducibility between different experiments at the same wind 
speed was fair, the average difference between corresponding entries 
being 13%, and there was no apparent dependence of the velocity of 
deposition on the dosage of activity passed. 

In all cases there was a fall-off in the velocity of deposition down-wind 
from the leading edge of the plate. At a given distance from the leading 
edge, the velocity of deposition increased with increasing wind velocity, 
approximately as the square root of the wind velocity. 


3.3. Effect of Laminar Boundary Layer 

The falling off in velocity of deposition down-wind from the leading 
edge of a flat surface is due to the building up of a laminar layer of flow 
over the surface. In the turbulent stream of air in the open tunnel, eddy 
diffusion causes the concentration of vapour to be effectively the same 
everywhere. When an absorbing surface is present the air close to it is 
denuded of activity, and the rate at which activity is brought to the 
surface depends on the characteristics of the air flow over it. Where the 
laminar layer is big, there is a considerable distance, of the order of a 
millimetre in these experiments, over which the transporting agency is 
molecular, not eddy diffusion. 

4F2 
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In the radioautograph of Plate 48 the effect is shown of breaking up of a 
well developed laminar layer over a flat surface by an obstacle, a screw 
inserted perpendicular to the surface. In the wake of the screw, a heavier 
deposition of activity has taken place, shown by the dark trail on the 
radioautograph, due to the turbulence created by the screw. 


3.4. Deposition and Boundary Layer Theory 
The results of table 1 can be expressed in terms of a dimensionless 
number N,,, analogous to the Nusselt number of heat transfer. 


Iz 


where J is the activity deposited per cm? at distance 2 from the leading 
edge of the plate, ¢ the dosage of activity in the free air, and D the 
molecular diffusivity of iodine vapour in air. D will be taken as 0-080 cm2/s 
at 20°c, following the results of Topley and Whytlaw-Gray (1927). 


Fig. 1 


SURFACE WIND SPEED SYMBOL. 
60b PAPER 38 CmY/SEC. 
. iu ke 


2 ns. 
50 COPPER : 


Nn 


ue) N= 0.332 0/5R2 


° 10 20 30 40 50 60 7 80 99 100 ite) 120 130 40 1SO. 
O/sR1/2 
Deposition of 1I vapour from airstream to flat plates. 


In fig. 1 the results are shown of plotting the values of NV, derived from 
the results of table 1 against 6! R1/2, where R is the local Reynolds 
number* and @ is the ratio of the kinematic viscosity of air to the 
diffusivity of iodine vapour. @ has the value 1-88 at 20°c. 


‘=e Aia.the Reynolds number characterizing the flow over the plate at 
v 


the point in question, and not the Reynolds number of flow in the tunnel as a 
whole. 


N 
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The points derived from experiments at different wind speeds lie 
approximately on the same curve, indicating that the deposition 1s con- 
trolled by aero-dynamical factors. ; 

Also in fig. 1 are shown the values of J, derived from experiments in 
which thin copper plate was used as the absorbing surface in place of 
filter paper. These points lie on the same curve as the filter paper points, 
except possibly at the upper end of the curve. 


Fig. 2 
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The straight line in fig. 1 is the line 
N ,=0°332 6/8 RY2_ PES ey he ES 

This is the theoretical (Pohlhausen) value of V,,, to be expected on 
laminar boundary layer theory on the assumptions : 

(a) that the vapour pressure of !°I at the surface is zero, 

(b) that the molecular diffusion is the transporting agency across the 
laminar layer. 

The agreement between theoretical and experimental values is good 
for small values of R'?, but as RV? increases the experimental values of 
N,,, increase more rapidly than R12, At the upper end of the velocity 
range the limiting value of R is approached at which the laminar boundary 
layer will break up and be replaced by a turbulent boundary layer. 


3.5, Deposition on Copper Cylinders 
Another series of wind tunnel runs were made with copper cylinders, 


varying in diameter from lin. to 0-028 in., and wind velocities from 
38 to 460 cm per sec. 
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The dimensionless number V,, is now given by 
Id 


N= cp? sees ts ls ec. (5) 


where I is the deposited activity per cm? averaged over the cylinder and 
d is the diameter of the cylinder. 

In fig. 2 WV, is shown plotted against 61/R1/2, In evaluating R, the 
effect of the obstruction caused by the larger cylinders in increasing the 
local wind velocity in the tunnel had to be allowed for. 


3.6. Comparison of Mass Transfer Numbers and Heat Transfer 
Numbers, for Different Surfaces 


In fig. 3 the graphs of V,, against 01/3R1/2 for 1811 deposition on flat 
plates and cylinders are compared with the analogous curves : 

(a) heat transfer by forced convection given by Hilpert, quoted by 
Goldstein (1937), 

(6) heat transfer and evaporation of water from spherical drops, 
given by Ranz and Marshall (1952). 


Fig. 3 
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In the heat transfer case the parameters NV, and @ are replaced by the 
corresponding heat transfer entities namely NV,, the Nusselt number 
and o the Prandtl number, which has the value 0-73 at 20°c. 

The curves are all quite close together, those for cylinders being between 
the straight line representing the sphere results and the curve for the 


flat plate. 
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The divergence of the experimental line for deposition on a flat plate 
from the theoretical line at moderate Reynolds numbers is seen as a 
tendency to approximate more nearly to the experimental curve for 
cylinders. 


$4. DEPOSITION ON SURFACES OF OTHER MATERIALS 


A number of experiments were made to compare the deposition on 
other surfaces with that on paper and copper. 

When flat plates 3in.xlin. were exposed to the same flow of 
‘ carrier-free ’ 131] vapour the deposition relative to that on a copper 
plate was as given in table 2. The level of activity on the copper was 
1500 d.p.m./em?. 


Table 2. 1211 Deposition on Flat Plates of Various Materials Relative to 
Deposition on Copper Plate 


Copper with 


Surface Paper | Steel he Glass | Perspex thin film 
ere paraffin wax 
Relative 
deposition 0-9 1 0-5 0-3 0-6 0-04 


The effect of slight differences in surface conformation, and position in 
the air stream, were such that these ratios have a possible error of about 
10%, and in particular, in the experiments with large plates, the copper : 
filter paper ratio was almost exactly unity, as shown in fig. 1. 

When additional inactive iodine carrier was added to the source at the 
rate of 20 micrograms per millicurie, and the surfaces again compared, it 
was found that the ratios of the deposition on the other surfaces to that on 
copper was reduced. Critical experiments to determine the saturation 
point in terms of number of iodine molecules per cm? have not been 
completed. 


§5. CoNCLUSION 


Copper and filter paper surfaces are found to act as ‘ perfect sinks ’ for 
the deposition of tracer amounts of carrier-free “I vapour. In these 
circumstances the deposition is controlled by aerodynamical factors, and 
is analogous to heat transfer by forced convection. The radioautograph 
technique enables fine detail of the deposition to be seen, and there is a 
possible practical application to the study of heat transfer problems. 
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° The idea that fission products deposition could be considered as a 
problem analogous to heat transfer was first suggested by Cockburn (1948). 
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By Wituiam Hume-Roruery, F.R.S. 
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[Received June 6, 1953] 


ABSTRACT 


The structures of the metallic carbides and nitrides of the type MX 
with the sodium chloride structure are discussed. The view that this 
structure is formed in order to provide the non-metallic atom with an 
octahedral interstice is shown to be improbable or only part of the truth 
and it is considered that the face-centred cubic arrangement with the _ 
metal atoms is favoured mainly because it provides the metallic atoms 
with mutually perpendicular bonds to six neighbours. The melting 
points are discussed, and the need for further work on the electrical 
properties of these compounds is emphasized. 


THE metallic carbides and nitrides with the Na-—Cl type of structure are 
of great interest in view of their hardness and high melting points. It has 
long been recognized that in these compounds there is some kind of 
homopolar bonding between the metallic and non-metallic atoms. 
These compounds were discussed by Brager (1941), who unfortunately 
did not take into account the effect of co-ordination number on the atomic 
radii. A more recent discussion by Rundle (1948) is of great interest, but 
appears to us to be only partly correct, and in the present note an attempt 
is made to remove some of the inconsistencies. 

The MX carbides and nitrides of high melting point have the sodium 
chloride structure, with the metal atoms arranged on a face-centred cubic 
lattice. The metallic elements concerned nearly all crystallize in the 
close-packed hexagonal, or body-centred cubic structures, and hence, 
as pointed out by Rundle, it is not justifiable to regard the MX carbides 
or nitrides as solid solutions of carbon or nitrogen in the crystals of the 
metallic elements themselves. The compounds may quite properly be 
regarded as interstitial compounds, with the small carbon and nitrogen 
atoms in holes between the larger metallic atoms, but some additional 
factor must be present to account for the taking up of a face-centred 
cubic structure by the metallic atoms. 

Rundle points out that, whereas in some cases (e.g. LaN) the metal- 
metal distances in the MX compound are almost the same as those in the 
metallic element, there are other cases where the former distance is the 
a 


* Communicated by the Author. 
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greater.* Rundle appears to regard this as in contradiction to the 
previously accepted views of these substances as interstitial compounds 
controlled mainly by the radius ratios, but it seems doubtful whether this 
criticism is valid. If the NaCl type of structure is regarded as composed 
of large and small spheres, the small spheres occupy the ‘ octahedral 
holes ’ between the large spheres, and a radius ratio of 0-41 corresponds 
to an exact fit. In the previous views of these compounds it was consid- 
ered that the NaCl type of structure was formed when the radius ratio 
lay between about 0-4 and 0-6, but it was not suggested} that the M—M 
distances were independent of the radius ratio. In the following table 


Table 1 
Interstitial diameter 
Atomic diameter Observed M—M for exact fit in octa- 
Metal for co-ordination distances in hedral holes of hypo- 
No. 12 compound thetical undistorted 
f. c. cubic metal 
Se 3°20 ScN 3°14 Test 
La 3-75 LaN 3-73 1-54 
Ce 3-64 CeN 3-54 1-50 
Pr 3-64 PrN 3-65 1-50 
Nd 3-64 NdN 3°64 1-50 
Ti 2-93 TiC 3°05 1-20 
TiN 2°99 
Zr 3°19 ZrC 3°32 1-31 
ZrN 3:27 
Hf 3:16 HfC (3-15) 1-29 
Th 3-59 ThC 3°75 1-47 
ThN 3°68 1-47 
V 271 Ae 3-03 LT 
VN 2:92 
Nb ; 2-94 NbC 3-16 1-21 
NbN 3°12 
Ta 2-94 TaC 3-14 1-21 
U ? 2-8-3:3 UC 3-50 1-15-1-35 
UN 3:45 


we show the atomic diameter for co-ordination number 12 of the metals 
concerned, together with the observed M—M distances in the carbide or 
nitride, and the atomic diameter of the sphere which would exactly fit 
the octahedral holes in the hypothetical undistorted face-centred cubic 
metals. These interstitial atomic diameters may be compared with the 
atomic diameters of C 1-54.4 and N 1-47 a from the normal co-valent bond 
length of the C-C and N-N linkages, which require uncertain corrections 


for co-ordination number and bond-number. As pointed out by Zener 
ne he ee ea ie ee eee 
* Rundle’s collected table contains no cases of a contraction greater than 
0-1 a, which is about the accuracy of the data. ' 
+ We do not, of course, imply that nobody suggested this, but it was certainly 
not the generally accepted view. 
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(1950), it is not justifiable to apply Pauling’s equation uncritically for 
corrections of this kind, because, quite apart from any theory of resonating 
co-valent bonds, a shortening of the apparent atomic diameter with 
decreasing co-ordination number is to be expected on grounds of energy 
alone. It is clear, however, that the cases of marked expansion in table | 
are those in which the ideal interstitial diameter is considerably smaller 
than the atomic diameter of carbon or nitrogen, and qualitatively the data 
are in agreement with the older views, except for UC and ThC where the 
expansion is abnormally great. The fact that uranium is an actinon and 
involves the complicated (7s, 6d, 5f) equilibrium may perhaps explain its 
abnormal behaviour. 

Rundle (loc. cit.) also pointed out the hardness and brittleness of the 
MX compounds, and suggested that this implied directed bonds, but this 
conclusion does not appear justified, except in so far as bonds may be 
directed in normal transition metals. Thus ruthenium (m.pt. 2550°C) is 
sufficiently brittle to be ground to powder in a percussion mortar at room 
temperature, and iron is known to become brittle at low temperatures. 
It is thus only to be expected that MX compounds with melting points 
of about 3000°c will be brittle at room temperatures. 


Table 2 
Supra conducting Sp. resistance at 
transition room temperature 
Substance temperature a Notes 


Cmpd. Metal Cmpd. Metal 


ScN — -- 3 —- 5°, impurity mostly 
oxide. «+ 

TiC 1-1? 1-1 1-8-2°5 0-8 0-28 impurity. 

TiN 2°6, 1-2 NER 1:3 0-8 2°% impurity. Oxide, 
alumina, carbon. 
The 1-2° suggests 
free Ti. 

vc _ ~ 1-5 1-6, 0-3 | Data for metal vari- 
able. 

VN 1:3 — 2 1-6, 0:3. | VINI2 O-1 or ryt 
002% Ni. 0:5% 

' insol. alumina. 

ZrN 3:2 — 1-6 0-4 ZrN 7-10% impuri- 
ties. 

NbC 10-1 8-2 1-75 ? | 0-1-0-2 | NbC. V impure. 

NbN — oo 2 0-1-0-2. | NbN no details. 

TaC 9-2 4:4 1-2 0-1 TaC fairly pure, traces | 


of Nb, Zr, Fe. 


1 

SeN and VC were shown not to be supraconductors at the lowest temperatures 
reached. NbN was not examined by Meissner and Franz, and the notes above 
refer to the paper by Friedrich and Sittig. 
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The MX carbides and nitrides are commonly said to exhibit normal 
metallic conductivity of the same order as that of the parent metals. 
This is probably true, but has not yet been proved conclusively for pure 
compounds of exact equiatomic composition. Thus, titanium exhibits 
supraconductivity with a transition temperature of 1-1°K, whilst in the 
work of Meissner and Franz (1930) TiC was shown to be non-supra- 
conducting, but to exhibit an inflection in the resistance-temperature 
curve at 1-1°K, which suggested strongly that free titanium was present. 
The data of Meissner and Franz are summarized in table 2, from which 
it will be seen that the position as regards purity and exact composition 
is most unsatisfactory. The notes in table 2 refer to the work of Meissner 
and Franz (loc. cit.) on supraconductivity. As far as can be seen, the 
values given by these authors for the specific resistances at room temper- 
ature are not from their own work. Some appear to be taken from the 
work of Friedrich and Sittig (1925) and there is seldom evidence that the 
exact composition MX was obtained. In spite of the higher melting 
points, the specific resistances of the compounds are greater than those of 
the constituent metals by factors which may be as high as 10. Inso far as 
the data are valid, one interesting point is the similarity between the 
specific resistances at room temperatures of all the carbides and nitrides 
between TiC and TaC in table 2. For all these the specific resistances lie 
between 1x 10-4 and 2x 10-4. 

According to Rundle, the carbon atoms give rise to mutually perpendi- 
cular bonds derived from one (sp) and two (p) orbitals, whilst the nitrogen 
atoms may exhibit the same behaviour, or more probably use three 
mutually perpendicular p orbitals. Rundle then regards the metal atoms 
as changing to the face-centred cubic structure, in order to provide 
‘ octahedral interstices ’ in which the carbon or nitrogen atoms will have 
six neighbours. This, however, appears incorrect, because the close- 
packed hexagonal structure itself contains ‘ octahedral interstices ’ to the 
number of one per atom, and so far as the environment of the carbon or 
nitrogen atoms is concerned, this could be obtained equally well in a close- 
packed hexagonal structure. The difference between the two structures 
lies in the environment of the metal atoms when the octahedral holes are 
filled by atoms of carbon or nitrogen. In the NaCl type of structure for 
the MX compounds, each metal atom has six interstitial neighbours 
arranged on three mutually perpendicular axes, whereas in an MX inter- 
stitial compound built up by filling the octahedral holes of a hexagonal 
structure of close packed spheres, only some of the MX bonds proceeding 
from an M atom are mutually perpendicular, and others are at an angle of 
about 70°. If the axial ratio of the close-packed hexagonal structure is 
slightly different from that for close-packed spheres [1/(8/3)= 1-633], the 
bonds to an atom in an octahedral interstice will not be exactly at right 
angles, but it seems improbable that this would invalidate the above 
argument. In any case, if we regard the environment of the carbon or 
nitrogen atom as the controlling factor, we shall be left with the problem 
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of why the metallic atoms prefer the arrangement of the cubical close 
packing rather than that of hexagonal close packing of spheres, parti- 
cularly since the latter could be obtained by such a slight alteration of the 
lattices of the actual close-packed hexagonal metals. We may thus 
conclude that, although the provision of octahedral holes for the carbon or 
nitrogen atoms may be one factor, the face-centred cubic arrangement of 
the metal atoms is favoured mainly because it provides the metallic atoms 
with mutually perpendicular bonds to six neighbours. This clearly 
suggests that d2sp? bonding from the metallic atoms is the main factor 
underlying the formation of these MX carbides and nitrides. 

We may now suggest the following qualitative picture of the MX 
carbides and nitrides with the NaCl type of structure. The interatomic 
distances make it almost certain that these involve covalent rather than 
ionic bonding, because the C++++ ion would be too small for an interstitial 
fit, whilst a C---- ion would be too large, and also too unstable. In the 
carbides we may regard each carbon atom as tending to acquire a share in 
an octet of electrons. Since the carbon atom itself contains four electrons, 
it must obtain on the average 2/3 of an electron from each of its metallic 
neighbours which give rise to resonating co-valent bonds. This is 
essentially the same view as that of Rundle, and in TiC, ZrC, and HfC, 
where the metal atoms are tetravalent, all the electrons of the metal atoms 
are needed to build up the octet, and the conductivity arises from the fact 
that only some of the hybrid (spd) orbitals are filled. On passing to VC, 
NbC and TaC, Rundle’s view is that the additional electron goes to form 
M-—M bonds, but the melting points in table 3 do not support this 
view. On passing from TiC—VC, there is a fall in melting point, and on 
passing from ZrC—NbC, and from HfC—TaC the melting point is unchanged 
(see below). 

In the nitrides, each nitrogen atom has five electrons and so requires 
half an electron from each neighbour to complete its octet. Unless 
charges are to be built up, the nitrogen atom will not contribute more than 
it gains, and consequently the bonding in the nitrides will be weaker than 
in the carbides, because the former involve three and the latter four 
electrons per atom. This is strikingly borne out by the data in table 3, 
which show that in every case there is a fall in melting point on passing 
from the carbide to the nitride of a given metal. At the same time we 
can understand why the M—M distance in the nitride is less than that in the 
carbide, because the atomic diameter of nitrogen is less than that of carbon. 
We can also understand why there are compounds ScN and LaN with the 
NaCl structure, but not ScC or LaC. Rundle’s view is that in passing 
along a sequence such as TiN-VN, the additional electrons are used to 
form M-—M bonds, until the latter are so strong that the metallic atoms 
are pulled into some other structure, but the melting point data for the 
nitrides do not support this view. The marked fall in melting point on 
passing from TiN + VN, and from ZrN -> NbN suggests that the additional 
electrons weaken rather than strengthen the bonding. 


Carbides and Nitrides of the Type MX 


We have next to consider why the MX nitrides and carbides are not 
formed by the elements of Group VI, and here it is interesting to note that 
WC and MoC have not structures in which the metal atoms considered 
alone are in close-packed assemblies. 
tures with the atoms arranged so that each metal atom is surrounded by 
six carbon atoms forming a trigonal prism, and each carbon atom is 
surrounded by metal atoms in the same way. The actual dimensions of 


the lattices are : 


MoC : a=2-90 
C=2-77 
cla=0-955 


Table 3 


Instead we find hexagonal struc- 


WC: 


cla=0-975 
The vertical side of the trigonal prism is thus nearly a square. 


Substance | Melting point °c Source 
TiC 3150 Goldschmidt 
Silsyy/ Emeleus 
TiN 2930 F. and S.* 
2947 Sidgwick 
2947 Emeleus 
ZrG 3500 Goldschmidt 
3532 Emeleus 
ZrN 2982 Emeleus 
ar 2982 Sidgwick 
2930 F. and S. 
HfC 3900 Goldschmidt 
Emeleus 
VC 2750 Goldschmidt 
2750 Emeleus 
VN 2300 Sidgwick 
2050 F. and 8. 
NbC 3500 Goldschmidt 
NbN 2300 Sidgwick 
2050 F. and 8. 
TaC 3900 Goldschmidt 
3877 Emeleus 
TaN 3087 Sidgwick 
3087 Emeleus 
2800 F. and S. 


(eS Se ee ee 


* F. and S.=Friedrich and Sittig (1925). 
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If the transition elements are considered as a whole, we know that on 
passing along any one Long Row of the Periodic Table, the (n— 1)d states, 
which are above the ns states in Group IA pass to relatively lower energies 
until at Group IIB the (n—1)d states are below the ns states. As we pass 
along the Period, we may therefore expect a tendency for the d part of 
the hybrid bonds to increase to a maximum, and then to diminish. We 
have suggested that in Groups IV and V, the MX carbides and nitrides 
result from octahedral d2sp? bonds, and it is of interest to note that the 
trigonal prism arrangement of six neighbours corresponds to d‘sp or d°p 
bonding (see Eyring et al. 1944), and we may suggest that this type is now 
involved. Each carbon atom has again six neighbours, and must absorb 
2/3 electrons per bond in order to complete its octet. If this view is 
acceptable we can perhaps understand why WC has a lower melting point 
than metallic tungsten. The latter crystallizes in the body-centred cubic 
structure in which each atom has six valency electrons to bind it to eight 
close neighbours, whilst, if our picture is correct, the cohesion in WC results 
from each tungsten atom using four electrons to bind itself to six carbon 
atoms. The fact that there is no analogous compound CrC would agree 
with the scheme of valencies proposed by Hume-Rothery, Irving and 
Williams (1951), who pointed out that in the chemistry of these metals, 
the higher valencies were favoured in the later Periods. 
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ABSTRACT 


Measurements have been made of the ferroelectric hysteresis and di- 
electric constant of small single crystals of barium titanate from 20°c to 
150°C. | 

In untwinned crystals, free from external constraint the results clearly 
indicate that the ferroelectric transition occurring in the region of 120°c is 
a thermodynamically first order transition. Over a small range of 
temperature above the zero field transition electric hysteresis is reintro- 
duced by the application of a sufficient electric field. The forms of the 
hysteresis figures in this temperature range are qualitatively explained 
in terms of Devonshire’s treatment of the ferroelectric effects. 

In twinned crystals, optical and electrical measurements show that it is 
not possible to fix accurately a transition temperature for the whole 
specimen, as different regions transform at different temperature. It is 
suggested that the behaviour of these specimens may be due to imper- 
fections in the crystal lattice, which give rise to local stresses in the region 
of the transition. 


4 
$1. [yTRopUCTION 


Untit recently, the ferroelectric transition which occurs in barium 
titanate in the region of 120°c has been considered to be a thermodynamic- 
ally second order effect (Forsbergh 1949, Megaw 1947), though it was 
recognized that the changes in the dimensions of the unit cell which occur 
immediately below this temperature, are more rapid than would be 
expected for such a 4 type transition. 

In a phenomenological treatment of the ferroelectric effects in barium 
titanate, it was suggested by Devonshire (1949) that the available data 
was more consistent with a first order change. Recently, several further 
experiments have been carried out which tend to confirm this suggestion 
(Kanzig 1951, S. Roberts 1952) but no evidence appears yet to have been 
produced of any direct measurement of the discontinuous changes which 
occur in the electric polarization in a first order transition. 

Over the past three years, we have made extensive measurements of 
optical and dielectric effects in both twinned and untwinned crystals of 
barium titanate (Cross 1952). As some of these results give direct 
amr co nt Wt her eae eerie te 
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evidence of the polarization changes in the region of the transition, which 
are not in agreement with those of previous observers, it seems desirable to 


place them on record. 


§ 2. EXPERIMENTAL DETAILS 


Previous optical measurements (Kay and Vousden 1949, Forsbergh 
1949, Cross et al. 1949), have shown that the internal twinning in the 
ferroelectric tetragonal phase of barium titanate is easily changed by 
mechanical, thermal or electrical treatment. To make any correlation 
between the electrical properties and the internal structure it is thus 
desirable to maintain a constant check on this feature. In these measure- 
ments this was accomplished by making simultaneous optical observations 
through semi-transparent electrodes. The individual twin components 
are uniaxial, and in the simpler specimens it is easy to determine the 
distribution and orientation of all the components by an examination in 
polarized light. 

Measurement of the electric hysteresis was carried out by a modified 
‘Sawyer and Tower ’ method (Von Hippel e¢ al. 1946), and the dielectric 
constant values were deduced from the capacity of the specimen measured 
on a modified Maxwell bridge. For small values of the exciting field, the 
specimen behaves as a linear condenser, and the bridge may be perfectly 
balanced. If the measuring field is increased, no complete balance may 
be obtained, as the curvature of the polarization field curve introduces 
higher harmonics in the specimen polarization. An approximate balance 
may, however, be obtained for the component of the polarization in phase 
with the fundamental bridge frequency. 

The value of the dielectric constant for weak fields we have called ep, 
and the values deduced for higher field strengths e,. 

The crystals used were small flat plates of average dimensions 
1mm x 1 mmx 0-05 mm, only the more perfect rectangular shapes being 
selected. Semi-transparent gold electrodes were deposited on the larger 
faces by evaporation at reduced pressure. Contact was made onto the 
plated faces by lightly pressing platinum probes, which in no case exerted 
sufficient pressure to effect any optically detectable change in the specimen. 

For work at elevated temperature, the specimen holder was placed in a 
well lagged viewing chamber, on the stage of a polarizing microscope, and 
heated by a stream of hot air. The temperature was recorded on a 
copper constantan thermocouple in contact with the specimen, and 
temperatures up to 150°c could be maintained with a regulation of better 
than -+-0-5°c. 


§ 3. OBSERVATIONS 


To allow direct intercomparison of the performances of different 
specimens, the electric hysteresis figures have been re-plotted from the 
original oscillograms with the scale suitably adjusted. Where the 
spontaneous polarization is given, this has been deduced by extrapolating 
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to zero field the linear regions of the hysteresis curve of a saturated 
specimen. 

3.1. Untwinned Crystals 


All the specimens used in these measurements had their tetragonal 
¢ axis perpendicular to the plane of the crystal plate. 


Fig. 1 
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Electric hysteresis in a typical untwinned crystal at increasing temperatures 
(a) 24°c, (b) 120-5°o, (c) 129°o, 


A typical hysteresis figure, taken at 24°c is shown in fig. 1 (a). On 
elevating the temperature, only minor changes occurred in the hysteresis 
figure up to the transition temperature. Figure 2 shows the spontaneous 
polarization obtained by extrapolating the hysteresis curves in this region. 

At temperatures from 1°c to 5°c below the transition temperature, 
bright birefringent bands usually appeared, due to transient internal 
twinning. From the hysteresis figures, however, it was evident that this 
had a negligible effect on the polarization. 


4G2 
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In the absence of an electric field, the crystal of fig. 1 (a) transformed 
suddenly at 120°c, losing all its birefringent areas at this temperature. 
Immediately below the transition there was a normal hysteresis curve. 
Just above 120°c, a peculiar kinking occurred at the coercive field, fig. 
1 (b), and this became more marked as the temperature was raised. At 
123°c, a loss free region began to develop in the centre of the curve, and at 
129°C a well marked almost linear region was evident (fig. 1 (c)). This loss 
free region continued to extend as the temperature was raised, until at 
135°C, with a field approaching the breakdown strength of the crystal, only 
the tips of the curve showed hysteresis. 


- a | 
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Spontaneous polarization P, as a function of temperature in untwinned crystals. 


An optical examination at any temperature in the above range showed an 
interesting series of changes as the electric field was increased. At zero 
field the crystal appeared completely black, and was obviously above its 
transition temperature. On increasing the field, no change was evident 
over the loss free region of the hysteresis curve, but immediately the field 
attained a sufficient value to excite hysteresis, bright birefringent bands 
flashed out, similar to those observed in the field free crystal just below the 
transition, and these were again cleared as the field was further increased. _ 

The variation of reversible dielectric constant ¢, with temperature, 
measured with an alternating field of 10 volts/em and without any steady 
biasing field, is shown in fig. 3. The curve of in-phase dielectric constant 
<, against temperature, for a measuring field of 1300 volts per em is shown 
in fig. 4. The manner in which ez varies with temperature when a steady 
biasing field is also applied has been determined, and shows an elevation 
of the transition temperature by electric fields of 0-9°c/kv/em, the 


graphs are not, however, reproduced here as similar results have already 
been published (K anzig loc. cit.). 
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3.2. Twinned Crystals 


It is difficult to summarize completely the very varied performance of 
the twinned crystals, but the following features were generally evident. 

(a) Even for a simply twinned crystal containing only a single twin 
boundary, the coercive field was very much greater than obtained from an 
untwinned specimen under similar conditions. 


Fig. 3 
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(b) The crystals were difficult to saturate electrically, but even where the 
optical properties showed that the tetrad axis was completely oriented in 
the field direction, the extrapolated spontaneous polarization was 
considerably smaller than for an untwinned crystal. 

(c) Optical observations showed that the transition in the region of 
120°C was not discontinuous, but often the specimen contained a patch- 
work of transformed and untransformed regions for a range of temperature. 
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_ Hysteresis effects in twinned crystals. 
In a simply twinned crystal (a) 20°c, (b) 109°c, (c) 112°c. 
In a crystal with more complex twin structure (d) 15°c, (e) 115°c. 


In some crystals, the change could be made to occur behind a moving 
boundary, and by maintaining a constant temperature, this boundary 
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could be halted and the two phases retained indefinitely in the same 
specimen. 


Some typical hysteresis curves for various types of twinned crystal are 
shown in fig. 5. 


§ 4. COMPARISON BETWEEN EXPERIMENTAL AND THEORETICAL VALUES 


It has already been indicated that the phenomenological treatment first 
applied to barium titanate by Devonshire (1949), assumes a first order 
transition in the region of 120°c. This type of change is a direct conse- 
quence of the expression chosen for the free energy of the unstressed crystal. 
In the original paper expressions are derived for the variation with temper- 
ature of spontaneous electric polarization, and of weak field dielectric 
constant in terms of four constants, 7’, the transition temperature, 7’, the 
extrapolated Curie temperature, P, the spontaneous polarization at 7',, 
and yx ' the reciprocal electric susceptibility at 7’. 

The effects of an external electric field may be considered by adding a 
further term to the expression for the free energy (Devonshire’s eqn. (6.1)). 
For a field H applied in the z direction the additional term will be —#.P,,. 
In the region of the transition temperature Px—=Py=0, and the necessary 
condition for a stationary value of the free energy may be shown to be 


B=a (t—4a?-+ 3a4), ee ae 
where B=H/P1xo', a=P,/Py, t=(T—T9)/(T1—T). 
From (1) the electric hysteresis figure at zero frequency may be 


computed for any value of temperature. 
Taking the values 


P,=15-4 » coulombs/cm?, 


C 120°c, 
from our measurements on untwinned crystals, and 
J = 1107 6, 


Xot=1-0x 10-3, 
from the measurements of 8S. Roberts (1949), the theoretical values of 
polarization and dielectric constant have been plotted on fig. 2 and fig. 3, 
and the relation between polarization and electric field for three temper- 
atures above the transition temperature in fig. 6. 

It is evident from a consideration of the energy, that in regions such as 
AB in fig. 6 (b) the free energy is a maximum. [If it is assumed that the 
system must have sufficient energy to surmount this barrier there will be 
irreversible polarization changes at A and B and the dotted hysteresis 
curve will be traced out. In practice it is probable that inhomogeneity 
in the specimen will facilitate a reversal of the polarization of lower energy. 


§ 5, Discussion or RESULTS 


Devonshire’s treatment clearly gives a satisfactory correlation of the 
polarization measurements taken on the untwinned crystals. With the 
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simplifications which are necessary in such a treatment exact agreement, 
with the experimental values is not to be expected, but the general form 
is well reproduced. In particular, the first order change in the polarization 
is confirmed, and the hysteresis above the transition temperature 
explained. 

The reintroduction of the ferroelectric state by a sufficient electric 
field above the transition temperature shown by the hysteresis effects is 
also evident in the optical and dielectric measurements. The sudden 


Fig. 6 
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Hysteresis figures above the upper transition temperature from Devonshire’s 
treatment (a) 121°c, (b) 124°c, (c) 129°o. 


increase in polarization in the ferroelectric region produces an electro- 
strictive distortion, transient internal twinning occurs, and the change is 
evident optically by bright birefringent regions which appear in the crystal 
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ry 


In dielectric measurements, the effect is most easily observed in measure- 
ments of the temperature variation of « at constant field, the sharp peak 
in €p which occurs at 120°c being moved to higher temperature. Using 
this method Kanzig (loc. cit.) finds an elevation of 1-2°c/kv/em, and 
similar measurements on our crystals give a value of 0-9°o/kv/em. 

Further confirmation of the discontinuous transition in the polarization 
is provided by the measurements of <;. The very high values of ¢, just 
below the transition are clearly associated with the reorientation of spon- 
taneously polarized ‘domains’, and the discontinuous change in «, is 
therefore further evidence of a discontinuity in the polarization. 

The continuous change observed in ¢p in the region of the transition 
temperature (fig. 3) is not in agreement with the postulated polarization 
change, but it must be remembered that even in the best single crystals 
there is transient internal twinning in the region of the transition. This 
inhomogeneity in the crystal must influence the dielectric properties, and 
is probably sufficient to account for the observed effects. 

The only previous measurements of polarization in untwinned crystals 
(Merz 1949) show a generally lower value of polarization than our results, 
with a rapid but continuous change in the region of 120°c. There is, 
however, some evidence in the shape of the temperature variation of 
polarization to suggest that saturation is incomplete, and in certain speci- 
mens hysteresis curves similar to our fig. 1(c) were obtained. 
Unfortunately simultaneous optical measurements were not carried out, 
and subsequent piezoelectric measurements (Caspari and Merz 1950) on 
similar specimens indicate the presence of some internal twinning. 

No complete explanation of the different properties of the twinned 
crystals will be attempted here. Piezoelectric and dielectric measurements 
suggest that the individual twin components are electrically neutral, and 
it is therefore unlikely that the twinning is a ‘ domain ’ formation necessary 
to eliminate free poles. Forsbergh (loc. cit.) has suggested that the 
twinning is due to lattice imperfection, and has demonstrated that in 
certain impure crystals, the spontaneous distortion in the ferroelectric 
region is utilized by suitable twinning to relieve strain caused by a local 
concentration of impurities in the specimen. 

This suggested association of the twin structure with defects in the 
erystal lattice provides a qualitative explanation for our results. It may 
be shown from Devonshire’s treatment that constraint applied to the 
specimen is likely, both to reduce the transition temperature, and by 
modifying the coefficient of the terms in P* in the free energy, to reduce the 
discontinuous change in polarization at this temperature. Thus to explain 
the observed effects, it is suggested that the lattice defects give rise to 
local constraints, which produce inhomogeneity in the dielectric 
performance in the region of the transition. 

Clearly, although this suggestion gives some explanation for the 
differences observed in the dielectric performance, quantitative calculations 
of the magnitudes involved could only be made after much more detailed 
information of the behaviour of the twin boundaries has been obtained. 
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It is hoped to obtain further information of the properties of twinned 
crystals, by examining the electric hysteresis at much lower frequency, 
where optical tracking of the twin boundary motion should be possible. 
Also a further examination of the hysteresis above the transition is being 
carried out, using measuring fields of higher frequency, to obtain inform- 
ation of the relaxation process between ferroelectric and unpolarized 
states. 
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ABSTRACT 


It is shown that commutation relations are derived consistently from 
the quantum dynamical principle, for systems obeying first order equations 
of motion. 


§ 1. INTRODUCTION 


THE author (Schwinger 1951) has developed a formulation of quantum 
dynamics which enables one to deduce all properties of a quantum 
mechanical system from a single dynamical principle. It has been 
asserted, in a paper published in this journal (Burton and Touschek 1953), 
that the scheme fails for systems obeying first order equations of motion. 
That statement is not correct. 

In order to clarify the origin of this misunderstanding, we shall present a 
brief discussion of such systems with a finite number of degrees of freedom. 
It should be realized that the employment of first order equations of motion 
is not a restriction but rather a standardization of systems that obey the 
conventional determinism principle—that knowledge of the values of a 
complete set of compatible physical quantities at a given time constitutes a 
state of maximum information which permits the determination of such 
states at any other time. The differential equations of motion for systems 
of this type must be of finite order, and can therefore always be presented 
as first order equations by suitable adjunction of variables. 


§ 2. Tom DyNAMICAL PRINCIPLE 


It is the object of quantum dynamics to construct transformation 
functions of the type 


(C1! ty [e2” te)=(¥ (Cy t,)" Y (Co" ta) sega yen at) 
where ¢, and ¢, are complete sets of commuting Hermitian operators. 
The fundamental dynamical principle is a differential characterization of 
such transformation functions. Any infinitesimal change can be repre- 
etl by 8(0y ty oo” to)=(t/h) (Cy ty | Wie oo” te); 
which is the definition of the infinitesimal operator 5W,,. The composi- 
tion and reality properties of transformation functions imply that 


OW 13=8 W12+8 Ws, 
5 Wi2'=8 W 40. 


* Communicated by the Author. 
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The basic postulate of the dynamical principle is that there exists a class of 
alterations for which the associated operators 5W,, are obtained by 
appropriate variation of a single operator, the action operator W x2; 


dW 12=8( W 19). 


According to the additivity property of the action operator, it must 
possess the form 


ty 
Was | dt L{t], 
to 


where L[t] is a Hermitian function of fundamental dynamical variables 
a,(t), in the infinitesimal neighbourhood of the time t. There is no loss of 
generality in requiring that the ~, be Hermitian. 

For a given dynamical system, that is, for a given form of the Lagrangian 
operator L, we can alter the dynamical variables x,, and the times ¢, and 
t,. This must correspond to the freedom in description of a given 
dynamical system, the possibility of introducing infinitesimal changes in 
the states occurring in the transformation function (1)—infinitesimal 
alterations of the complete sets of operators, or of the times to which they 
refer. These variations are localized at the times ¢, and ¢,, and are 
described by infinitesimal unitary transformations of the eigenvectors 


Ou (C2” ty)=— (a/h) Gs abe (f,” ty), 
Ont (¢y' ty)'=(i/h) ye (fy t,)" Gy, 


where G, and G, are infinitesimal Hermitian operator functions of the 
dynamical variables at ¢, and t,, respectively. Thus, for a given system, 


bW2=G4,—4 > i ee eee 
which is the quantum form of Hamilton’s principle. It asserts that 
W,, is stationary with respect to variations of the dynamical variables in 
the interval between ¢, and ¢,, since G, and G, only contain variables 
referring to t, and tz. This action principle implies equations of motion 
for the dynamical variables and yields expressions for the generators of 
infinitesimal unitary transformations, G, and G,. 
Note that two Lagrangians differing by a time derivative describe the 
same system. Thus 


L=L—(d/dt)W, W=W(z, t), 


yields Va Waa (Vl): 

Accordingly, the action principle for Wis is satisfied if it is obeyed by Wo. 

That is, W,,and W,, imply the same equations of motion, but 
8W1.=G,—G,, 

where 5W,=G,—G,, dW.=G4,—G,, 

defines new generators of infinitesimal transformations at t, and fp. 


The latter equations are a special case of (2), in which all operators refer 
to the same time. 
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In order that the Lagrangian yield first order, explicit equations of 
motion, it must possess the general form 


at, dz, 
1 d. d 
=5 (x4 s — 74°) —H (x,t), 


in which the time derivative term has been symmetrized with respect to the 
operation of integration by parts, corresponding to the freedom of adding 
a time derivative to the Lagrangian. If LZ is to be a Hermitian operator, 
this property must apply to H, the Hamiltonian .operator, and the 
numerical matrix A must be skew-Hermitian. 

The action operator is 


ty 
Wee | [eed de—de An)—Hdt] 
ts 


in which form the limits of integration are objects of variation. We can 
introduce an auxiliary variable 7, which is such that the variations S¢, and 
dt, are represented by a change in the functional dependence of ¢ upon vz. 
This enables us to vary the x, and ¢, with fixed limits 7, and 7,. The 
auxiliary variable 7 need not be written explicitly, since it is not subject 
to variation. Hence 


5W4.= J [6vA du—dx Adxu—dH dt+dH St] 
+ fd[4(vA da—du Ax)—HSt]. 
The action principle now asserts that 


dxA dx—dxA 6u=6H dit—dH 5t, 
or 


dH dx da 


and supplies the information that 
G=4 (xA da—dx Ax)— Hot, 


Adz, Se es) 


at times ¢, and ¢,. 

The Hamiltonian operator is an arbitrary function of the x,. If its 
variation is to possess the form (3), with dz, appearing only on the left 
and on the right, these variations must possess elementary operator 
properties, characterizing the class of variations to which the dynamical 
principle refers. Thus, we should be able to displace 6x, entirely to the 
left, or to the right, in the structure of 6H, 


SH —(0H/at)8t=8a(d,H/dx)=(0,H/Ax) Sx, 


which defines the left and right derivatives of H. In view of the complete 
symmetry between left and right in the multiplication process, we infer 
that the terms in (3) with 5a on the left, and on the right, are in fact equal. 
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We thus obtain dH/dt=0H/dt, 2 SS rr ae 
and 


2A (dx/dt)=0,H/ox, 
—(dx/dt)2A = —2A* (dx/dt)=0,H/ox, 
where the latter must be equivalent forms of the equations of motion. 
Similarly, 
G=—dxAx— Hot 
and 


G=xAdx—Hdt—(A*x)dx—Ht 


must be equivalent’ forms of the infinitesimal generator G. 

If the matrix A completely decomposes into two submatrices, with a 
corresponding decomposition of the dynamical variables into two sets, we 
speak of the two sets of variables as being kinematically independent, 
since the derivative term of the Lagrangian is additively constructed from 
the two sets, as is the infinitesimal generator G, apart from the term 
—H6&t. (Should the Hamiltonian also decompose in this manner, the 
two sets are dynamically independent.) Now the matrix A can always be 
written as the sum of an anti-symmetrical, real matrix a, and of a 
symmetrical, imaginary matrix s. In view of the manner in which A and 
the transpose of A enter in equivalent versions of the equations of motion, 
and of G we anticipate that a consistent theory can be constructed only if 
the expression of A as 

A=a-+s 
is a complete decomposition. Accordingly, there exist two kinematically 
independent sets of dynamical variables ; variables of the first kind, 2,, 


associated with a, and variables of the second kind, €,, associated with s. 
The Lagrangian now becomes 


b= 340, «at+3|6 °F | —H(ax, &€, t). 


The equations of motion appear in two sets, 
2a(dx/dt)=0,H/dx=0,H/dx, Ry Pid ca. (tH 
and 
2s(dé/dt)=0,H/0é=—0,H/d€, “<cghiniy Mesa ae 
while the infinitesimal generator reads 
G= —da(ax)—8&(sE)— H8t 
=—(ax)du-+ (s€)5E—H8t. 
Concerning the variables of the first kind, we see fron. the identity of the 
left and right derivatives of H that the 5x, must commute with all 


operators. This is consistent with the equivalent expressions for G. 
Hence 


[d2:;, &;|=[82;, é,J=0, 
which is +o say that the 5x, are infinitesimal multiples of the unit operator. 
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On comparing the two forms of G, we infer that the 5, anti-commute with 
the variables of the second kind, 


10g, é e}=0. 
The opposite signs of the left and right derivatives of H with respect to 
variables of the second kind implies that H must be an even function of the 
é,, but it remains an arbitrary function of the x, provided that 

[6 ce v ]=0. 
We have, in these operator properties of 5x; and 5€,, the precise specific- 
ation of the class of variations to which the fundamental dynamical 
principle refers. 

We shall assume, for simplicity, that the equations of motion can be 
solved explicitly for the dz,/dt and dé,/dt which requires that the matrices 
a and s be non-singular. In order that this property holds for the anti- 
symmetrical matrix a, the number of variables of the first kind must be 
even (2n,). We shall later present an argument indicating that the 
number of variables of the second kind must also be even (2n,). The 
separation of the 2n variables of either type into two sets of ” variables 
leads to the canonical version of the formalism. 


§ 3. VARIABLES OF THE First KIND 
We employ the fact that, by means of real linear transformations of the 


variables, the anti-symmetrical real matrix 2a can be brought into a 
diagonal form in which the n, diagonal elements are the two-dimensional 


matrices 
== 
1 On x 


If the dynamical variables associated with the kth two-dimensional 
matrix are called g, and p,, the relevant time derivative term in the 
Lagrangian becomes 


22, (312 wii ae : 


the equations of motion read 

dq,,/dt=0H/0p,, —4dp;,/dt=o0H/0q;,, Jn hh SOOT AN UG 
while the contribution to G arising from variations of the variables of the 
first kind is 


4 


Ga, x2 ee (Pr 9%%—S8P x Tn) ee Liners Ao) 


Less symmetrical but more convenient forms are obtained by adding a 
suitable time derivative to the Lagrangian with a corresponding alteration 
of the generator G. Thus, with 


W=—424{px, Ux}; 


k 
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the derivative term in L becomes 
1 dq, \ 
2 5) {Pe a? 
and the contribution to the infinitesimal generator is 
Grae Pr Ce 


If we reverse the sign of W, these quantities are 


and 
Gy= ye SHE a semi OP pe | 


The advantage of these forms is the unambiguous interpretation attached 
to the infinitesimal generators G, and G,. Evidently G, describes that 
infinitesimal unitary transformation in which the operators q; are 
changed by 8q;, ({,={—54;,), With no alteration of the operators p;, and 
similarly for G,. The expression of this for an arbitrary function F of the 
canonical variables, q;, and p,, is 


(L/n)[P, Gale (OF /09;) 9x: 
and (L/th)[F, G,J=2 (OF /ep,) Op, 


or, remembering that the 6q;, and dp;, commute with all operators, 

(L/h) [F, py]=OF/09;,, 

(1/th) lq, FJ= oF /Op;,, 
from which we obtain the commutation properties of the canonical 
variables 


[Vx VJ=lPx: pi jJ=9, 


[Vis Pil=h Oxi 
The evident interpretation of 
G j= — Host 


as the generator of the infinitesimal unitary transformation from variables 
at time ¢ to variables at time ¢+-d¢ supplies the important verification of 
the consistency of the information concerning equations of motion and 
generators of infinitesimal unitary transformations. We have 

(1/th) LF, G,J=—(dF/dt— oF /dt)8t, 
or 

dF /dt=oF /ot+-(1/th)[F, H]. 
By successively placing F=H, q;,, p;,, we deduce (4) and the canonical 


equations of motion (7) from the properties of infinitesimal 
transformations. 


We now remark that the generator (8) is 
Ga o=3 G, 2 Gs 
which clearly produces the transformation in which the q, and p;, are 
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changed by 469, and }5p,, respectively. Since this statement refers 
symmetrically to the two sets of canonical variables, the same interpre- 
tation applies to the generator which employs no special form of the 
matrix a, 

G,,=— 6x az, 


namely, that it describes changes of x by 48x. The equivalent 
commutation statement is 
(1/th)[2ax, F|=dF/ox, 

which yields equations of motion in the form (5), and exhibits the 
commutation properties of the variables of the first kind in the single 
equation : 
(v7, @ == + (a) 3. eg mele h fone Te 9) 

The oversight that occurs in the paper of Burton and Touschek (1953) 
is the treatment of the 2n, variables x as canonical coordinates, and the 
linear combinations of these same variables, —ax, as canonical momenta. 
This results in a commutator which lacks the factor 1/2 of (9). 


$4. VARIABLES OF THE SECOND KIND 


It is convenient to anticipate that the real, symmetrical matrix (1/ih)s 
must be positive—definite. Hence there exists a real, linear transformation 
of the variables which reduces this matrix to the unit matrix. The 
generator of changes in the variables of the second kind then reads 


ema IoC es hy 2) (10) 


which shows that the various €, are kinematically independent. The 
operator expression of this statement is 


es &, ) &,/=0, =D; 
which implies that 
Leas E,}=0, asp, 


in virtue of the anti-commutitivity of 5€, with €,. 

Now, if the number of variables of the second kind were odd, the product 
of all these operators would commute with every €,, and would therefore 
be a multiple of the unit operator, contradicting the independence of the 
various €,. An alternative version of the argument stems from the 
requirement that the set of operators €, be algebraically complete, in the 
sense that an operator possessing the properties of 5£, can be constructed 
from them. For an even number of variables, the product of all the €, 
anti-commutes with each €,, but no such operator exists for an odd number 
of variables. We conclude that there must be an even number of 
dynamical variables of the second kind. 

To introduce the canonical formulation for variables of the second kind, 
we observe that, in virtue of the anti-commutitivity property of the 6€,, 
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a pair of terms in (10) can be combined in the following way : 
ih (E, 8E,+£ So) =thF [(E1—tE a) 8(E1 +182) —8(E1 —1F 2) (Fr +089) 

=3[pdq— pq]. 
Here p and q are non-Hermitian operators, which are related by 

p=thq', 

if we choose g as £,+7£. On pairing the 2, terms of (10) in this manner, 
we obtain the canonical form of the contribution to G associated with 
variables of the second kind, 


Gest 2 (7, 80.— 80¢ Ja) + as Bae 


The derivative term in the Lagrangian and the equations of motion, 
expressed in canonical variables, are 


ly, 1 qa, dp, 
22 \3| 2° a atc 
and ag, |dt-=2 HO. oi ee 
On adding suitable time derivatives to the Lagrangian, we obtain 
modified forms of the time derivative term and the associated generator, 
namely 
a 31 Pas dq,/dt}, 


G = 2 Da 84x 
and —2' 3[dp,/dt, q,], 
i G,= ae SDy, Qa==' Va SDs 


According to the evident interpretation of the infinitesimal generators 
G, and G,, we have 

(L/AR)LF, GJ=% (0,F/09) 89a 
and (1/ii)LP, G,]=E 5p,(,F /Ap,)- 
[fF contains an odd number of variables, of the second kind, these relations 
imply that 

(l/h), p.j=0,F/0q,=0,F/09,, 

(1/th) {93 F}=0,F /0p,=0,F /dp,, 


while, for operators constructed from an even number of second kind — 
variables, 


(1/ch)[F, Pa]=0,F/0q,.= —0,F/0q,; 
(1/th) Mas F)= 0,F/0p,—=—d rl/Op,. 


From the first set, we obtain the commutation properties of canonical 
variables of the second kind, 


{Jas I= {Px Pe=9, 
{dar Dp} =th B yp. 
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As applications of the second set of relations, we deduce that variables 
of the first kind commute with variables of the second kind, and verify 
that the equations of motion for second kind variables follow from the 
properties of infinitesimal transformations, thus ensuring the consistency 
of the theory for both types of dynamical variables. 

Since the generator (11) produces the transformation in which the q,, 
and p, are changed by 46, and $6p,, respectively, the same interpretation 
applies to the generator which employs no special form for s, 

G-=— 8€ 8, 

namely, that it describes changes of € by $5. Hence, for operators F 
which are, respectively, even and odd functions of the variables £, we have 

(1/ih)[28§, F]=0,F/0é=—0,F/0€, 
and 

(1/ih) {2s€, F}=0,F /0é=0,F/0€. 

The first relation, applied to the Hamiltonian, yields the equations of 
motion in the form (6), while the commutation properties of the Hermitian 
variables of the second kind, €,, emerge from the second equation, 


{&,, &g}=ih4(s) yp. sel beyioebare ooh seems 9) 
Incidentally, the structure of the left side, in this result, provides the basis 
for the statement that (1/i%)s must be positive—definite. In their dis- 
cussion of this type of system, Burton and Touschek regarded the 2n, 
variables € as canonical coordinates, and the linear combinations of these 
same variables, s£, as canonical momenta, which produced an anti-com- 
mutator that lacked the factor } of (12). 

We shall carry the discussion no further here, since our objective was to 
show that the dynamical principle yields consistent commutation pro- 
perties for systems obeying first order equations of motion. The reader is 
referred to the second paper of a series by the author (1953), ‘‘ The 
Theory of Quantized Fields ”’, in, which fields obeying first order equations 
are considered. 
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CXXV. Schwinger’s Dynamical Principle 


By W. K. Burton and B. F. TouscHEK 
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ABSTRACT 


A comparison is made of the methods applied by Schwinger and by 
the authors to the derivation of commutation relations for systems whose 
equations of motion are of the first order in the time derivatives. 


In the foregoing note Schwinger (1953) criticizes our paper (Burton and 
Touschek 1953) on the commutation relations in Lagrangian quantum 
mechanics. A misunderstanding as to the scope and the contents of this 
paper seems to have arisen. 

We did not wish to suggest that Schwinger’s dynamical principle as 
such was wrong ; in fact we used it in our treatment of first order systems. 
Our intention was to point out an inconsistency arising from an uncritical 
application of this principle in the form in which it was published by 
Schwinger (1951), and to show how this difficulty could be avoided. Our 
method consisted in the consideration of variations of the variables 
induced by variations of the sources, a method originally introduced by 
Peierls (1952). In this way we obtained by applying Schwinger’s principle, 
and without recourse to the canonical formalism, the correct commutation 
relations for the systems under consideration. In field theory, this 
method gives the covariant commutation relations. 

In his note Schwinger gives a discussion of a generalized form of the two 
examples investigated in our paper, showing that by using canonical 
variables the correct commutation relations can be recovered from his 
formalism. However, the necessity for introducing these variables was 
not suggested in the original (1951) form of the theory. In the case of 
variables of the second kind for which no correspondence with classical 
mechanics exists, the introduction of canonical variables is by no means 
trivial. In fact, previous to the publication of our paper we discussed 
at a conference at Glasgow in July, 1952, the introduction of the variables 
d= +42 and q+=q,—iq, (Schwinger’s canonical variables of the second 
kind) in place of the Hermitian variables g, and q., but we discarded this 


procedure later in favour of the Peierls method, which we felt to be more 
cogent and powerful. 


a a ee 


* Communicated by the Authors. 
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The combination of the Schwinger principle and the Peierls method 
seems to us to have the advantage that it is not necessarily restricted to 
Lagrangians which allow a canonical formulation, so that our method 
seems more in keeping with the spirit of Schwinger’s dynamical principle 
than the recourse to canonical variables. 
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COMMENT ON THE ABOVE NOTE 


There is no objection to deriving commutation relations from the 
consideration of source variations. Indeed, the use of external sources has 
always been an integral part of the technique associated with the dynamical 
principle (see (1951) and (1953) which, incidentally, was submitted before 
the paper of Burton and Touschek became available). But one cannot 
ignore the commutation properties implied by the variations of the 
dynamical variables. The theory must be consistent, and this requirement 
applies to the consequences of all the variations to which the fundamental 
dynamical principle refers. 
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Measurement of the Elasticity of Solid Argon by U ltrasonic Methods 


By J. R. Barxer, E. R. Doses and G. O. JoNEs 
Department of Physics, Queen Mary College, University of London 


[Received July 16, 1953] 


THE importance of a knowledge of the elastic constants of the solidified 
inert gases for the development of the theory of interatomic forces has 
been recognized by many authors, but no measurements have previously 
been reported because of the difficulty of producing and handling crystals 
at low temperatures. We have carried out measurements of the velocity 
of longitudinal ultrasonic waves in solid argon leading to preliminary 
estimates of the compressibility at two temperatures. 

The freezing of liquid argon usually produces a loose heap of small 
crystals (de Smedt and Keesom 1925), but we have found that slow 
deposition from the vapour at temperatures below the triple point 
(83-8°K) results in a compact, transparent, polycrystalline solid. For 
the measurement of the ultrasonic velocity it was possible to deposit 
the solid directly on to one face of a gold-plated, ‘ X-cut ’, quartz crystal, 
whose other face was in contact with a copper rod maintained at the 
desired temperature of the experiment. When the crystal had been 
covered by a layer of several millimetres of solid argon, a metal piston 
whose surface was parallel to that of the quartz crystal could be pressed 
slowly into the argon to produce a flat reflecting surface for the 
longitudinal ultrasonic waves—the excess solid being caused to. 
vapourize by maintaining the piston at a slightly higher temperature. 

The quartz crystal was excited by a 3 Mc/s sine-wave generator, 
producing standing waves in the solid argon which could be recovered 
each time the piston was moved (by a micrometer screw) a half wavelength 
further into the solid. The following preliminary values were obtained 
for the velocity (v) of longitudinal waves in an infinite solid : 1600 m sec-1 
at 60°K ; 1300 msec"! at 78°k. Results previously obtained with the same 
apparatus for the liquids ether and carbon tetrachloride at room 
temperature agreed with the published values to within 1°%. Because 
of severe attenuation, however, the indications of resonance in the 
polycrystalline argon were relatively weak, and the accuracy of these 
results cannot be so high. Attenuation at the grain boundaries probably 
accounts also for the complete failure of our first attempts, carried out 
at 15 Me/s using a pulse technique, because at this frequency the 
wavelength is about 0-1 mm—of the same order as the grain size of the 
solid argon as estimated from x-ray diffraction photographs. 
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Since we have not yet measured the velocity of shear waves, we can 
give an estimate of the compressibility («) only if we assume a value for 
Poisson’s ratio (c), Assuming c=} and the density (4) of solid argon 
to be 1-63 g cm~3 (Simon and von Simson 1924, de Smedt and Keesom 
1925, Clusius and Weigand 1940), we estimate from 


v=/{3(1—o)/[(1 +0) ]} 

that « is 0-36 x 10-1° cm? dyne~! at 60°K and 0-54 10-1° cm? dyne? at 
78°K. The chief causes of uncertainty in these values are: (a) the 
necessity at present to assume a value of Poisson’s ratio (a change of o 
from 4 to 0-45 would decrease x by 25°); (b) the possible existence of 
up to 0-5% impurity in the argon; (c) an uncertainty of perhaps 2% 
in the assumed value of density ; (d) the possible existence of a small 
temperature gradient in the argon; (e) an effect of dispersion due to 
scattering at the grain boundaries. : 

Herzfeld and Goeppert-Mayer (1934), using an interatomic potential 
of the form Ar-*+Be-"?, have set up an equation of state for solid 
argon, which has been re-calculated by Kane (1939) with the use of 
improved empirical data. We have evaluated the compressibility for 
our temperatures from Kane’s p—V isothermals for two values of p. 
Another estimate may be obtained from the work of Rice (1944, 1951), 
who calculated the potential energy of the lattice for various values 
of the lattice parameter by subtracting the Debye thermal and zero-point 
energies from the experimentally determined total energy (we have added 
a small correction to « to take account of the thermal pressure). These 
theoretical estimates are listed in the following table for comparison with 
the experimental results :— 


«X 10'° (cm? dyne~*) 


T' (°K) After Kane (1939) 
Experimental |_——-————_—————_|_ After Rice 
(p=0-2091) (p=0-345) (1951) 
60 0-36 0-50 0-85 0-15 
78 0-54 T 7 1-6 0-21 


Although, strictly, the experiments yield values of the adiabatic 
compressibility, the differences between the isothermal and adiabatic 
compressibilities are much smaller than the experimental uncertainty 
and we have not distinguished between them. It will be seen that the 
experimental values lie in the range covered by these theoretical estimates. 
Because of inconsistencies in the data used by the authors referred to, 
and because of the present experimental uncertainty, we do not wish 
at present to comment in greater detail on the agreement between the 
experimental and theoretical estimates of compressibility. It is interesting 
to note, however, that the observed compressibility of argon is higher 
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than that of its neighbouring alkali metal, potassium, in agreement with 
the suggestion of Bridgman (1949) that the peaks in the curve of com- 
pressibility of the elements against their atomic numbers would be 
occupied by the solidified inert gases. We hope to carry out further 
measurements with refined techniques, using both pulse and resonant 
methods with shear and longitudinal waves. 


We wish particularly to thank Professor Kathleen Lonsdale, F.R.S., 
for her advice and for the help given in carrying out the x-ray measure- 
ments in her laboratory. We are grateful to Professor H. R. Robinson, 
F.R.S., for his friendly encouragement throughout the work, and to 
Dr. R. O. Davies for helpful discussions. We acknowledge with thanks 
the assistance of the Central Research Fund of the University of London 
and of the Department of Scientific and Industrial Research. 
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Determination of the Widths of the Levels of p*! from the Reaction 
S7AL (as pie sk 


By R. R. Roy and C. GopEAu 
Université Libre de Bruxelles, Laboratoire de Physique, Bruxelles* 


[Received July 28, 1953] 


THE reaction ?7Al (x, p) 2°Si has been investigated by many workers. It 
is known that the resonance reaction takes place corresponding to the 
incident energies of «-particles H,—4:0 Mev and 4-4 Mev which produce 
excitations of 14-3 Mev and 14-8 Mev in the compound nucleus p?!—above 
its ground state. At these excitations the compound nucleus decays with 
the emission of four groups of protons of Q-values 2-26 Mev, —0-02 Mev, 
— 1-32 Mev and —2-49 Mev. In the present experiment we have investi- 
gated the widths of the levels of p*! for energies of «-particles H,=4-0 Mev 
and 4:4 Mev, for the group Q=-—- 1-32 Mev. 

A diffusion chamber was used for the detection of proton tracks. The 
dimension of the chamber was 24 em in diameter and 15 em in height, and, 
filled with air, it was operated at atmospheric pressure, at a gradient of 


temperature of 7°c/em. The vapour source was a mixture of 70% 
a ee 
* Communicated by the Authors. 
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methyl alcohol and 30°% water. The cold plate was in contact with solid 
CO,. The aluminium target was bombarded by «-particles from 
polonium of strength 10 me. 

Photographs, 5 000 in all, were taken by a vertical camera through a glass 
window fixed on the top hot plate. The background drops could be 
reduced considerably simply by rubbing the glass window with a woollen 
cloth ; this produced an electric field across the chamber which removed 
the previous ions. The operation of rubbing the glass"plate was repeated 
every five minutes during which period 20 photographs could be taken. 
Photographs obtained under these conditions were very sharp. The 
accompanying plate (Plate 49) shows the track of a proton in which the 
phenomenon of * depletion ’ can clearly be seen. The depletion of vapour 
was caused by the earlier passage of a cosmic-ray particle across the 
proton track—a passage marked by very diffused ions. 


Fig. 1 


Number 


9.7 1073 109215 lea 
Range in ¢m. ———> 
The distribution of dispersion in ranges of the 107 proton tracks for 
E,,=4:0 Mev, Q=— 1-30 Mev. 

Considering the dimensions of the chamber, the group for Q=— 132 Mev 
for incident energies 4:0 Mev and 4:4 Mev could be most conveniently 
studied : the number of proton tracks measured was 107 and 85 respec- 
tively. The actual Q-values for the energies 4:0 Mev and 4-4 Mev were 
found to be — 1-30 Mev and —1-37 Mev. The proton tracks of a particular 
group, emitted at a definite angle with respect to the direction of the 
x-particles, should show a dispersion in their ranges due to the width of the 
level. Since the protons are of about the same energy, error owing to 
straggling of protons in the target is small. As the number of protons 
observed at any particular angle was not large, we have reprojected all 
the tracks at 0° with respect to the incident direction of «-particles. 

Figures 1 and 2 give the distribution of the ranges of protons for 
E,=4-0 Mev and 4-4 Mev respectively. The widths deduced from these 
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curves are 0-36--0-06mev for Q=—1-30 Mev, E,=4:0 Mev and 
j 7 for Q= A = 4-4 Mev. 

0-30-+0-06 Mev for Q=— 1-37 Mev, Ee | 

The two earlier experiments to determine the widths of the levels A 
the reaction 27Al (x, p) 2°Si, are due to Chadwick and Constable (1932) 
and Duncanson and Miller (1934). Their results are given along with ours. 
in the following table. 

Fig. 2 


(o) 


22 


Number 
fo) 


12,1 125 (aio SSS 
Range in cm. — => 


The distribution of dispersion in ranges of the 85 proton tracks for 
Ey,=4:4 Mev, Q=—1-37 Mev. 


Resonance levels 


Authors corresponding to Width of the level 
E. (Mev) 
Duneanson and Miller 5:8 0-3 +0-2 Mev 
Chadwick and Constable 5-2 0-12< r<0-32 Mev 
4:8 0-20 Mev 
Roy and Godeau 4-4 0:30 +0-06 Mev 
4-0) 0:36 +0-06 Mev 


The values obtained in the present experiment cannot strictly be: 
compared with theirs as we have confined ourselves to the protons of one 
group only, that is Q=—1-32 Mev, whereas their values refer to all the 
proton groups for the incident energies studied. 
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Temperature Change in Paramagnetic Resonance Absorption of Solutions of 
Organic Free Radicals 


By D. Bis and A. C. Rosz-Innes 
Clarendon Laboratory, Oxford* 


[Received August 10, 1953] 


We have investigated solutions (about 0:01 molar) of three organic 
radicals in different solvents (compare table 2) between room temperature 
and 90°K at about 10 000 Me/s and in magnetic fields of about 3000 oersted. 
The substances studied were : 

(1) Banfield and Kenyon’s radical (8-(phenyl nitric oxide)-8 methyl 
pentane-6-one oxime N-pheny] ether) 

ON (C,H;5) C (CH3), CH, C (CHs) N (C, H;) O ; 

(2) Di-phenyl nitric oxide 4, NO; | | 

(3) Di-p-anisyl nitric oxide (p-An), NO. 

At room temperature in sufficiently dilute solutions we always found 
three practically equidistant, not completely resolved, absorption lines. 
The separation, 4H, of the central line and one of the outer lines is listed 
in table 1. 


Table 1 
g (solution) 
g AH. 1 102A Av 
Substance (solid) Low Room _ | (oersted)|(cm~!)} (Mc/s) 
tempera- | tempera- 
ture ture 
is NO = 2-0066 2-0066 9-5 0-88 | 39-541 | 


B. and K.’s radical | 2-0052 | 2-0050 2-0052 11-0 1:03 | 46-2+1 


(An), NO 2-0072 | 2-0048 | 2-0037 | 10-4 | 0-97 | 44-6+1 


The accuracy of all g-values is about +-0°0005. 


Originally we hoped to improve the resolution of the absorption lines by 
cooling the solution, but at low temperatures the three lines merged into 4 
single line. This effect was studied in detail in solutions of B. and K.’s 
radical, but occurred as well in solutions of the other radicals. In table 2 
are listed the different solvents used and the temperature range in which, 
on cooling, the transitions occurred. The results were reproducible. 

Further details of the transitions are given in fig. 1 (Plate 50), where we 
reproduce some absorption lines observed as traces on the screen of a C.R.O. 
In one group of solvents (marked ‘ A ’ in fig. 1 and table 2)the pattern is 


* Communicated by Professor F. Simon, F.R.S. 
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symmetric at the highest and lowest temperatures, but is clearly asym- 
metric in between. In n-heptane the transition occurs below the melting 
point of the solvent. Inthe other group of solvents (marked ‘ B’) the pattern 
- is symmetric at all temperatures, and the transitions occur above the 
melting point of the solvent. Further analysis is made more difficult by 
the incomplete resolution of the lines. Assuming, however, that the line 
shape is the same for all lines at all temperatures, it was possible to re- 
construct the experimental traces by composition of three curves of 
different width and peak intensity. This procedure led us to the following 
qualitative conclusions. 


Table 2 
Tempera- 
Solvent Type ture Remarks 
range (°K) 
Methyl alcohol aa ee No change in H.F.S. 100° to 155°K 
lc Cl, HAN No change in H.F.S. 120° to 250°K 
ene ‘A’| 155-180 
ence “B’|} 150-185 
| CS, Be 180 Sudden change 
arsine “B’| 150-185 


Pet.-ether (b.p.<40°c)| ‘B’ | 141-155 


If the temperature is lowered in both cases ‘ A’ and ‘ B’, the width of 
all lines increases and, moreover, the central line becomes more intense 
relative to the outer lines. In case ‘ B’ the distance between the lines 
becomes smaller ; the same effect probably occurs in case ‘ A ’, 

The pattern observed at the higher temperatures strongly suggests a 
hyperfine structure caused by the 14N nucleus in the NO-group. This is 
corroborated by the hyperfine structure observed in solutions of a radical 
of the same class as the radicals we investigated, the nitroso disulfonate ion 
(SO3),NO~ ~ (Townsend et al. 1953). We therefore will interpret the 
observed lines as a hyperfine structure. This H.F.S. must be virtually 
isotropic and therefore can be described by the spin Hamiltonian 


H =gB (S.H)+A (8.1), 


where f is the Bohr-magneton, H is the magnetic field, gy and A are con- 
stants, and S and J are the spin operators of the unpaired electron and the 
nitrogen nucleus respectively (S224) f=). 

Table 1 contains the observed g-values (both for solution and solid), the 
separation, 4H, of the central line and one of the outer lines, the value of 
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A and the splitting in zero field, Jv. The values of g and 4H were 
independent of the solvent; the magnetic field was calibrated with 
a, «-diphenyl f-picryl hydrazyl (g=2-0036+0-0002). 

It is very unlikely that the change in relative intensity of the central 
line and the outer line during cooling of the solution can be explained in 
terms of a hyperfine structure. It is quite feasible, however, that on 
cooling one of the components of the system crystallizes out until a possible 
eutectic point is reached, below which both components would crystallize. 
This would probably make the observed pattern a superposition of a 
hyperfine structure of decreasing intensity and a single line of increasing 
intensity, coinciding (for B. and K.’s radical) with the central line of the 
hyperfine structure. 

The decrease in hyperfine structure is probably a more interesting effect. 
as it may throw some light on details of the interaction between. the 
unpaired electron and a nuclear magnetic moment in a free radical. 
Further experiments on this phenomenon are in course of preparation. 


We would like to thank Mr. R. Foster, B.A., Dyson Perrins Laboratory, 
Oxford, and Dr. A. Burawoy, Chemistry Department, College of Techno- 
logy, Manchester, for providing us with samples. One of us (D.B.) is 
indebted for the award of a Pressed Steel Company, Ltd. Research 
Fellowship ; the other (A.C.R.-I.) for a maintenance grant of the Depart- 
ment of Scientific and Industrial Research. 


REFERENCE 
TOWNSEND, J., WEISSMANN, S. I., and Paks, G. E., 1953, Phys. Rev., 89. 606. 


On the Reaction 18O(p, «) &N 


By R. R. Roy, A. Lagasse and Miss M. J. Decock 
Université Libre de Bruxelles, Laboratoire de Physique, Bruxelles 


[Received July 31, 1953] 


THE reaction 180 (p, «) &N was established by Burcham and Smith (1939) 
and by Freeman (1950); since then, the resonance reaction has been 
reported by Mileikowsky and Pauli (1950) and Seed (1951). The first 
resonance, according to Mileikowsky and Pauli, occurs for energy of proton 
E,,=680 kev and according to Seed at 640 kev: the latter value is now 
generally accepted. At this incident energy the compound nucleus AY 
receives an excitation of 8-59 Mev above its ground state, and then decays 
with the emission of «-particles. The nucleus !N is similar to the nuclei 
of 7Li, 4B and 1°F in that it can be regarded as composed of «-particles 
minus one nucleon. The lowest levels reported so far for 7Li, 11B and ¥F 
are at 0-478 Mev, 2-14 Mev and 1-4 ev respectively ; on the other hand, 
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the first low level observed in 45N is at 5-28 Mev from the reaction 
14N (d, p) &N (Hornyak et al. 1950). 

We have investigated the possible existence of low lying excited levels 
in N using the reaction 18O (p, «) &N with energies of protons up to 
700 kev. 

The 4-7°%, enriched 18O target was in the form of water which was 
deposited as ice on a copper disc inside an evacuated glass cylinder of 
pyrex, while the proton beam, the current of which was measured by an 
integrator connected with the target assembly, bombarded the target after 
it had been deflected by the magnetic analyser. The voltage scale of the 
generator was calibrated in accordance with the standard reaction 
7Li (p, «) Be. For the detection of the y-rays which might be emitted 
from an excited level of 1°N in the process of de-excitation, a scintillation 
counter unit with a NaI-thallium enriched crystal was employed. 


—-~ Yo wo Ff oO HAN @O O 


400 65 700 


Excitation curve of the reaction !80 (p, «)!®5N for energies of proton up 
to 700 kev. 


The figure shows the result obtained: the ordinate N represents the 
number of counts per x-Coulomb and the abscissa the energy of the 
proton. Between 600 kev and 700 kev the number of counts, which varies 
linearly with the incident energy of protons, was noted at intervals of 
10 kev. The width of the resonance level, as can be inferred from the 
curve published by Seed, is about 45 kev, corresponding to the proton 
energy 640 kev. If there had been an excited level in °N for the incident 
energy 640 kev, we should have expected a sudden increase in the number 
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of counts, starting from about 618 kev. The absence of such an increase 
leads us to conclude that, in the first place, no excited level is formed in 
1N for the reaction 180 (p, «) ™N at the resonance energy of proton 
640 kev and, secondly, that the compound nucleus °F decays with the 
emission of one group of «-particles, leaving the final nucleus 1=N in the 
ground state. It should be mentioned that, with an ice target from pure 
distilled water, a curve similar to that in the figure was obtained. 

Our results might be interpreted in the following way. Protons falling 
upon the ice target produced x-radiations ; the higher the energy, the 
greater the depth of penetration of the protons into the target which gives 
rise to the linear increase in the number of counts. 
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The Cross Section for 1°1Ta (y, n) 1%Ta at 17-6 Mev 


By J. H. Carver and H. J. Hay 


Research School of Physical Sciences, Australian National University, 
Canberra* 


[Received July 28, 1953] 


THE apparent ‘ resonance ’ character of many photo-nuclear reactions is 
now well known. In the main these processes have been studied with the 
aid of the bremsstrahlung produced by betatrons and synchrotrons and the 
reaction identified by activation techniques (a number of references has 
been given by Carver, Edge and Wilkinson 1953). One of the reactions 
which has been most thoroughly investigated by this technique is the 
181Ta (y, n) 18°Ta reaction. The figure shows the cross section for this 
reaction as determined by Haslam, Smith and Taylor (1951). Using the 
gamma-rays from the Li (p, y) reaction, Carver, Edge and Wilkinson 
(1953) showed that competing reactions did not destroy the ‘ resonance ” 
character of the gamma-ray absorption. They found the ratio of the 
cross sections for 181Ta (y, 2n) 17°Ta at 17-6 Mev to 181Ta (y, n) 18°Ta at 
14-6 mev to be 0:29+0-11.f Taken at their face value these measurements 
indicate that, at gamma-ray energies of the order of 18 Mev and greater, 


* Communicated by Professor E. W. Titterton. 

+ This result agrees with the neutron yield measurements obtained by 
Halpern, Nathans and Mann (1952) and Whalin and Hanson (1953) using 
betatron bremsstrahlung and with the indirect estimate made by Eyges (1952). 
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the gamma-ray absorption occurs almost entirely by means of direct 
interaction rather than by the immediate formation of a compound 
nucleus. The (y, 2n) cross section at 17-6 Mev (Carver, Edge and 
Wilkinson 1953) is only slightly larger than the (y, n) cross section at 
17-6 mev (Haslam, Smith and Taylor 1951) whereas calculations based on 
the statistical theory (Blatt and Weisskopf 1952) indicate that the (y, 2n) 
cross section at 17-6 Mev should be some 16 times the (y, n) cross section 
at 17-6 Mev. 

In order to determine whether or not the high energy tail in the 
181T'a (y, n) 18°Ta cross section is indeed as large as is indicated by the 
bremsstrahlung results we have measured the ratio of the cross sections 
for this reaction at gamma-ray energies of 14-6 and 17-6 Mev using the 


12 


lO 


CROSS SECTION ARBITRARY UNITS 


8 10 12 14 16 18.) 20) Ff22 
GAMMA-RAY ENERGY MEV 


The cross section for 18!Ta (y, n) '8°Ta at 17-6 Mev compared with the curve 
obtained by Haslam, Smith and Taylor (1951). 


radiation from the 7Li (p, y) reaction. The principle of the method was 
the same as that used by Carver, Edge and Wilkinson (1953) except that a 
measurement of the induced 1°Ta activity (half life 8-15 hr) replaced the 
total neutron yield measurement. 

The proton bombardment of lithium produces gamma-rays of 17-6 Mev 
(negligible width) and 14-6 Mev (2 Mev width) and the relative intensity of 
these gamma-rays depends on the proton energy. At a proton energy of 
500 kev (thick target) the ratio of the intensity of the 17-6 Mev line to that 
of the 14-6 MeV line is 1-70+-0-20 ; at 1-15 Mev (target thickness 250 kev) 
this ratio Is 0-45+-0-06 (Stearns and McDaniel 1951). These ratios apply 
to radiation emitted at an angle of 75° to the proton beam. We have 
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irradiated a tantalum Geiger counter under these conditions and deter- 
mined, at both proton energies, the relative induced 1°Ta activity per 
incident quantum. Throughout the irradiations (which were for periods 
of about eight hours) the gamma-ray yield was monitored continuously 
by a Nal (TI) crystal which was biased so that none of the ~0-5 Mev 
quanta, due to inelastic scattering at the higher proton energy, was 
detected. Several irradiations were performed at both proton energies 
and corrections were made for fluctuations in the gamma-ray yield 
assuming the half-life of 18°Ta to be 8-15 hr. It was found that the ratio 
of the 1°!Ta (y, n) 18Ta cross section at 17-6 Mev to that at 14:6 Mev 
was 0-26+0-06.* As shown in the figure this result agrees with the 
bremsstrahlung results and therefore confirms the conclusions drawn by 
Carver, Edge and Wilkinson (1953) concerning the importance of direct 
effects at gamma-ray energies above the ‘ resonance ’. 
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Introduction of Exchange in the Thomas—Fermi Model for Atoms 


By N. H. Marcu 
Department of Physics, The University, Sheffield 


[Received August 4, 1953] 


In a recent paper, Scott (1952) has considered the corrections to be 
applied to the customary formula due to Milne (1927) for the binding 
energy of the Thomas—Fermi (T.F.) atom, and has obtained results in 
good agreement with experiment. Scott concluded that two corrections 
were necessary, the first arising from the so-called ‘ boundary effect ’, due 
to the failure of the T.F. method in the neighbourhood of the nucleus, 
and the second from the effect of exchange. 


* This result applies to a monochromatic gamma-ray of 14-6 Mev and not to 
the broad line used in the experiment. We have corrected our result for the 
width of this line (and for the 0-6 Mev increase in gamma-ray energies at the 
higher proton energy) by assuming that near 14-6 Mev and near 17:6 Mev the 
shape of the (y, n) cross section is as given by Haslam, Smith and Taylor (1951) 
but no assumption as to the relative magnitudes of the 14-6 and 17-6 Mev cross 
sections has been made. 

+ Communicated by the Author. 
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The present work has been carried out with the object of examining 
the validity of Scott’s method of evaluating the binding energy (without 
the boundary correction, with which we shall not be concerned _ here) 
when exchange is included and investigating the possibility of modifying 
the usual T.F. electron distribution to take some account of exchange. 
For it is well known that if one attempts to use Dirac’s method directly to 
calculate a modified electron distribution for an atom (Dirac 1930) 
serious difficulties arise, as there is no solution of the Thomas—Fermi— 
Dirac (T.F.D.) equation representing an electron cloud which extends to 
infinity. In fact the only kind of ‘atom’ in the model is one of finite 
radius with a discontinuity of charge density at the boundary and this is 
clearly very unsatisfactory.* Scott avoided the difficulties associated 
with the T.F.D. model by using the ordinary T.F. electron density to 
calculate the energy including exchange. 

Scott’s choice suggested to us the use of the variational function 


n(7) =A n (Ar). 22 ee 


for the electron density, n,(r) being the usual T.F. density. A° is included 
in (1) simply for normalization purposes. Such a scale factor X was first 
introduced by Fock (1932) as a device to enable him to prove the virial 
theorem in the T.F. theory, and was subsequently put equal to unity. 
We propose here to use A as a proper variational parameter. If H, is the 
total energy corresponding to the electron density (1), and 7',, U, and A, 
are respectively the kinetic, potential and exchange energies calculated 
using the unmodified density n,(7) then it follows that 


#H,=T,4 U,+4, ] 
Pas aie a) ene). e 
=AT,+AU,+AA, | (2) 


Hence the value of A, say A, for which #, is a minimum is given by 


(U,+A)) 

A= =H 27a . . . . . . . . . (3) 

Thus by choosing A to have the value A given by (3) instead of 1 as in 

Scott’s calculation we find a lower energy and hence an improved charge 
distribution. 

If E,(=T',,+U,) is the usual T.F. energy neglecting exchange, we can 

write, from the work of Milne (1927), Baker (1930) and Scott (1952) 


E,=—fZ'8, U,=—2pZ78, 


5 


T =p, A,=—yZ58, (4) 
where B=0-7687 and y=0-2208, both in atomic units. 
ee ee ee 

* It should be mentioned here that Plaskett (1953) has recently criticized 


Dirac’s treatment and proposed an alternative equation. This modified 
equation has associated with it all the difficulties of the T.F.D. treatment. 
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Hence we obtain from (3) and (4) 


A=1+(y/26)Z-2° , 
=1+0-1436Z-2/3 (5) 
Equations (1) and (5) completely define a new charge distribution in 
which some account is taken of exchange. 
It is instructive at this stage to write down the separate expressions for 
the kinetic, potential and exchange energies. They are easily found to 
be given in terms of 8, y and Z by the following equations : 


T 4 =BZ"P + yZ°8+ (P/48)Z, 


U ,=—2BZ78—yZ52, (6) 


A = —yB— (7 /28)Z. j 
It follows that the energy #, in this approximation is given by 
Hi PAF —yZ 8 —(y2/48)Z, 2 «ws. (7) 


which agrees with Scott’s result apart from the term in Z. Since he 
neglects terms of order Z*? it follows that his binding energy formula is 
left unchanged. 


20 


D(r) in atomic units. 
no FF D © 
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O-15n 6-28 0560-4) 05 606-0-7,,08%095° 1-0 
r in atomic units. 


Radial charge distribution for argon. 
(1) Without exchange. (2) With exchange. 


We see then that whilst the effect of modifying the electron density to 
take account of exchange is to introduce terms of order Z° in the kinetic 
and potential energies, these terms exactly cancel in the total energy and 


1196 Correspondence 
Scott’s result is effectively regained. It is further interesting to note 
that from (6) 
2T ,+U,+A,=0, .% apes Bee Pos 
and hence the virial theorem is satisfied in our approximation, whereas 
it does not, of course, hold in Scott’s treatment. 

To illustrate the effect of modifying the electron distribution in the 
manner considered here, we have plotted in the figure the radial charge 
distribution D(r) (defined as density x 47r?) for argon, with and without 
exchange. The contraction of the charge cloud which is known to occur 
when exchange is introduced is well shown. (For larger values of r than 
are shown in the figure, curve 2 of course lies below curve 1.) 

The method outlined above can also be applied to positive ions, for which 
a direct application of the T.F.D. equation again leads to difficulties. 
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Event No. 5 (table 2). Example of anomalous scattering. 
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Probable example of anomalous scattering. 
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Slip lines on (111) face of ZnS, 
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A typical example of a proton track obtained in the diffusion chamber, showing 
the phenomenon of ‘ depletion ’. 
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Paramagnetic resonance absorption in solutions of B. and K.’s radical 
at various temperatures. 
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